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colliders 
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This thesis is concerned with the theory and the phenomenology of rapidity gap 
processes. We perform perturbative calculations of energy flow observables in jet-gap- 
jet processes, which consist of resummed primary emission calculations specific to the 
soft gluon geometry at HERA and an estimate of non-global (secondary emission) 
effects in clustered energy flow observables. The resulting predictions agree well with 
HI data. We also study hard diffraction and use a factorised model, with a Monte Carlo 
event generator, to make detailed predictions for gap-jet-gap events at the Tevatron. 
We find that we can describe the data in a natural way by using HERA parton densities 
and a gap surivial factor consistent with theoretical estimates. 
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Chapter 1 
Introduction 



In this thesis we will use the theory of the strong force, Quantum Chromodynam- 
ics or QCD, to calculate detailed predictions for rapidity gap processes and perform 
comparisons to experimental observation. 

Chapter 2 is devoted to a survey of the fundamental ideas of QCD, and equips 
the reader with some of the tools that are used in the rest of the thesis. We start 
with an introduction to QCD, explain the crucial idea of asymptotic freedom, and 
proceed with a discussion of factorisation. The colour mixing matrices of QCD are 
then described, with an example calculation for a quark process and a gluon process, 
and then we outline the basic ideas of Regge theory, in preparation for our later 
studies of diffractive processes. We continue with brief sections on Monte Carlo event 
generators and rapidity gaps, and finish with a summary. 

In chapter 3 we study the diffraction of hadrons at the Tevatron. In these pro- 
cesses, diffracting hadrons produce a central dijet system which is separated from the 
intact hadrons by rapidity gaps. These gaps are attributed to pomeron exchange. 
We describe the factorised model of Ingelman and Schlein, which views the process 
by double pomeron exchange, and use the diffractive Monte Carlo event generator 
POMWIG, coupled with pomeron parton densities from HERA, to produce a set of 
diffractive predictions. We test these predictions against Tevatron data and find, by 
using a gap survival factor consistent with theoretical estimates, we can naturally 
describe the experimental observations. 



Chapter 1. Introduction 



19 



Chapter 4 is a survey chapter examining resummation, which is a consequence 
of the factorisation properties obeyed by QCD cross sections, and applies the idea 
to rapidity gap processes. We start by writing down the factorisation properties of 
the cross section in a specific region of phase space, using jet functions to describe 
approximately coUinear quanta and a soft function to describe soft gluon emission, and 
develop the resummation formalism. The result is an expression for the cross section 
in which the large logarithms in the jet and soft functions are resummed. We develop 
the application of resummation to rapidity gap processes, focusing our attention on the 
soft function, which describes the soft, wide angle emission of gluons into a restricted 
region of phase space. 

The theme of chapter 5 is non-global logarithms in interjet energy flow observ- 
ables. The resummation formalism to describe interjet energy flow, discussed in chap- 
ter 4, fails to include the effects of secondary gluons, which are radiated outside of the 
rapidity gap and subsequently radiate into it. We study the effect of these secondary 
gluons on 2 jet cross sections in the presence of a clustering algorithm, at leading order 
and at all-orders, and make predictions for the impact of non-global logarithms based 
on an overall, gap dependent suppressive constant. We find that, compared to the 
non-clustered case, the use of a clustering algorithm reduces, but does not remove, the 
suppressive effect. 

In Chapter 6 we draw on the analyses of the previous chapters and apply our 
ideas of the resummation of interjet energy flow and non- global observables to gaps- 
between-jets measurements at HERA. We include primary interjet logarithms using 
the resummation formalism of chapter 4 by making detailed soft gluon calculations for 
the specific gap geometry at HERA. Non- global logarithms are approximately included 
by an extension of the work in chapter 5. We find that our calculations are consistent 
with HI data, and we make predictions for the ZEUS gaps-between-jets analysis. 

Finally, in chapter 7 we draw our conclusions. 



Chapter 2 

QCD at the frontier 



2.1 The theory of the strong force 

We are concerned with the theory of the strong force - Quantum Chromodynamics, or 
QCD. This theory attempts to describe the fundamental constituents of hadrons using 
point-hke quarks and gluons; the former making the matter content of the hadrons, 
and the latter mediating the colour force which binds the quarks together. 

Baryons and mesons were suggested to have a composite nature in the early 
sixties, resulting in the colour degree of freedom being introduced to maintain the 
fundamental link between spin and statistics, and providing the baryons with an anti- 
symmetric wavefunction. Feynman then continued the theme of hadron constituents 
with his high energy parton model, which he used to explain scaling properties of DIS 
and Drell-Yan cross sections. Today the theory of quarks and gluons, enjoying the 
status of a Yang-Mills non-abelian gauge theory, is a cornerstone of the "standard" 
model of particle physics. There are six quarks, each of which is an SU(3) triplet, 
interacting through the gluons, which are the SU(3) gauge bosons. 

We shall now give an overview of some of the tools of QCD we will be using in 
the rest of this thesis; we shall resist the impulse to list the QCD Lagrangian and the 
other trappings of modern gauge theories and assume the reader is familiar with the 
basic ideas. For further details of this, and the rest of the material in this chapter 
see [1-4] . In all of this work, we will use dimensional regularisation with d — A — 2e 



Chapter 2. QCD at the frontier 



21 



and use the Feynman gauge, unless otherwise stated. Some of the figures in this thesis 
were produced using the package axodraw [5]. 

The main property that allows a perturbative approach to QCD is the feature of 
asymptotic freedom. This means that the coupling constant, as, is a function of the 
scale Q and decreases as the scale increases. For a one-scale problem we write our 
(dimensionless) observable as a perturbative expansion in ag, 

n \ r* / 

where /i is the scale used to renormalise the theory. We now note that the observ- 
able a cannot possibly depend on the precise choice of /i, for we can (in principle) 
measure a in our laboratory and a change in a theoretical scale cannot affect such a 
measurement. Hence 

This equation, and many like it, are a vital part of a particle physicist's toolkit and 
will play a central role in the analysis of this thesis. The solution of this equation 
indicates that the scale dependence, or running, of ag is given by the solution of the 
so-called renormahsation group equation (RGE), 

Q'^,^P{as), (2.3) 

where the QCD /^-function is calculable using QCD. By restricting ourself to a one-loop 
solution we find. 



where P{cks) — —ba^ and^ 



33 — 2nf . ^. 

iZTT 



We have written this equation in terms of the experimentally determined parame- 
ter Aqcd, at which the coupling diverges and perturbative calculation breaks down. 
The most commonly used value of Aqcd is the five- flavour QCD scale Aq^^'^, where 
^Uf is the number of flavours. 
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we calculate the contributing one-loop Feynman diagrams in the modified minimal 
subtraction scheme. The world average value for the coupling at the mass is 

a^(Mz) = 0.118 ±0.002, (2.6) 

for which A^^'^'^ is deduced to be 

Agg^ = 20811 MeV. (2.7) 

The non-abelian gauge group of QCD ensures a negative /3-function (assuming n/ < 16) 
and hence a reduction^ of the coupling with scale. The strategy of perturbation the- 
ory is to expand the observable in powers of a^, equation (2.1), and hope that, by 
the smallness of q;^, it is sufficient to calculate just the first one or two terms of the 
expansion. 

In general, the perturbative expansions of QCD two-scale observables, generally 
denoted i?, are littered by large logarithmic enhancements of the form 

n m 

where V is the ratio of the two scales (normally the hard scale and a softer scale) . The 
leading logarithmic (LL) set is the set of terms with the most number of logarithms 
for a given cts; note that the definition of the LL set is observable dependent, and 
can be up to two logarithms per a^. If these large logarithms overcome the smallness 
of ttg, then it is insufficient to calculate the first one or two terms in the perturbation 
series, because all the terms are potentially large, and all orders must be considered. 

The physical origin of these large logarithms is the soft and/or coUinear limit of 
Feynman diagrams. In these regions of phase space, the denominators of some internal 
propagators vanish and these regions are logarithmically enhanced. Examples of ob- 
servables which include large logarithms include the thrust (T) distribution as T — > 1 
in electron-positron cross sections and cross sections in the vicinity of partonic thresh- 
old, where the partonic system has just enough energy to produce the observed 
final state of mass Q. In the latter example, large logarithms occur in the limit 2; — > 1, 

^QED, conversely, has an abelian gauge group and a positive /3-function. The result is 
that weakly rises with scale. 
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where z — Q"^ js. The theoretical machinery to include the effect of terms of all-orders 
is known as rcsummation. In this thesis we will be primarily concerned (beginning 
with chapter 4, once we have conchided our study of diffractive processes in chapter 3) 
with the effect of large logarithms which arise from restricting soft particle emission 
in restricted regions of phase space, where the large terms arise from an incomplete 
real and virtual diagram cancellation. An important tool for such an analysis is the 
use of quantum mechanical incoherence, to which we now turn. 



Factorisation is a statement of the quantum mechanical incoherence of short and long 
distance physics, and plays a central role in the predictive power of QCD. In this section 
we will outline the statements of factorisation that are of most use to us, namely the 
"standard" factorisation theorems, which write cross sections as convolutions of long 
and short distance functions, and the refactorisation theorems of the short distance 
function. 

Factorisation is the QCD generalisation of Feynman's parton model. Colliding 
hadrons, in the centre-of-mass frame, are highly Lorentz contracted and time-dilated 
and the interaction probes a frozen configuration of partons. The interaction thus 
proceeds by one parton from each hadron undergoing a QCD hard scattering event; 
the spectator partons cannot interfere with this process as the interactions between 
these take place at longer, time-dilated scales. The unscattered partons go on to form 
the hadron remnants. The inclusive cross section is written as a convolution of a 
long distance function describing the dynamics of partons in the hadrons with a short 
distance function describing the hard event. 



where is a factorisation scale separating the long distance dynamics from the short 
distance dynamics, fi/h^{x,fj,) denotes the distribution of parton i in hadron A with 
momentum fraction x (known as a parton distribution function, a parton density or 



2.2 Factorisation and refactorisation 
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a PDF) and Q is the scale of the process. The parton densities are non-perturbative 
and arc required to be determined experimentally, whilst the short distance function 
is calculable using perturbation theory. 

A remarkable consequence of factorisation is that measuring a parton density 
at one scale /i allows us to predict the parton density at another scale /i', provided 
that /X, ^ A. This result, known as the evolution of parton densities and structure 
functions, is a powerful predictive tool in perturbative QCD (pQCD). This evolution 
is most transparently expressed using a set of integro-differential equations, 

^^'^J^/kA^,^^')= E fd^P^^{J^<^s{n')^fJ/HAtn'), (2.10) 

which are known as the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions and are one of the most important sets of equations in pQCD. The evolution 
kernels (or splitting functions) give the probability of finding species i in species j, 
and are calculable as a power series in a^. The combination of factorisation of ob- 
servables into short distance functions and non-perturbative parton densities, and the 
subsequent evolution of the parton densities using the DGLAP equations is central in 
the impressive success of QCD as the gauge theory of the strong force. 

We can now perform a further refactorisation on the short distance function, for a 
specific class of observables in a particular limit of their final state phase space. In this 
region, we are interested in a QCD hard process at scale Q, which is only accompanied 
by soft radiation up to the soft scale Qs- This region of phase space is known as the 
threshold region, and is relevant if we make a specific restriction on the energy of gluon 
emission or, for example, in the production of heavy quarks near threshold. The soft 
radiation is described by a function S and we write the cross section as the product 
of a hard and a soft matrix, 

o^^'^iQ, Qs, as) = E HiLiQ, asifif))SuiQs, «.(/^/)), (2.11) 

I,L 

where /x/ is a new factorisation scale. The proof of this statement follows standard 
factorisation arguments [6,7]. The soft gluon emission is sensitive to the colour state 
of the hard event and hence we have written the soft and hard functions as matrices 
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in the space of possible colour flows of the process. Now all the dynamics at the softer 
scale Qs arc described by the soft matrix, with all higher energy dynamics described by 
the hard matrix. For example, we can apply this factorisation to a dijet process, and 
restrict the interjet radiation to the scale Qs] in these rapidity gap processes, the soft 
function describes soft radiation into the interjet gap. We will use this factorisation in 
chapter 4 and in chapter 6, where we will apply these ideas to rapidity gap processes 
and exploit the factorisation to resum large interjet logarithms in energy flow processes 
at HERA. 

2.3 Colour mixing in QCD 

The refactorisation properties discussed in the last section will be exploited in later 
chapters to resum large QCD logarithms. An important tool in these calculations is 
the mixing of the basis of colour tensors^, over which the hard and soft matrices of 
the previous section are expressed, by quantum corrections. In this section we will 
calculate these mixing matrices for a quark process and a gluon process at one loop. 
The full set of mixing matrices for all processes is in appendix F and appear in [8,9]. 

The physical importance of the decomposition of an observable into its possible 
colour flows becomes clear when we note that the emission of soft radiation in the 
QCD process is sensitive to the colour state of the hard interaction. This colour 
coherence effect means that the soft gluon emission pattern depends on the overall 
colour charge of the parent system. This dependence of the radiation on the colour 
state can be understood by consideration of a QCD quark-antiquark scattering process, 
for example, in the large Nc limit. If the quarks interact by exchanging colour, then 
the outgoing quarks will be colour connected, and a colour dipole will be stretched 
between the outgoing partons. Hence the region between the quarks will be filled with 
gluonic radiation. However, if the outgoing quarks are not colour connected then there 
will be no dipole stretched between them. Therefore the soft gluon radiation pattern 
is sensitive to the colour state of the hard scattering, and we are required to consider 
^The full set of colour bases (or tensors) used in this work is in appendix D. 
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the decomposition of an observable into its different colour flows. 

We shall consider the colour flow of QCD 2—^2 scattering, which is described 
by 4 colour indices. The initial state particles will be labelled A and B, the final state 
particles will be labelled 1 and 2 and we will use lower case roman indices for internal 
lines. In this notation, the standard colour algebra for qq scattering with a t-channel 
gluon would then be written as ti^tg2- In order to decompose the colour flow of a 
matrix element, we need to specify a basis of colour tensors, linking the 4 indices, which 
describe the possible underlying colour flows. For example, the process qq qq will 
have a two-element basis consisting of elements which have the physical interpretation 
of singlet or octet colour exchange exchange, 

Cl = SaiSb2, 

1 1 

C2 = --^SaiSb2 + ^SabSi2, (2-12) 

where we denote elements of the basis as q. Note that C2 in this basis is interpreted 
as the colour flow for a t-channel gluon. The 2 — > 2 scattering amplitude can then be 
decomposed over this basis, 

M = (2.13) 

i 

where the Aii coefficients encode the amount of colour tensor q in Ai. The calcula- 
tion of these coefficients will allow the successful decomposition of an arbitrary QCD 
amplitude over an appropriate basis. 

The colour tensor basis set will be mixed into itself by higher order diagrams. 
For example, a diagram with singlet colour flow will become a diagram with octet 
colour flow by the addition of a virtual t-channel gluon. To compute this effect, which 
is the aim of this section, we dress a colour tensor with a virtual gluon connecting 
two external legs and, by considering the colour content of the resulting diagram, 
express the result in terms of elements of the basis set. A diagram illustrating the 
addition of a virtual gluon to the tensor Cj is shown in diagram 2.1. Therefore the 
virtual gluons will cause the bare colour tensors to mix into each other, with a matrix 
describing this mixing; we refer to this matrix as the colour mixing matrix. Therefore 
if the undressed basis set is denoted (01,02)'^ and the basis set which has mixed under 
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quantum corrections is denoted (c/^jCy^, then 




(2.14) 



This is just standard operator mixing under quantum corrections and we produce a 
process and basis dependent matrix describing how the colour tensors mix. 

Note that although it is the dynamics which cause the colour tensor mixing, we 
are only interested in the resulting colour structure in this section. We will start with 
a detailed example for a quark-only process and then describe the complications in 
the presence of external gluons. Appendix A contains a set of SU(3) group identities, 
which are used in this section. 



2.3.1 Quark-only processes 

The colour mixing matrix for quark-only processes is found using the fundamental 
identity 

tpli = l^uSkj - ^JiM (2.15) 

where denotes a SU(3) matrix in the fundamental representation. For the process 
q{A) q{B) g(l)g(2) we choose the basis of equation (2.12), which encodes singlet 
and octet exchange in the t-channel. Note that any other choice that completely spans 
the colour space is acceptable; the benefit of a singlet-octet choice is that the basis is 
orthogonal, and the tensor product of different elements of the basis is zero. We can 
attach the virtual, soft gluon to any two of the external legs, which gives six possible 
attachments with three classes of colour structure. Let C^^^ represent the colour 
decomposition obtained from dressing the colour tensor with a soft gluon connecting 
the i and j external legs. Each class of diagram has an associated dynamical piece, but 
we are only interested in the colour structure at this stage. It is these dynamical pieces 
that will ultimately undergo the mixing, through their associated colour structure. We 
have illustrated the dressing of a colour tensor in figure 2.1, where we showC,^^^^ We 
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Figure 2.1: The diagram obtained by dressing the colour tensor q with a virtual gluon 
a which connects external legs A and B, for the process ^'(^4) q{B) — > g(l) q{2). In the 
notation in the text, this diagram is denoted C\^^\ Note that in this section we are 
only concerned with the colour flow, and not the dynamics, of this diagram. 



can now work out the colour decomposition. Starting with cj"^^^ we get the following. 



(AB) 



^Bn ^mAi 



(2.16) 



where m and n are dummy indices. Note that we have suppressed the colour indices on 
the colour tensors Cj, in accordance with equation (2.12). By contracting the indices 
with the 5-functions we find that cj^"^^ 
similar manipulations we find that 



C2, which also means that C[^^^ = C2- By 



^(Al) 



(2.17) 
(2.18) 



For the first of these results we have used the fact that 



(2.19) 



f 12) 

Now we turn to the cases of diagrams involving C2. Starting with Q we have 



C. 



(12) 



+0, ±a 
^Bn t-mA) 

^'bu 2 {^mn^l2 — J^^ml^n2 



mA' 



(2.20) 



r» o 1 
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If we expand the terms, then we can use equation (2.19) on the first term and equa- 
tion (2.15) on the second to obtain 

Cf = ^c.+ fc.--L)c.. (2.21) 



2Nc \ 2N,^ 
A similar calculation for the other classes of diagram yields, 

er - Cr^-i-^c^^^c.. (2.22) 

We can now write down the mixing matrix for this process in terms of the three classes 
of diagram, and we use a, (5 and 7 to denote the associated dynamical piece, C,^'^^\ of 
the relevant diagram, 

^ = ^(Bi) + ^(A2)_ ^2.23) 

We will use this notation, and discuss the origin of the dynamical pieces, in chapter 6. 
Recall that if the undressed set is denoted (ci, 02)-^ and the set which has mixed under 
quantum corrections is denoted {€[,€'2)'^, then 

(2.24) 









U J 







(2.25) 



where the mixing matrix, (7^^^*^*, is given by 

- - I Cf(3 ^(« + 7) 

ya + T CfQ;- 2|^(q; + /3 + 27) 

This matrix describes how the colour tensors, and the associated dynamical pieces, mix 
under quantum corrections. The mixing matrix for the other quark process, qq — > gg, 
can be calculated in a similar way and appears in appendix F. 

2.3.2 The addition of gluons 

The presence of external gluons increases the complexity of the colour mixing cal- 
culations, due to more involved group theory. In this section we will describe the 
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calculation of the mixing matrix for the process q{A) g{B) — > q{l)g{2), which will 
introduce the majority of the tools needed for the gluon processes. The basis for this 
process is 

Ci = 6aiSb2, 

C2 — 0?B2c(^f)iA' 

C3 = ifB2c{t'F)lA, (2.26) 

where the first element describes singlet exchange, and the second and third mem- 
bers describe symmetric and antisymmetric colour exchange respectively. The con- 
stants f^^'^ (^A, B,C = 1, . . . ,8) are the antisymmetric structure constants of SU(3). 
They are antisymmetric under the exchange of any two indices and satisfy the Jacobi 
identity, 

/aSe/eCD + fcBsfAED + foBEfACE = 0. (2.27) 

The constants d^^*^ (^4, B,C — 1, . . . ,8) are the symmetric structure constants of SU(3). 
They are symmetric under the exchange of any two indices, obey (Iabb — and satisfy 

fABEdECD + fcBEdAED + foBEdACE = 0. (2.28) 

From these, and the group theory identities in appendix A, from which we mainly use 

(2.29) 

we can deduce the following basis relationships. 



.a.b _ 1 



J^^ab^ik + {dabc + ifabc) tik 



C3 X Ci - C3, (2.30) 
C2 X Ci = C2, (2.31) 

C3XC2 = -T^2 + ^^C3, (2.32) 
ca X C3 = ici + - ^C3, (2.33) 

which are useful when we see that adding a virtual t-channel soft gluon to the colour 
tensor Cj is equivalent to taking the tensor product of this colour tensor with 03^. These 
^This is true for the t-channel basis set we are using for this process. 
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identities are straightforward to prove; for example, if we label the two exchanged 
gluons as c and d, the virtual quark as i and the virtual gluon as j, we get the 
following group algebra for C3 x C3, 



(2.34) 



where we can see the meaning of the x notation. By using equation (2.29) and 
equations (A.10-A.16) we can prove equation (2.33). We can now use these identities 
and deduce that 





C3, 


M^B) _ 
^2 — 


C3 X C2 


M^B) _ 
^3 — 


C3 X C3 


Hence 




cm ^ ^ 


n(AB) _ 



1 

2' 



'c 



( 









1 



4 

We -4 

4JVc 



Turning now to C^^^\ we see that 



(Al) 



1 \ 

2 

4 



(2.35) 
(2.36) 

(2.37) 



(2.38) 



(2.39) 



and that 



^2 



— dB2ct'^jitiAiij: 



dB2c fji 2 
1 



2N, 



2N, 



dB2c ^lA) 



■C2, 



(2.40) 



where we have expanded the pair of SU(3) matrices t^^ tfj and used the fact that = 0. 
Finally, working in a similar way. 



2K 



C3- 



(2.41) 
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Note that these last three results are not the same for C^^^) . For this case the group 
algebra gives (where we note that = Ca), 



p{B2) 



— t'^Adijd'i'fBic'i'fj2C — ~2^'^'' 

— fijd fsic fj2C — -;rCs- 



(2.42) 
(2.43) 

(2.44) 



For the first of these results we have used the fact that fBcdf2cd — NcSb2, for the 
second of these results we have used the fact that djdi fiBc fc2j — —dB2d and for the 
third of these results we have used the fact that fjdi fiSc fc2j — —^fdb2- Therefore we 
find that 



C 



{B2) 



( 



\ 






Ca 












\ 





2Nc ^ 
-2^ / 







(2.45) 



f y 



(2.46) 



After similar, albeit longer, calculations we find 



.(^2) 



^2 



1A2) 



^C^i2c fsjd tu tiA - 
— fj2cfBjdtiitiA 



- 4 

1 



(2.47) 
(2.48) 

(2.49) 



and hence the result 



,(A2) 



v 






1 



4 



4 

4 / 



(2.50) 



4iVc 



and so we find the colour mixing matrix for the process qg — > qg, in basis 2.26, is 

( CpC^^^^ + CaC^^^^ -\{oc + i) ^ 

X -f(« + 7) ' (2-51) 



C 



-(a + 7) 



4Nc 



(a + 7) 



X 
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where 

X = ^(«-7)-^C(^^^ + %^"^^. (2.52) 

In this equation, as in the previous section, (^'^^^ denotes the dynamical piece of the 
relevant diagram. We can proceed in a similar way and calculate the mixing matrices 
for any partonic process that may interest us. The bases and the mixing matrices for 
all the processes used in this thesis can be found in appendices D and F and we will 
use these results when we study the resummation of rapidity gap events in chapters 4 
and 6. 



2.4 Regge theory and the pomeron 

Regge theory [10] provides a framework to describe scattering amplitudes in the Regge 
limit s ^ —t. Whilst in this thesis we are primarily concerned with the perturbative 
description of rapidity gap processes, we will use some of the ideas of Regge theory in 
chapter 3 to examine difTractivc processes and so in this section we will outline some 
of the basic concepts. Further details can be found in [11-13]. 

The whole idea of Regge theory, starting from the analytic properties of the S- 
matrix, is to extract the high energy behaviour of scattering amplitudes in a model 
independent way. The fundamental result is that the dominant contribution to the 
high energy scattering amplitude has a general form, which has the interpretation of 
the exchange of objects known as reggeons in the t-channel. It can be shown that if 
the high energy cross section in the Regge region can be described only in terms of 
reggeon exchanges, then the total cross section will behave hke 

atot oc (2.53) 

where a{0) is known as the intercept of the Regge trajectory (the Regge trajectory is 
a straight line in the spin- mass squared plane of the exchanged mesons). The fits of 
Chew and Prautschi [14] to meson data showed that a{0) — 0.55, and this observation 
of a{0) < 1 was found to be true for other exchange particles. Hence the total cross 
section should decrease with s, and vanish asymptotically. However when the total 



Chapter 2. QCD at the frontier 



34 



cross section was measured experimentally, it was found that it slowly increased at 
high s. Therefore reggeon exchange is not the whole story, and a new trajectory is 
needed to describe the data, with an intercept greater than unity, a{0) > 1. This new 
trajectory is known as the pomeron trajectory, or the pomeron. 

Prom fits to total pp and pp data, Donnachie and Landshoff [15] extracted the 
pomeron intercept and found 

a(0) ~ 1.08. (2.54) 

This is normally known as the "soft" pomeron intercept. We will use the ideas of 
Regge theory in chapter 3, when we study diffractive processes at the Tevatron and 
we will describe the diffractive interaction of two hadrons by the exchange of both 
pomerons and reggeons. 

2.5 Event generators 

Monte Carlo event generators have many uses in particle physics, for example estimat- 
ing the backgrounds to measured processes and calculating production rates of exotic 
particles. In this thesis we will use the Monte Carlo event generator HERWIG [16, 17] 
on two occasions as a calculation tool and in this section we will give an overview of 
the theoretical methods used. Our discussion in this section on Monte Carlo event 
generators will hence refer to the operation of HERWIG. 

The idea of a Monte Carlo event generator is to give a simulation of a particle 
physics event, starting from the initial interaction of the colliding particles and cul- 
minating in the angular distribution and energy of (colourless) final state particles. 
This is still a level removed from what is seen in particle detectors, and the inclusion 
of detector effects is possible, but it is largely sufficient for the research work in the 
following chapters. To simulate the event, HERWIG starts with a hard subprocess 
and evolves the resulting off-shell partons into colourless final state particles; to do 
this HERWIG is separated into a number of distinct phases. 

For the process of two hadrons interacting to produce a number of hadronic jets, 
HERWIG starts by calculating the hard subprocess, producing two final state partons 
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using the QCD factorisation theorems with a leading order calculation and the parton 
distribution functions of the hadrons. At this point wc have selected a QCD 2 — > 2 
subprocess, but any perturbatively calculable process of the standard model or any of 
its extensions is permissible. The result of the calculation is two initial state and two 
final state partons, with their associated 4-momenta, and the total cross section for 
the process, which is calculated by evaluating the integral, 



where pt is the transverse momentum of the outgoing partons. The cut pt,min is 
necessary to prevent the integrand approaching divergences in the matrix element and 
the Landau pole in the coupling. 

Once the momenta of the final state partons produced by the hard event have 
been generated, the initial and final state partons are evolved through perturbative 
branching, and all emitted partons are themselves evolved, until all partons reach an 
infrared hadronisation scale which characterises the onset of non-perturbative physics. 

Final state parton shower 

The off-shell partons produced by the hard subprocess subsequently evolve with the 
emission of QCD radiation. The exact calculation of the matrix element for processes 
with many final state particles is not possible, and we must use a branching algorithm 
which correctly takes into account the regions of phase space in which QCD radiation 
is enhanced. These regions are associated with kinematical configurations in which the 
relevant matrix elements are enhanced, and correspond to emission of a soft gluon, or 
when a gluon or a quark splits into two coUinear partons. The aim of a parton shower 
algorithm is to identify and sum up the leading behaviour in these regions. 

The parton showering is controlled by a set of Sudakov form factors, which en- 
code the probability that a parton with a virtual mass scale t will evolve without 
resolvable branching to the lower virtual mass scale ti. If branching does occur, split- 
ting functions are used to compute the momentum fraction and virtual mass scale 
of the products. The Sudakov form factors take into account coUinear enhancements 




(2.55) 



TTl 
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and also soft enhancements using the coherent branching formalism and the parton 
shower is terminated when all partons reach the cut-off scale Iq. The virtual mass 
scale variable, t, which is evolved to to^ is a combination of parton energy and branch- 
ing opening angle; this choice of the evolution variable is the correct one to include 
both soft and coUinear enhancements. There is an excellent explanation of the parton 
shower formalism, and the application to Monte Carlo event generators, in [1,18]. 

Initial state parton shower 

The final state parton shower is called a forward evolution algorithm, as in every 
step the partons move to a lower virtual mass scale. For the initial state parton 
shower, it is more convenient to start with the (most negative virtual) partons which 
particpate in the hard event and evolve backward to the lower virtual mass scales 
of the partons in the interacting hadrons. In this backward evolution formalism, the 
parton distribution functions are used as part of the input, to "guide" the evolution 
to the correct distribution of partons. 

Hadronisation 

Using the parton shower formalism, all partons are evolved until they reach the cut- 
off scale to- It is at this point that we enter the long- wavelength, non-perturbative 
regime where the partons group themselves into the colourless hadrons we observe; 
this process is included in HERWIG using a phenomenological model. 

A property of the parton showering is that the flow of momentum and quantum 
numbers at the hadron level tends to follow the flow at the parton level. This hypothe- 
sis is known as local parton-hadron duality (LPHD), and underlies the cluster hadroni- 
sationm.ode\ [19,20] used by HERWIG. In this model, clusters of colour singlet partons 
form after the parton shower, inheriting the momentum structure of their constituent 
partons. The process occurs by the splitting of gluons into quark/antiquark pairs, 
and then neighbouring quarks and antiquarks form colour singlet objects. Unstable 
hadrons are then allowed to decay according to experimentally determined branching 



Chapter 2. QCD at the frontier 



37 



ratios and this process, which tends to disfavour heavy mesons and baryons, provides 
a good model of observed final states in e"'"e~ collisions. 

HERWIG also includes a model of the underlying soft event, in which the hadron 
remnants undergo secondary soft interactions. Once the final state has been deter- 
mined, the experimental cuts for the analysis of interest can be applied and observables 
calculated. It is also possible, as we briefly mentioned earlier, to add a further phase 
of detector simulation, where the experimental signature of the produced final state 
is simulated. Such a process allows for detailed comparisons to be made with experi- 
mental observation. 

We will use HERWIG in chapter 3, where we are especially interested in the 
effects of parton showering and hadronisation and again in chapter 6, where we will 
be interested in the computation of the hard event. 



The aim of this thesis is to examine the physics of rapidity gap processes, and make 
predictions to compare to recent experimental data; in this section we will give a brief 
overview of rapidity and rapidity gaps in particle collisions. The scattering of two 
hadrons provides two beams of incoming partons, with a spectrum of longitudinal 
momenta described by hadronic parton densities. The parton centre-of-mass frame is 
boosted with respect to the incoming hadron frame, and it is natural to describe the 
final state in terms of variables which have simple transformation properties under 
longitudinal boosts. To do this we describe the final state in terms of transverse 
momentum pti rapidity y and azimuthal angle 0, and a general 4- vector is written 



2.6 Rapidity gaps 



p 



{E,Pl,P2,P3), 



{rrit cosh y, pt sin 0, pi cos 0, rrit sinh y) , 



(2.56) 
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Figure 2.2: The experimental signature of a "gap-jet-gap" event. The dark circles 
denote two jets, and the gap regions are separated from the central dijet system by 
the dashed lines. 



where we define the transverse mass for a particle of mass m hy rrit = \/j^ + rn?. 
Rapidity is defined as 

y = ^log 

= - log 

2 ^ 

and is a measure of polar angle and speed. 

Pseudorapidity is often more useful in practice, which is defined 

?7 = -logtan . (2.59) 

Rapidity and pseudorapidity coincide in the massless limit, but pseudorapidity is more 
convenient in that the polar angle 6 can be measured directly from the detector and 
we require no knowledge of the particle mass. Rapidities (and approximately pseudo- 
rapidities) are additive under z-axis Lorentz boosts, and (pseudo)rapidity differences, 

Ar] = \r]i - r]2\, (2.60) 

are boost invariant. 

A rapidity gap event is defined as an event producing jets, with the region in 
pseudorapidity between any two jets. At], being devoid of particle activity. This 
interjet region, known as a gap, is often symmetrical in azimuthal angle 0. The precise 



E+P3 

1 + V COS 6 
1 — V cos 6 



(2.57) 
(2.58) 
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7] 

Figure 2.3: The experimental signature of a "jet-gap-jet" event. The dark circles 
denote two jets, and the gap region is bounded by the dashed lines. 

definition of the gap depends on the experimental geometry, and is experimentally 
defined by a cut on the total transverse energy flowing in the gap region. In this thesis 
we will consider two kinds of rapidity gap processes: 

• "gap-jet-gap" processes in diffractive hadron collisions. In this class of process, 
the final state consists of a central (in rapidity) dijet system, separated from the 
intact colliding hadrons by a rapidity gap on each side. Hence the "gap-jet-gap" 
experimental signature, which is illustrated in figure 2.2. We shall study these 
processes in chapter 3, where the rapidity gap is produced by the exchange of a 
pomeron between the hadron and the central system. 

• "jet-gap-jet" processes in photoproduction. In this class of process we see a 
rapidity gap between the two dijets. We shall study these gaps-between-jets 
processes in chapters 4, 5 and 6, where we will make perturbative calculations of 
the rate of gap-event production and make comparisons to experimental data. 
The experimental signature is illustrated in figure 2.3. 

For both classes of rapidity gap process, we will find that the predictions of the methods 
that we use will give a good description of the data. 
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2.7 Summary 

In this chapter we have outhned some of the tools we shall be using in this thesis. The 
features of asymptotic freedom and factorisation give QCD immense predictive power, 
and one of the central aims of our work is to examine the consequences of factorisation 
for the perturbative calculation of rapidity gap processes at HERA. An important 
element of our calculations will be colour mixing, and in section 2.3 we derived a set of 
process-dependent QCD mixing matrices for a fixed basis. We will use these matrices 
in chapter 6. We then described the basic ideas of Regge theory, event generators 
and rapidity gaps, in preparation for the use of these ideas in our calculations. So in 
conclusion, the research work in this thesis deals with the calculation of rapidity gap 
processes in QCD and builds on the fundamental tools presented in this chapter. 
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Chapter 3 

DifFractive dijet production 



3.1 Introduction 

In this chapter we will study hard difTractive processes in hadron-hadron collisions. To 
be specific we will make detailed calculations for the process pp — > p + JJX + p, where 
JJX denotes a centrally produced cluster of hadrons, containing at least two jets. The 
process is theoretically diffractive in the sense that both the initiating hadrons remain 
intact in the collision and experimentally diffractive in the sense that the initiating 
proton remains intact, whilst the antiproton either remains intact or dissociates to 
a low-mass system. In both cases the initiating hadrons sufi^er only a small loss of 
longitudinal momentum, and the process is hard in the sense that the central subpro- 
cess takes place at high momentum transfer. The central (in rapidity) jet-producing 
system is separated by a rapidity gap from each of the interacting hadrons, giving the 
experimental signature of "gap-jet-gap" events and the term hard double diffraction. 

Space-time arguments [21] suggest that hard events are well localised in space 
and time, and therefore the effect of the incoming particles is to act independently of 
the hard event. This has led, in analogy to conventional QCD factorisation theorems, 
to the concept of a diffractive parton density and the hard diffractive factorisation 
theorems. Diffraction factorisation has been proven for diffractive DIS, but for hard 
diffraction in pure hadron collisions counterarguments exist which predict that the 
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factorisation will be violated at the Tevatron [21-24]. The manifestation of this viola- 
tion is the invalidity of using diffractive parton densities, obtained in DIS experiments 
at HERA, in hadron-hadron collisions at the Tevatron and we should expect to see 
an overestimation of predicted cross sections. The extent of the breakdown of fac- 
torisation is still under investigation and has been embodied in so-called gap survival 
factors, which will have direct relevance in this work. 

The standard factorised model of double diffraction draws its inspiration from 
Regge theory. In the Ingelman-Schlein (IS) model [25], diffractive scattering is at- 
tributed to the exchange of a pomeron, which is defined as a colourless object with 
vacuum quantum numbers. The model assumes that each of the diffracting hadrons 
"emits" a pomeron, which then collide to produce the central dijet system. This is 
illustrated in figure 3.1. Implicit in the model is Regge factorisation, which writes a 
diffractive parton density as the product of the pomeron parton density and a pomeron 
flux factor. The double diffractive process then proceeds by the QCD interactions of 
the emitted pomerons. This kind of event mechanism is known as double pomeron 
exchange. 

In this chapter we will use the Ingelman-Schlein model of double diffraction to 
make predictions for gaps-between-jets processes at the Tevatron. The CDF collabo- 
ration has recently presented analyses of such events [26] . We will include the effects 
of fragmentation using the diffractive event generator POMWIG [27], and use the 
pomeron parton densities from HERA DIS experiments [28]. Previous calculations of 
this type of event [22,29] have led to an overestimation of the dijet production cross 
section by several orders of magnitude. In this chapter we will show that a combi- 
nation of the HERA parton densities and a gap survival factor [30] consistent with 
theoretical estimates [31,32] is sufficient to describe the recent CDF data. Central 
dijet production within a factorised model has been studied before [22,29,33,34] and 
the research work in this chapter is published in [35]. 

This chapter is organised as follows. We start by introducing the experimental 
analyses in section 3.2, and section 3.3 is a detailed discussion of the IS model and the 
concept of diffractive parton densities. Section 3.4 then discusses the diffractive event 
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Figure 3.1: Dijet production in the factorised model. The zig-zag hues denote the 
exchanged pomerons, and the double-hned arrows represent the pomeron remnants. 

generator POMWIG, and how it is constructed from HERWIG, and section 3.5 gives 
our results. Finally, in section 3.6, we draw our conclusions. 

3.2 Diffractive dijet production at the Tevatron 

The first observation at the Tevatron of dijet events with a double pomeron exchange 
topology was published in 2000 by the CDF collaboration [26] . In this section we will 
describe the experimental signature of these events, the observables and cuts used, 
and finish with a description of the out-of-cone corrections implemented by CDF to 
"undo" the jet broadening effects of fragmentation. 

3.2.1 Diffractive events and the CDF analysis 

The diffractive data used in this work was collected by the CDF collaboration when 
the Tevatron was colliding protons and antiprotons at a centre-of-mass energy y/s = 
1800 GeV. The events were obtained by making a further analysis of the single 
diffractive (SD) data set; SD meaning that only one interacting hadron undergoes 
diffraction. In this analysis, events with a rapidity gap on the antiproton side were 
selected by requiring a fractional momentum loss of the antiproton, ^p, to satisfy 
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Figure 3.2: The experimental signature of a "gap-jet-gap" event. The dark circles 
denote two jets, and the gap regions are separated from the central dijet system by 
the dashed lines. 

0.035 < < 0.095, by detecting the antiproton in a forward Roman Pot Spectrom- 
eter (RPS). The requirement was then made of at least two centrally produced jets 
of transverse energy i?™™ > 7 GeV and a 4-momentum transfer from the antiproton 
squared of \t\ < 1 GeV^. This SD sample (of 30,439 events) was then used to search 
for double diffractive events by demanding a rapidity gap on the proton side. There 
is no RPS on this side of the CDF detector and the events were selected by directly 
demanding a region devoid of particle activity in the pseudorapidity region between 
the central dijet system and the outgoing proton. This method is not a very accurate 
way to determine whether a given event is indeed a genuine double diffractive event, 
and consequently gives a large uncertainty on the the fractional momentum loss, C,p, 
of the proton. As a result, CDF found that the C,p distribution was concentrated in 
the region 0.01 < < 0.03 and used this region for the cut. The subset of the SD 
sample satisfying this criteria is the double diffractive sample (denoted DD, in which 
two interacting hadrons undergo diffraction). The double diffractive experimental sig- 
nature is illustrated in figure 3.2. In summary, the CDF cuts used to select the double 
pomeron exchange events were 

• The antiproton fractional energy loss C,p satisfies 0.035 < < 0.095. Such a 
condition corresponds to a rapidity gap on the antiproton side. This cut is made 
by tagging the antiproton. 
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• The proton fractional energy loss satisfies 0.01 < < 0.03. Such a condition 

corresponds to a rapidity gap on the proton side. Note that the proton is not 
tagged at CDF and the cut is made by looking for a rapidity gap on the proton 
side. 

• \tp\ < 1 GeV^. \tp\ is the four-momentum transfer from the antiproton squared. 

• The event must have two or more jets in the pseudorapidity region 



and at least two jets must have a minimum transverse energy of = 7 GeV 
or 10 GeV. The CDF collaboration define the outgoing proton direction as the 
positive rj direction, also known as "East" . 

The CDF analyses presented the total dijet cross section at > 7 GeV and 

^mm y GeV. They also presented event distributions in ^p, ^p, the mean jet 
transverse energy E^, the mean jet rapidity r)*, the azimuthal angle separation of the 
dijets and the dijet mass fraction Rjj. The last quantity is defined as the mass of 
the dijet system evaluated using only the energy within the cones of the two leading 
jets, Mj°^^, divided by the mass of the central system, — \J s^pCp- 

3.2.2 Jet finding and the cone algorithm 

The classification of a hadronic final state into one or more "jets" of hadronic particles 
is a task performed by a jet algorithm. Modern algorithms have their roots in the 
Sterman- Weinberg jet definition [36], and today several difi^erent algorithms exist. 
The method used by the CDF collaboration in their double diffractive dijet production 
analyses was the cone definition. In this algorithm, a jet is defined as the concentration 
of all particles inside a cone of (77, 0) space, of radius R in the (77, 0) plane. 



The jets are reconstructed to place as much of the transverse energy as possible inside a 
cone, with a minimum energy of e'^"* GeV. Overlapping cones are handled by excluding 



-4.2 < 77 < 2.4, 



(3.1) 



R = V(A?7)2 + (A0)2. 



(3.2) 
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Figure 3.3: A 3 jet event, defined using the cone algorithm with R = 0.7, from the 
CDF experiment. The clusters of particles defined as a jet are ringed by the cone 
radius. This figure is taken from [38] with the authors' permission. 



jets with more than a certain fraction, normally known as the overlap parameter, of 
their energy coming from other jets. By defining R in terms of Arj (other choices 
could be in terms of A^, for example) we obtain a jet measurement invariant under 
longitudinal boosts. Note that any cross section is dependent on the choice of the 
cone radius R. In this work we will follow CDF and set R = 0.7. Figure 3.3 shows 
the application of the cone algorithm to a 3 jet event at the Tevatron. The clusters of 
particles defined as a jet are ringed by the cone radius in this figure. 

Once a concentration of particles has been identified as a jet, the transverse energy 
of that jet is given by the sum of all the transverse energies of the constituent particles, 

4^* = Et- (3.3) 
Other kinematic variables have similar definitions, of which 



jet _ E^gjct-^T^ 
, jet _ ^igjet Et 



1 J,* 



is one possible form. 
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A [GeV] 


B 


C [GeV-'] 


0.4 


22.999 


0.915 


0.00740 


0.7 


8.382 


0.846 


0.00728 


1.0 


3.227 


0.832 


0.00817 



Table 3.1: The out-of-cone correction parameters, as derived from Run I Tevatron 
data [39]. We have shown the out-of-cone correction parameters for three different 
values of the radius, although only R — 0.7 is relevant to the current work. 

3.2.3 Out-of-cone corrections 

The CDF collaboration have applied so-called out-of-cone (oc) corrections to the 
diffractive dijet data sample. These corrections are applied as the jet clustering algo- 
rithm may not include all of the energy from the initiating partons, and some of the 
partons generated during fragmentation may fall outside of the jet algorithm cones. 
Therefore the out-of-cone corrections add transverse energy to the jets to account for 
this missing pt and "correct" the particle-level jet energies back to the parton level. 
This procedure (which is difficult to justify) means that CDF no longer presents a true 
observable. These corrections depend on the parton fragmentation functions, and are 
totally independent of the CDF detector. The out-of-cone corrections take the form 
of a simple addition to the jet pt, and were derived from Run I data. The additional 
energy is given by 

Pt = A{R) (1 - S(i?)e-^(«)^*) , (3.5) 

where pt is the particle-level jet energy and the functions A{K),B{R) and C{R) depend 
on the cone radius (R) and have the values given in table 3.1. The impact of the out- 
of-cone corrections to the work in this chapter is that it is necessary to implement the 
out-of-cone corrections in the Monte Carlo simulation and compare the hadron level 
Monte Carlo cross sections with the corrections to the CDF data, and not (as it may 
first appear) compare the data to the hadron level cross sections. Nevertheless, the 
use of Monte Carlo event generators allows the impact of the fragmentation process 
to be assessed. 
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3.3 Diffractive factorisation and the Ingelman-Schlein 
model 

We will now examine the mechanisms for hard double pomeron exchange dijet pro- 
duction. In this work we will use the factorised Ingelman-Schlein model, which takes 
its inspiration from Regge theory. We will also describe the non-factorised model of 
lossless jet production [23,40-42], but we will not make any calculations using this 
model in this thesis. 

3.3.1 Hard diffraction 

The pomeron is the object thought to be responsible for rapidity gap processes. De- 
fined as an exchanged object carrying the quantum numbers of the vacuum, it is 
postulated to exist in many forms - from the soft pomeron of Regge theory to the 
perturbative QCD pomeron embodied by the famous Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation. Close to 20 years ago Ingelman and Schlein proposed that it was 
possible to probe the content of the pomeron by looking at the diffractive production 
of high pt dijets. In so-called hard diffraction the momentum transfer across the rapid- 
ity gap is small, with the gap corresponding to the exchange of a colour singlet object, 
but a high process occurs between the exchanged object and the other hadron. 
This is be contrasted with diffractive hard scattering, in which the colour singlet ob- 
ject is exchanged with a high momentum transfer across the gap. Therefore the basic 
idea of Ingelman and Schlein [25] is that the exchange mechanism for the pomeron is 
the same for hard diffraction and soft diffraction, and the pomeron interacts through 
its partonic constituents. These ideas took shape when looking at single diffraction, 
p + p ^ p + X , and led to the following picture: 

• The diffracting antiproton will emit a pomeron with a small momentum trans- 
fer t. The pomeron is sometimes said to "float" out of the "parent" antiproton. 

• The pomeron, carrying a partonic content, then interacts with the partons in 
the proton at high Q^. 
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Figure 3.4: An illustration of Regge factorisation. Figure from [43]. 



The standard Regge single diffractive cross section (single diffraction with a diffracted 
proton) is 

^ = /p/p(e,t)afV,t), (3.6) 

where s' = is the centre-of-mass energy of the pomeron-antiproton system, o"f^(s', t) 
is the total pomeron-proton cross section {a{pF X)) and /p/p(^,t) is the pomeron 
flux factor. We have illustrated this Regge factorisation of the cross section in fig- 
ure 3.4. This equation can be treated as a definition of the pomeron fiux, and allows 
us to write^ the single hard diffractive cross section for the production of dijets as, 



SD 



dxpdtdxidx2 



d'^asD 



d^(T{pF jj + X) 
dxidx2 



(3.7) 



_a{p¥ X) dx^dt 

where Xp is the momentum loss fraction of the antiproton (i.e. momentum fraction of 
the pomeron to the antiproton), t is the 4-momentum transfer-squared and cr{pF 
jj + X) is the proton-pomeron hard-scattering cross section. The pomeron flux factor 
(in square brackets) depends on xp and t. The pomeron-proton hard scattering cross 
section depends on the momentum fractions, Xi and X2, of partons in the parton and 
the pomeron respectively. Therefore this differential cross section is given by 

d^a{pF^jj+X) ^ n'2\f ( ri2N^^(u ^ dijets) 

dx,dx2dQ^ = ^ ^ ^ ^ dQ^ ' ^'-'^ 

where fi/p{xi,Q^) and fj/p{x2,Q'^) are the parton densities of the proton and the 
pomeron respectively. In the original IS model it was assumed that the pomeron was 
a purely gluonic object with a parton density that was and t independent, and two 



-"^In this section we use ^ and xp interchangably. 
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simple forms where proposed, 

xfg/w{x) = &x{l-x), (3.9) 
xfg/r{x) = 6(l-x)^ (3.10) 

the so-called hard and soft gluon densities respectively. However, these parton den- 
sities (combined with the (soon to be described) DL flux factor) failed to describe 
Tevatron data. In this work we will use modern pomeron parton densities obtained 
from HI fits to single-diffractivc structure function data [28]. 

The most commonly used pomeron flux factor was suggested by A. Donnachie 
and P.V. Landshoff [15], now known as the DL pomeron flux, 

^M^p) = 4^[^iWf (-) ' (3-11) 

where Fi{t) is the proton form factor, given by 

, ^ Ami- 2.8t f 1 \ , , 

where rrip = 0.938 GeV is the proton mass and /?o = 1-8 GeV~^ is the pomeron-quark 
coupling. The pomeron trajectory ap{t) has the form Q;p(i) = af{0) + a'^t, where ap{0) 
is known as the pomeron intercept. The DL model was fltted to diffractive data, which 
gave 

a'-p = 0.25 GeV-^ 
q;p(0) = 1.08. (3.13) 

This form of the flux was found to poorly describe Tevatron data and in this work 
we use the POMWIG [28] parameterisation of the flux, with a harder pomeron inter- 
cept of Q;p(0) — 1.20, coming from more recent HERA flts. Apart from this harder 

intercept (which is important when we come to compare to the Tevatron data), this 
parameterised form of the flux is approximately the same as the DL flux. 



3.3.2 Diffractive structure functions 

The conventional way to write the single diffractive dijet production cross section is 
in terms of diffractive structure functions. As we have discussed in the last section. 
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this SD cross section can be written as 



^'^"^ = /p/p-(e, t) Y: Ui^, Q')fj/M g^) ^^fa -;fj^^^) , (3.14) 



where we have relabelled the momentum fraction of a parton in the proton as x, and 
the momentum fraction of a parton in the pomeron as /3 = Xp/xp. We have also 
denoted the momentum fraction of the pomeron with respect to the antiproton as 
The diffractive structure function F^- of the antiproton is 

K%{^,t,(3,Q') = fF/p{U)fi/AP,Q'), (3-15) 

which is the product of the pomeron flux and the pomeron parton density (this level of 
factorisation, inspired from Regge theory, is unproven in hard diffraction and is known 
as Regge factorisation [25]. It is illustrated in figure 3.4). This structure function is 
sometimes known as F^p'*^^ ■ The SD dijet cross section is then written 

This expression is known as the (hard) factorisation in the diffractive structure func- 
tion. The t-integrated diffractive structure function is known as F^^^^ with definition 

= /p/p-(0/.7p(/3,Q')- (3.17) 

where fw/piO is the t- integrated pomeron flux. Performing the ^ integration gives us 

the definition of F.^P\ 
t/p ' 

<f (^'Q') = d^MOfM^,Q'). (3.18) 

This final quantity, which now only only depends on the momentum fraction of the 
parton relative to the diffracting hadron (x) and Q^, is often called the diffractive 
parton density in analogy to the non-diffractive parton densities. 
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3.3.3 Double pomeron exchange 

The double diffractive process in hadron-hadron coUisions is formulated as doubly 
occurring single diffraction. Therefore the process can be modelled using (a copy of) 
the IS model for both the proton and the antiproton. The diffractive event therefore 
proceeds by both the proton and the antiproton emitting a pomeron with a small 
momentum transfer, and these two pomerons interacting with each other through their 
partonic content at high . This QCD process produces two centrally produced jets, 
with a gap on each side of rapidity ( "bridged" by the pomerons) to the intact parent 
hadrons. This type of event topology is known as double pomeron exchange or DPE 
(in such terminology, single diffractive events proceed by single pomeron exchange 
and only a single rapidity gap is seen in the final state). The DPE event topology is 
depicted in figure 3.1. 

The DPE cross section can be written 



J dijet 
DPE 

drj^drj^dp 



y{pp^ p+ JJX +p) ^ J d^p J d^pFp/p{QFp/p{Q 



E/5pV(/^p)/5p-/./p(/3p-)^(^J - kl), (3.19) 

where Fp/p{^) = Fp/p{^) is the pomeron flux factor, /5 is the fraction of the pomeron 
momentum carried by the parton entering the hard scattering and fi/p{f3) is the 
pomeron parton density function for partons of type i. The rapidity of the outgoing 
partons are denoted 773 and 774, their transverse momentum is p± and ^^^{ij — ^0 
denotes the QCD 2-to-2 scattering amplitudes. The parton transverse momentum, 
is equal to the jet transverse energy Et at the parton level. In this model, the par- 
tonic content of the two pomerons which does not participate in the hard event forms 
pomeron remnants in the final state; hence DPE event topologies are inclusive. In 
this work we will use the ^-integrated flux parametrisation of POMWIG (discussed in 
section 3.4) with a pomeron intercept of q;p(0) = 1.20 as found by HI [28] in their fits 
to ^2^^^^ data, instead of the softer ap(0) = 1.08. 

We should also briefly mention the non-factorising models [23,40-42] for diffrac- 
tive central dijet production. In these models, all the momentum lost by the diffracting 
hadrons goes into the hard event, giving an exclusive event topology and no pomeron 
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Figure 3.5: Dijet production in the non-factorised modeL This model is characterised 
by an exclusive event topology. 

remnants. This topology is depicted in figure 3.5^. A significant non-factoring con- 
tribution would manifest itself as a peak (at unity) in the distribution showing the 
available centre-of-mass energy which goes into the jets, known as the dijet mass 
or Rjj] a feature which is absent in the CDF data we consider in this thesis. However, 
this need not be the case for the Tevatron Run II due to higher luminosity, and this 
ought to be a good place to look for such evidence of a non-factorised contribution. 
Consequently this thesis does not consider such models. 

3.3.4 Breakdown of factorisation and gap survival 

The notion of a parton density in QCD is only useful if it has the property of univer- 
sality. This means that once a parton density has been extracted from one process, it 
should be able to predict others processes with great accuracy. The current pomeron 
parton densities have been extracted from HERA DIS experiments, where the struc- 
ture function Fg^*-'^^ was measured. However when it came to testing HERA parton 
densities at the Tevatron, for hadron-hadron collisions, it was found that they overes- 
timated the data [22,29]. This violation of factorisation has been understood in terms 
of simple models, which indicate that the rapidity gaps will be filled with secondary 

^In [42] the possibility that there is additional radiation into the final state from the 
gluons in figure 3.5 is considered. 
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interactions, from spectator partons in the beam. These models also indicated that 
the process of spoiling the gap can be approximated by a simple overall multiplicative 
factor [30-32,35,41,42,44], with a weak process and a/s dependence. For some 
processes at HERA, this factor has been estimated to be ~ 0.6 [45], and some Teva- 
tron processes have been estimated to have a 5"^ ~ 0.1 [31,34]. In recent times various 
models have been postulated for the calculation of the gap survival factor [30], and 
seem to support this latter figure when applied to Tevatron observables. Therefore, 
as this work involves using HERA parton densities at the Tevatron, we expect that 
factorisation will be violated and we will include in our predictions an overall gap 
survival factor, S'^. We will extract this factor by fitting our overall cross sections to 
data. 

3.4 POMWIG - a Monte Carlo for Diffractive pro- 
cesses 

The event generator POMWIG [27] is used in this work to study the hadronic diffrac- 
tive process. The modifications to the HERWIG Monte Carlo event generator [16, 17] 
to study diffractive events are very simple once it is noticed that pomeron exchange 
events in hadron-hadron collisions look very much like the resolved part of photopro- 
duction in electron/proton collisions. In the latter type of event, the electron radiates 
a quasi-real photon according to a flux formula. This photon is then considered to 
have a partonic content which interacts, through QCD processes, with the partonic 
content of the proton. Therefore we can study pomeron exchange events in POMWIG 
by replacing the photon flux with a pomeron flux and by replacing the parton density 
of the photon with that of the pomeron. This philosophy is illustrated in figure 3.6. 

POMWIG is available as a modification package to HERWIG, and adopts the 
philosophy of minimal changes to the HERWIG structure. The generator uses the 
factorised IS model at the parton level, with pomeron fluxes and parton densities 
we will discuss later in this section, and allows the final state partons to evolve into 
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Figure 3.6: The similarity between the resolved part of photoproduction and pomeron 
exchange events. Figure from [27]. 

hadrons using a parton shower and hadronisation model, as discussed in the introduc- 
tion to this thesis. The inclusion of the fragmentation physics through a Monte Carlo 
simulation allows the phenomenological impact to be assessed. 

In POMWIG, the pomeron and reggeon fluxes are parameterised as 



where ap(t) = ap(0) + a'^t and a;M(t) = ^^(O) + a'^t. The constant is chosen 
so that ^2^^^"* matches the HI data at xp = 0.003 and similarly the constant is 
found from data. The flux parameters used are those found by the HI collaboration, 
assuming there is no pomeron/reggeon interference contribution to ^2^*^'^^ The default 
parameters are given in table 3.2. 

In this work we use the HI leading order (LO) pomeron fits to ^2^*'^'' measure- 
ments [28]. The measurement of Fg^*'^^ is quark dominated and gluon sensitivity enters 
only through scaling violations. Hence the gluon density has a large uncertainty of 
around 30% in the relevant region. This is important for the gluon dominated DPE 
process. The gluon densities are illustrated in figure 3.7. The mean value of (3 relevant 
at the Tevatron is in the region of 0.3 to 0.4. Fit 4'^ has a gluon content that is heavily 
■^For details of the fits, see [28] . 




(3.20) 
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Quantity 


Value 




1.20 




0.57 


Op 


0.26 


"r 


0.9 


Bp 


4.6 


Br 


2.0 


Cr 


48 



Table 3.2: The default POMWIG parameters. 
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FigTirc 3.7: The gluon densities in the pomeron as fitted by HI. The fits arc evaluated 
at = 50 GcV^, the typical scale of the jet transverse energy at CDF, and are 
multiplied by the pomeron flux factor at = 0.1. 
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suppressed relative to fits 5 and 6 (hence fit 4 is quark-dominated) , and fit 6 is peaked 
at high (3. Fits 5 and 6 arc now the favoured fits to describe HI data [28]. 

POMWIG also allows us to include the effect of non- diffractive contributions 
through an additional Regge exchange, which we refer to as the reggeon contribution. 
This is expected to be important in the region of explored at the Tevatron. Following 
HI, we estimate reggeon exchange by assuming that the reggeon can be described by 
the pion parton densities of Owens [46]. This contribution is added incoherently to the 
pomeron contribution. For further details of the implementation we refer the reader 
to [27]. 

3.5 Results 

In this section we compare our results using the IS model and POMWIG with the 
Tevatron DPE data. We start by testing POMWIG against an independent IS cal- 
culation, which we performed by the direct integration of equation (3.19) using the 
HI pomeron fiux and parton density. All results using POMWIG are then shown at 
the parton level, at the hadron level and at the hadron level with out-of-cone correc- 
tions implemented. To assist in the interpretation of the results, we then comment on 
results at the parton shower level. 

3.5.1 Hadronisation and parton shower effects 

Figure 3.8 demonstrates the agreement between POMWIG and an independent cal- 
culation for the total DPE cross-section as a function of the minimum jet transverse 
energy i?™". These curves were produced using HI fit 5 and contain only the pomeron 
exchange contribution. More interestingly, we see that the curve which describes the 
total cross section after the effects of fragmentation (hadronisation and parton shower- 
ing) have been included shows a significant reduction relative to the naive parton level 
calculation. This suppression effect may be understood by observing that the effect 
is a shift in the cross-section by /\Et = 2 GeV. This is a direct consequence of the 
broadening of the jet profile by the parton shower and hadronisation. The reduction 
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Figure 3.8: The total cross-section at both parton and hadron level, plotted as a 
function of the minimum transverse energy of the dijets. 



is lower for the quark dominated fit (HI fit 4) since quark jets tend to have a narrower 
profile. In practice, this hadronisation suppression will reduce the cross section by 
around a factor of 5. 



3.5.2 Total cross section 

We now turn our attention to the comparison of our theoretical predictions to the 
measured total cross section. The analysis we perform includes the processes in which 
the exchange particles are either both pomerons or both reggeons. We do not include 
the case where one is a reggeon and the other is a pomeron, nor do we include in- 
terference contributions. Whilst the latter may be small, the former will not be if 
the pure reggeon contribution is not negligible. This limitation arises since pomeron- 
reggeon interactions are not yet included in POMWIG. The results are presented in 
table 3.3 for E^^ > 7 GeV and in table 3.4 for E^^ > 10 GeV at = 1.8 TeV 
(for the data, the first error is statistical and the second is systematic). The parton 
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Parton level [nb] 


Hadron level [nb] 


Hadron level + oc [nb] 


CDF Result 






43.6 ± 4.4 ± 21.6 


P fit 4 


6.4 


2.2 


7.3 


P fit 5 


859.3 


190.6 


661.8 


P fit 6 


886.7 


230.8 


702.1 


P 


184.7 


13.2 


244.1 


P+P fit 4 


191.1 


15.4 


251.4 


P+P fit 5 


1044.0 


203.8 


905.9 


P+P fit 6 


1071.4 


244.0 


946.2 



Table 3.3: The total DPE cross-sections for a i?™"^ cut of 7 GeV at ^/s = 1.8 TeV. 
The contributions from pomeron and reggeon exchange are shown separately. 





Parton level [nb] 


Hadron level [nb] 


Hadron level + oc [nb] 


CDF Result 






3.4 ± 1.0 ± 2.0 


P fit 4 


1.2 


0.4 


1.4 


P fit 5 


138.6 


31.2 


110.9 


P fit 6 


174.7 


42.9 


141.7 


P 


13.2 





13.2 


P+P fit 4 


14.4 


0.4 


14.6 


P+P fit 5 


151.8 


31.2 


124.1 


P+P fit 6 


187.9 


42.9 


154.9 



Table 3.4: The total DPE cross-sections for a Ef" cut of 10 GeV at ^ = 1.8 TeV. 
The contributions from pomeron and reggeon exchange are shown separately. 

level result is found by considering the two partons produced from the hard event as 
final state jets, the hadron level result is found by applying the cone jet algorithm to 
the produced hadrons, and the out-of-cone corrections result is found by applying the 
CDF out-of-cone corrections to the hadron level result ^. 

Using fit 5, the overall cross section that we predict for a E^™ cut of 7 GeV 
is 203.8 nb at the hadron level. When we apply out-of-cone corrections to this hadron 
level result we obtain a total cross section of 905.9 nb; this is close to the parton 

^We have noticed that there is a considerable difference between our results for the 
Tevatron Run I, at s/s = 1.8 TeV, and results for the Tevatron Run II, at i/i = 1.96 TeV; 
for example the pomeron only, fit 5 total cross section for the Tevatron Run II is 1036.6 nb 
at the parton level and 254.4 nb at the hadron level, compared to 859.3 nb at the parton 
level and 190.6 nb at the hadron level from table 3.3. 
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level result and indicates that the out-of-cone corrections, as derived from Run I data, 
are approximately achieving what they were intended to do and sufficient amount of 
transverse energy is being added to the jets to "correct them" back to parton level. 
As we discussed in section 3.2, it is the out-of-cone corrected data that we need to 
compare to the CDF experimental data. 

This (out-of-cone corrected) predicted total cross section of 905.9 nb is in excess 
of the experimental value of 43.6 nb. A similar excess is present with a cut 
of 10 GeV. However, we can match our results to the data if we assume an overall 
multiplicative gap survival probability of around 5%^. The large reggeon contribution 
implies a non-negligible pomeron-reggeon contribution and naively estimating this as 
twice the geometric mean of the pomeron-pomeron and reggeon-reggeon contributions 
would push the gap survival factor down to around 3%. Given that the systematic 
error on the CDF cross sections is high, that the uncertainty in our knowledge of 
the gluon density directly affects the normahsation of the cross section, and that the 
size of the reggeon contribution is also uncertain it is not possible to make a more 
precise statement about gap survival. In any case the value we obtain agrees well 
with the expectations of [31,32]. Both fits 5 and 6 can describe the data in this 
way, although measurements of diffractive dijet production at HERA suggest that 
fit 5 is favoured [47]. The ratio of the fit 5 to the fit 4 cross sections is of the order 
of 100, which we can understand from the ratio of the gluon densities, illustrated in 
figure 3.7. This ratio is 10, which becomes 100 when we consider the gluons in both 
pomerons. Note that the relative size of the reggeon contribution compared to the 
pomeron contribution is not small. 

The suppression of the total cross section, relative to the naive parton level result, 

resulting from the parton shower phase in POMWIG has been looked at in [35] . Not 

surprisingly, the parton shower phase of POMWIG is responsible for a large part of 

the suppression relative to the naive parton level prediction. 

^ These estimates are extracted from the = 7 GeV predictions but approximately 

apply to all cuts. 
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3.5.3 Event distributions 

In figures 3.9 and 3.10 we show distributions in number of events of the mean jet 
transverse energy, E^, the mean jet rapidity, rj*, the azimuthal separation of the 
jets, A0, and the dijet mass fraction, Rjj: 

~ ^A- (3-21) 

The sum in the numerator is over all particles in the dijets. Some of these distributions 
have also been examined in [33]. In figure 3.9 we compare the data to results at 
the parton and hadron level, and we show the reggeon contribution separately^. In 
figure 3.10 we show results at the hadron level for the three different HI pomeron 
parton density functions. 

We urge caution when comparing the data and theory as is done in figures 3.9 
and 3.10 since the data are not corrected for detector effects''. The existence of the long 
tail to low angles in the A0 distribution illustrates the dangers: there is no possibility 
to produce such a long tail in a hadron level Monte Carlo simulation. This tail is 
caused by misidentifying the two highest Et jets in a three-jet event and defining A0 
as the azimuthal separation between two neighbouring jets; this cannot occur in a 
Monte Carlo simulation that has not been corrected for detector effects. Hence we are 
unable to draw any strong conclusions until the corrected data become available. 

3.6 Conclusions 

The work of this chapter has been focused on the so-called double pomeron exchange 
(DPE) process recently measured at the Tevatron [26]. This process is characterised 
by a double rapidity gap separating the intact, diffracted interacting hadrons from a 
central dijet system. We have extended previous calculations [22,29] by including the 
effects of parton showering and hadronisation and found that they lead to a suppression 

^All curves except the reggeon are area normalised to unity. The reggeon is normalised 
relative to the total. 

^Primarily because of the low of the CDF jets. 
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Figure 3.9: Comparison of theoretical predictions at the parton and the hadron level. 
Also shown is the separate contribution from the reggeon (normalised relative to the 
total) . 
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Figure 3.10: Comparison of theoretical predictions for the different HI fits to the 
pomeron parton density functions at the hadron level. 
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of the cross-section relative to the naive parton level by a factor of around 5. We also 
found that, in the kinematic region probed by the Tcvatron, the effect of non-diffractive 
(reggeon) exchange is probably important. The out-of-cone corrections used by CDF 
in their experimental analyses have been included in our cross section predictions 
and we found that, if we use a gap survival probability [30] of around 5%, we can 
describe the data in a natural way. At the present time the issue of gap survival, 
which embodies a violation of diffractive factorisation, is not very well understood. 
However, it is encouraging that we are able to describe the data with a gap survival 
probability which is consistent with previous theoretical estimates. 

The Tevatron Run II has produced around 10,000 DPE events, which are in the 
process of being analysed. Once this has occurred, and theoretical studies hke this 
one have been carried out, we will be in a far stronger position to understand double 
diffractive processes. 

Finally, we can claim that one can use HERA partons at the Tevatron. When 
the parton level calculations of [22, 29] were carried out, the very large theory-to- 
experiment ratio meant that a description of the data was impossible, even with gap 
survival. However, a harder (and more appropriate) pomeron intercept, new parton 
densities in the pomeron and a gap survival factor combine to explain this excess. 
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4.1 Introduction 

In this survey chapter we shall discuss the factorisation of cross sections in specific 
regions of phase space [6]. We shall make a statement of factorisation specific to colour 
exchange processes in QCD [48], and develop the consequences: the resummation of 
the soft and the jet functions [48-52] . The final form of the cross section is derived, 
with the leading logarithms of the jet function and the next-to-leading logarithms 
of the soft function resummed, the latter being in terms of so-called soft anomalous 
dimension matrices [48,53,54]. The application of these techniques to rapidity gap 
processes is then discussed, where we are interested in soft, wide angle radiation [53] 
into a restricted angular region. In these applications we are interested in colour evolu- 
tion and not in the resummation of the colour-diagonal jets; hence the soft logarithms 
become the leading logarithms. In this chapter we refer the reader to the literature 
for the proofs of the factorisation theorems we use [6] , and our aim is a survey of the 
consequences of factorisation and the phenomenological applications to rapidity gap 
processes. 
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4.2 Factorisation 

In this section we will describe the factorisation properties of the QCD cross sections 
we are interested in. The resummation formalism developed by Collins, Soper and 
Sterman [6,7,48,52] (known as CSS) depends on the factorisation properties obeyed 
by a cross section, with a hard scale Q, in a particular limit of its final state phase 
space. In such regions of phase space, partons which are off-shell by order are 
described by a hard scattering function H, and its complex conjugate H*. On-shell 
particles with momenta ~ Q fall into "jets" of coUinear particles, which will have the 
interpretation of hadronic jets if they lie in the final state and as parton densities if 
they lie in the initial state of a hadron-initiated process. We assume that the final state 
of this process is in the elastic limit and all finite energy particles are concentrated 
in the jet functions; hence the factorisation is valid at an "edge" of phase space for 
a given observable. This part of phase space, where there is just enough energy to 
produce the final state jets and very little else, is known as the threshold region^. 
In addition, the cross section involves the emission of soft particles, represented by a 
function S. In this work we are interested in hadron-initiated processes to jets, where 
it has been observed that there is no unique way of defining colour exchange in a 
finite amount of time. This is due to the fact that even very soft gluons carry colour 
and we can exchange an arbitrary number of these long-wavelength gluons in the 
period of time that is just after the short distance (short time) hard event. Therefore 
we expect the functions from which we construct the cross section to be written as 
matrices in the space of possible colour fiows - specifically the hard and soft functions, 
as the jets themselves contain only coUinear partons and are incoherent to the colour 
exchange. We can understand the need for a matrix structure by noticing that as the 
factorisation scale changes, gluons that are part of the hard function may need to be 
moved to the soft function (or vice versa), changing the colour content of both. The 
matrix structure of the soft function will be discussed further in section 4.7, where we 
discuss the construction of eikonal cross sections. 

^The resulting resummation formalism, valid in this region of phase space, is consequently 
known as threshold resummation. 
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The elastic hmit of the cross section is isolated by introducing an appropriate 
weight [7] for each final state. The weight, which is a linear function of all final 
state particles if it is possible to resum the observable, determines the contribution of 
each state to the cross section and vanishes in the elastic limit. In effect the weight 
parameterises the distance in phase space from the elastic limit. We shall assume that 
the contributions to the weight of particles in the soft function, cus, and from particles 
in the jet functions, Ui{i = 1 . . . n), are additive. 



0J« 



(4.1) 



for a n jet cross section, and the corresponding statement of factorisation near the 
elastic limit is written as a convolution over the weights [48], 



cr[uj 



H 



(0) / Pi 



n 



X 



X 



(0) f Pi-^i ^iQ 



fx 



where we have suppressed the functional dependence on as and i — 1 . . .n. Note 
that the arguments of H^^^ and S^^^ depend on all possible variables with i subscripts. 
For example, for n = 2 we would have H^^\pi/ ix,p2/ ^,^1,^2)- This expression is our 
starting point and is the cross section for a hard event H, described by Hjj, with n 
light-like jets of particles, with momenta pi. The factorisation scale is denoted /i. 
The kinematic details are determined by the precise process under consideration; for 
example in e+e" annhilation {n — 2), each pi labels an outgoing jet and we define 



(4.3) 



The jets are described by the functions Jj. The indices / and J label the colour 
structure of the hard scattering and each matrix has two indices - one for the amplitude 

and another for its complex conjugate. The vectors are used to define a jet in terms 
of a matrix element (just like a gauge-fixing vector) and hence are arbitrary, and the 
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vectors Vi define the direction of the appropriate jet. The function Sjj describes the 
dynamics of partons with energy less than the soft scale uJsQ\ we will consider its 
construction later in the chapter but it is enough to say that the soft dynamics of 
this function exactly match those of the full theory. Hence resummation of the soft 
function is equivalent to the resummation of the soft dynamics of the full cross section. 
We have denoted bare, unrenormalised quantities by a (0) superscript, and deal with 
the issue of renormalisation in section 4.3. The proof of the factorisation statement 
follows standard proofs of factorisation [6], and shall not be discussed in this thesis. In 
this chapter we are interested in the consequences of factorisation, and so will consider 
a simplified cross section with only 1 jet. Therefore = and we write 

= Hf] ^^^^ sf> s(. - .J - 

\ll J J UJj UJs \ 1^ 1^ J Va* / 

(4.4) 

where we have suppressed the functional dependence on ctg. The convolution in weight 
unravels after the transformation, 

POO 

a{N) = / dwe-^^'aicu), (4.5) 
Jo 

where the moment variable is conjugate to cu, and the leading behaviour in the 
small CO limit is determined by the large behaviour. We therefore obtain, in moment 
space, 

.W = <'(^,,)i<«.(^,-|)^-(-|,.«), (4.6) 

where we define the transform of the soft function by 

with a similar definition for the transform of the jet function. Noting the kinematical 
relation p'^ = Qu^ and writing the soft scale as Q/N = Qs we arrive at 

a(N) = HfJ (« . . s) J<») (9^, ^) Sf! (^, . . e) , (4.8) 

where we have shghtly changed the notation. Equation (4.8) is our statement of 
factorisation, in moment space, for a schematic 1 jet cross section. 



(4.7) 
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4.3 Renormalisation group properties 

In this section we demonstrate the renormahsation of the contributing functions [7,48- 
50] appearing in the factorised cross section. We start from the factorised expression 
discussed in the last section, written with n jets and in matrix notation, 

n 

a = Tr(^(°)£(°))n</f , (4.9) 

i 

in terms of the bare quantities. For the remaining part of this chapter, we shall drop 
the "tilde" notation for quantities in moment space; it will be clear from the context 
what is meant. We assume that Jj renormalises multiplicatively 

Jf^ ^ Z,{,,)U,,) (4.10) 

where Zi{^) is an UV renormalisation constant, Jj(/i) is the finite (renormalised) jet 
function and is an arbitrary renormalisation scale (different to the factorisation scale 
used in the last section). By observing that the bare jet function must be independent 
of the renormalisation scale, 

.|4°'^0. (4.11) 
we arrive at the renormalisation group equation for Jj, 

= -liJi- (4.13) 

In the last equality we have defined the (so called) jet anomalous dimension. 

The renormalisation of comes from the UV limit of the virtual corrections. We 
assume that the amplitude and its complex conjugate renormalise multiplicatively. 

Again, is a UV renormalisation constant, ^ is the finite (renormalised) soft function 
and II is an arbitrary renormalisation scale. Exploiting the invariance of the bare soft 
function against changes in /x, 

= 0. (4.15) 
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means we arrive at the RGE for the soft function 

S-S 



z 



a 11= 



(4.16) 
(4.17) 



In the last equahty we have defined the soft anomalous dimensions. The scale- 
independence of a now allows us to derive the renormalisation properties of H^. Insert- 
ing the renormalised expressions for the bare soft and jet functions into the factorised 
expression for the cross section gives 



^ ' i 

= Tr(^£)n'^^' 

i 

which implies that H^^^ renormalises multiplicatively, 



(4.18) 
(4.19) 



(4.20) 



We can now write down the hard RGE, 



(4.21) 



by noting that the bare hard function is renormalisation scale independent. 

The discussion so far has been for an observable in the threshold region of its final 
state phase space, where there is no extra energy available for gluon emission. This 
restriction that the gluon emission must be very soft means that this factorisation 
is also valid for rapidity gap observables. In such cross sections the gluon emission 
into a measured gap region is restricted, whilst being fully inclusive elsewhere, so the 
soft function depends on radiation into the gap fl, at scale Qs- To be more general, 
the soft function also depends on radiation out of the gap, the region Q, at scale M. 
Therefore we can write ^ = S,{Qs/ A*, M/ jj) and trade the ji derivative for Qg and M 
derivatives to obtain 



(4.22) 
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In [55] the M derivative has been interpreted to produce non-global logarithmic en- 
hancements in 2, starting at a^\og'^(Qs). Configurations producing these logarithmi- 
cally enhanced terms arise from radiation being emitted into Q from secondary gluons 
in Cl and as a result resolve the colour structure of the secondary gluon. Therefore such 
contributions cannot be absorbed in an anomalous dimension matrix for 2 — > 2 scat- 
tering in a simple way, and we drop such non-global terms from our soft RGE. We 
will include such effects (non-global logarithms) in a new way in chapter 5. Therefore 
our soft RGE (correct to primary leading log level) is 

- Pi9.)^ ]a = -^^&- ^r.- (4.23) 



' dQs dg. 

We will solve the soft and jet function RGEs in the next section. 

4.4 Resummation 

In this section we discuss how to use the factorised expression to produce resummed 
expressions for the various contributing functions in the cross section. The resumma- 
tion of the jet and soft functions parallels that of [48]. 

4.4.1 Resummation of the jet functions 

The expression for the factorised 1 jet cross section (in moment space) , introduced in 
section 4.2, is 

a{Q, Qs, as) = Hjj u ■ ^, asin)^ J (^-j^^ ^> «5()")^ Sji v ■ ^, ^^(/x)^ . 

(4.24) 

In this equation, we have restored the ag dependence of all the contributing functions. 
We will now define the scaled gauge parameters, r, such that 

r^u-^^p-^^Qr. (4.25) 

This factorisation cannot depend on the exact choice of the quantity r, 

r^a{Q,Qs,as)^0, (4.26) 



A A T~i 
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and any changes in the function J must be compensated by changes in the functions 

S and H. Therefore (suppressing function arguments) 

T^imJ-'rlj = -Tr {(r|:S£} " Tr {^(r|:£)} • (4.27) 
If we assume that this equation is true element-by-element, 

we obtain (restoring the function arguments) 

~ or \ II II ) or \ii 



d 

— I r 



■a{j^,r,a,{ii)^^a-'- (4-29) 



dr'' 

The log-derivative of J may depend on the hard scale Q or the soft scale Qs- This 
variation can be split into two functions: one function which depends on the variables 
held in common by J and namely Qr/fi and a.sifJ'), and one function which depends 
on the variables held in common by J and S^, namely Qs/a* and as{ii)- Therefore we 
can write 

1 J-V|: J(^, ^, a.iii)) = G «,(/.)) 1+k{^, «,(/.)) 1, (4.30) 



where we define 



H-'r^H(^^,r,as{ii)^ = -G (^^, «,(/.)) 1, (4.31) 



^ n ( „ _ / \ \ \ n— 1 ( Qs 



r^£(^^,r,a,(/i)JJ^-i = -K \^,a,{ii) j 1. (4.32) 

The auxiliary function G depends on the hard scale Q and contains all the short- 
distance physics; conversely the function K depends on the soft scale Qs and contains 
all the long-distance physics. Both of these functions depend on the gauge parameter r. 
Applying the RG-invariance operator /i(9^ to this equation, swapping the order of the 
differential operators, using the RGE of J and noting that 7j is r independent implies 
that the combination G -\- K is RG-invariant, 



d 



G [ — , as{ii) \ + K ( — , asin) 



0. (4.33) 
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Separating variables we obtain 



/^l-cf— ) = +7/f (4.34) 



5/i \ fX 

( — ,as{fx) ) = --fxiasifx)), (4.35) 



dfjL \ IjL 

where the separation constants (Sudakov anomalous dimensions) only depend on the 
common variable as{fi). To solve these equations, we postulate the solution (which 
can be checked by differentiating) 

K asii^)) =K{ ) + I ^^^{asifJ.')). (4.36) 
\ A* / J// A* 

Now we set the upper hmit of the integral to the scale Qs and choose the constant 
K{ ) (which denotes a //-independent K) to satisfy the RG equation. Doing this we 
obtain 

K(^,asi^^)) = K{l,as{Qs))+ r^lKMtJi')). (4.37) 
G f ^, a,{^i)] = G(l, asiQr)) + T ^7^(a,(/.')). (4.38) 

\ J Jqt 

In order to combine these two expressions, we shift the evaluation on the coupling in 
equation (4.37) from Qg to Qr, 



K{l,as{Qr)) = K{l,a,{Q,))+ da, ^ ' (4.39) 

= X(l, as{Qs)) + r ^P{as{n'))^^^P^. (4.40) 
We have used the QCD /3-function in the last line, 

du' 

as^as{n')^^as^ (5{as) — , (4.41) 
where I3[as) — ^'das/d^' . Combining the K and G functions, we now obtain 

G{9L^as{^i)^ + K(^,as{ii)^=G{l,as{Qr)) + K{l,a,{Qr)) (4.42) 

^ ^ {lK{a,{^')) + /5(a.(/.0)^^^') . (4.43) 
This expression can be written in the compact form 

G f^, «,(/.)) + K f ^, asiii)] = - r ^A{aM) - A'{as{Qr)), (4.44) 
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where we define 



(4.45) 

A' {as) = G{l,as) + K{l,as). (4.46) 

We now have two evolution equations for J, from RG invariance and the invariance 
of J against changes in the gauge parameters r, 

/x^log J ( — , — ,as(/^) J = -7j(as(/^)), 

The first of these two equations has the solution 

'Qr Q 



(4.47) 
(4.48) 



log J 



(4.49) 



,as{i^)) = log J (^,l,as{Qs)] - [ ^'yj{as{n')), 
which must also obey the second equation. This gives 

r|: log J {^^, 1, asiQs)^ = G (^|^, a,(g,)^ + (1, a,(g,)) , (4.50) 
which has the solution (recall that both G and K depend on the gauge parameter r) 



\ogj{^,l,as{Qs)^ = logJ(l,l,a,(g,)) 



r dr' 

JQs/Q ^ 

= log J (1,1, ^.(g^)) 

r dr' 

+ / -J- 



G (^, «,(//)) +i^(l,a,(g,)) 



or, after a change of variable (r' = /x'/g), 

log J ^^,l,Q;,(gs)^ = logJ(l,l,Q;s(g,)) 



+ 



d^ 



'Q 



' |/l(a.(0)-^'K(/^')) 



(4.51) 

(4.52) 
(4.53) 
(4.54) 



(4.55) 



. (4.56) 



Now we can write down an expression for J which organises both the Qr depend 
and the N dependence (through , 

log J f^,^, «,(/.)) = logJ(l,l,a,(g,))- r^7j(«.(/^')) 



JQ. ^ I Jo. 



^ df 
, -^AiasiO) 
iQs s 



+ ^'(a,(A))i. (4.57) 
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Extending the ^ integration range using the ©-function, 

1= r^Aia,{OM\-0, (4.58) 

jQs ? 

allows us to switch the order of the ^ and A integrals, and obtain 

f^'di,, ,,,, f^" d\ f^' di ,,,,, fQr\ 



Therefore we obtain our final form of the resummed expression for J, or in general for 
the i*^ jet, of 

log Jj ( — , — ,as(/^) ) = log Ji (1, l,a5(Q3)) 



-7.(«s(A)) (4.60) 
I {^(a.(0) log + A'(a.(0)} . 



This can be generalised to an arbitrary number of jets by repeating the reasoning for 
each jet. Doing so we obtain 

nJ.(^,^,«.(/^)) =nJ.(l,l,«.(Q.))e-'5^ (4.61) 
where we define the jet Sudakov exponent by 

Sj, = T (7.(«.(A))) + £ I {A(a.(0) log + ^'(«.(e))} ■ (4.62) 

The resummation is controlled by the exponents 74,^4' and the set 7^. The latter are 
known as anomalous dimensions, which are expressed as series in as- This completes 
the resummation, in moment space, of the jet function. 



4.4.2 Resummation of the soft function 

The RGE for the soft function at scale ji is (suppressing ag dependence) 

/^|^£(a^) = -^m - m r^. (4.63) 

To solve this equation, we need to use the path-ordered exponential, which allows us 
to expand exponentials of non-commuting objects (in this case the non-commuting 
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objects are matrices) 

Pexp f r *^r,(,')) = 1 + /" *ir.(.o + r itir.(M. 



■ ... + 



+ ^ . . . ^r,(/.„) . . . r^(/.o . (4.64) 

We have introduced a new scale /Iq {< /i). In this definition of the path-ordered 
exponential, we arrange the matrices in order of increasing scale from left to right 
(increasing /i). Therefore at a new order of the expansion, we gain a new inner 
integration of a new variable from /in-i to /Iq, and place the new rj^iin) on the 
left of the existing matrices. The anti-path ordered exponential, denoted P, is similar, 
except that we arrange the matrices in reverse order. 



+ ... + 



+ r ^ . . . r ^r^\n,) . . . r,t(;,„). (4.65) 

The effect of applying jid^ to a path-ordered object can be seen by inspection to be 

fid^Fexp(^J ^r^(/^')^ = -Pexp^y ^r,(/i')^ r^(/i), 

/xa,Pexp^'^°^rt(;x')) = -rt(;x)Pexp^'°^r,t(^'))- (4.66) 
Therefore it is easy to see that we can solve the soft RGE by 

£(/.) = P exp y £l(/^')) exp ■ (4-67) 

We have picked the scale of £ on the right-hand side so that the equation is true 
when ji — jiQ. Note that if Ts{ij) commutes with Ts{iJi') the path-ordering becomes 
irrelevant and we recover the usual exponential form. For the soft function defined in 
this chapter, ^ is a function of the soft scale ratio Qs//i, and so the RGE is given by 

4fi(^)--£i£(f)-fi(^)E.. (4.68) 

This equation is solved by 

where we have picked the arbitrary scale i^o — Qs- 
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4.4.3 The resummed expression 

We can now present (in matrix notation) the resummed cross section in moment space 
for our 1 jet cross section. We have chosen ji — Q throughout, and so the hard function 
becomes H^iQ/n) — ^ ^(1) and the integrals in the path ordered exponentials of the 
soft function and the jet Sudakov exponent are now between Q and Qs- We have 
made the choice of r = 1 for the gauge parameter. We obtain 

aiQ, Qs, a,) = Tr |^(1)P exp (^j'^^ ^^w) ^ exp (^f^ ^El.w) } (l' ^) ' 

(4.70) 

where we have not written the full expression for the resummed jet function, which 
can be read off from equation (4.61). 



4.5 Diagonalisation 

In this section we develop the solution, equation (4.69), of the soft RGE. We do this 
by choosing a basis in which the I\ matrices are diagonal [53,55], and hence are given 
by the eigenvalues of the matrix , 

(f ,) = X^^^5,^, (4.71) 

where denotes the eigenvalues of F^. We will denote quantities in the diagonal basis 
with tildes. If e^"^^ is the J^^ eigenvector of Tg (with eigenvalue A'^'^^) then the matrix 
whose columns are the eigenvectors, 

Rij^eS'\ (4.72) 

will diagonalise F^. We therefore write the anomalous dimension matrices in the new 
diagonal basis using a similarity transform, 

(r;(Ar7, n)) = A^'^)^^;, = R-i (r.(Ar7, fi)),^ Rjp. (4.73) 

To do this we need to explicitly expand the anomalous dimension matrices to one-loop 
(work at leading order in a^), 

r,(A) = ^^a. (4.74) 
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where ^ is a constant (i.e. scale independent) matrix. This imphes that different 
scales will commute with each other (all scales have the same constant matrix Q 
inside of the path-ordering, and the matrix can come outside of the integral on the 
left hand side. Hence 

Pe.(/;^a.)^e.(a/;f^). (...) 

Now by inserting RR~^ = / as required, we rewrite the path-ordered exponential as 



(^R-'Pexp(^J^^^UX)^R^ =5,.exp(^|''^^A(^)(A)), (4.76) 



yr^Wj^j =5jjexp^ 
where the eigenvalues have been expanded to one-loop. 



^ ^A(^)'\ (4.77) 

TT 



and we have moved the constant piece of the eigenvalues back inside the integral 
over A. Transforming the hard and soft matrices to the new basis, 

I = R^S=E, (4.78) 
H = R-^S^R-^\ (4.79) 

allows us to write down the resummed expression for the cross section (with i jets) in 
the diagonal basis, 

a{Q, Qs, as) = Tr |l(l)|(l) exp ^ (A*(^) + A«)) | H J^ (l, ^) • (4-80) 

By introducing the following combination of the eigenvalues, 

^iu) ^^(i),i*^^iJ),i^ (4.81) 

using equation (4.77), we write the resummed expression as 

a{Q, Q., a.) = IV {l(l)|(l) exp (^E^^'^ £^ } E (l' f ) ■ (4-82) 

Inserting the one-loop definition of the running coupling. 
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where /3o — 11 — ^Uf and A ~ 0.2 GeV, allows us to rewrite the soft Sudakov expo- 
nential as 

We have, therefore, written a resummed soft function in moment space in which all 
leading logarithms of Qs/Q have been rearranged (or resummed) in terms of the 
eigenvalues of the soft anomalous dimensions, which are matrices in the space of 
possible colour flows of the system. Prom the factorisation arguments earher in this 
chapter, the resummed soft function contains the same soft logarithms to all-orders 
as the full cross section. 



4.6 Resummation and colour flow in rapidity gap 
processes 

The CoUins-Soper-Sterman resummation formalism discussed in this chapter takes 
place in moment space, and an inverse Mellin transform is required to undo the 
transformation. The two main ways to do this are a fixed order expansion of the 
exponents [51], where the inverse Mellin transforms are trivial, and by a numerical 
inversion prescription [56] . In the latter case the cross section depends on an integral 
in moment space over all scales, including the Landau pole, and a prescription must 
be used to avoid this region. The disadvantage of this procedure is that the resummed 
cross section becomes prescription dependent. For a discussion and further references 
on the inversion procedure see [56,57]. 

In this thesis we are interested in the dynamics of wide-angle soft radiation into 
a restricted angular region Q and the role of colour flow in the hard scattering. The 
colour-diagonal jet functions are incoherent to the soft function colour structure and 
only include the dynamics of coUinear soft radiation. Therefore in the study of rapidity 
gap processes, we only consider the soft function dynamics and the interjet logarithms, 
which arise from soft, wide-angle emission are the leading logarithms. In general, the 
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differential distribution in tfie gap energy Qg is given by 

n>l 0<m<n-l ^ / + 

where we see logarithmic enhancements oiQ/Qs, Qs < Q, at each order in the pertur- 
bation series. We have only shown the leading behaviour for Qs <^ Q (there are other 
non-singular terms at each order in perturbation theory) and the coefficients Cnm are 
numbers. The terms in the plus distribution diverge as Qs — > 0, singular but integrable 
in this limit. Here we recall the definition of the plus distribution for the function g{z), 



1-z 



9{z) \ f \ > 9{z) 



where 



. Qs 

z — 1 — — . 

Q 



(4.86) 



(4.87) 



This definition of the plus distribution can be expressed as an integral with some 
smooth function f{x), 



dx f{x) 



9[x) 
1-x 



dxlf{x)-f{l) 



9{x) 
1-x 



/(I) / dx 



9[x) 
1-x 



(4i 



The case of m n — 1 is the leading logarithmic set. The corresponding integrated 
distribution in momentum space is 



^Q,Qra..)= --—dQs = l + y2as{Qr V Cnmi^rlog 

•^0 n>l 0<m<n-l \m + i 



Q 

(4.89) 

We can take the Mellin moment of the differential cross section, equation (4.85), using 



dz z 



N-l 



log™(l - z) 



1-z 

which generates terms like 



^log'"+'^ + 0(log"^A^), 
m + 1 N 



(4.90) 



a{N) = l + J2as{Qr Yl C^m (^—^og^^' l/N + 0{log"' N)Y (4.91) 

n>l 0<m<n-l ' 

where we only show the leading logarithms in moment space. Note that, in the limit 
of large the Mellin transform used in this section is equivalent to the Laplace 
transform used in section 4.2. The soft resummation procedure organises terms like 
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this (at LL) and hence amounts to a knowledge of the LL coefficients to all-orders. 
Comparing this with the integrated distribution in momentum space, equation (4.89), 
we see that the coefficients of the LL term in momentum space correspond to the soft 
resummed terms in Mellin space. Therefore, to LL level, our resummation of the soft 
function in Mellin space corresponds to a resummation of all leading soft logarithms 
in momentum space. Therefore we can write down the factorised cross section for the 
production of 2 jets at fixed rapidity interval as 

= Hil{s, i, II, as)SLi{Qs, fJ', as), (4.92) 

where S^j is a soft matrix in momentum space, s and t are the Mandelstam variables 
and we write the soft scale as Qs- The resummation of the soft function in moment 
space presented in this chapter therefore allows logarithms of Qs/Q to be resummed 
in da/dArj, the cross section in momentum space. 



4.7 Eikonal cross sections 

In this section we will briefly describe the interpretation of the soft function, which 
appears in the factorised expression of the cross section. This topic has been studied 
in [50,53]. 

This function, which describes the emission of soft radiation in the cross section, 
is constructed by representing the initial and flnal state partons by Wilson lines. 
The soft radiation pattern of these Wilson fine constructions, so-called eikonal cross 
sections, exactly matches that of the full theory. We represent the initial and final 
state particles by an ordered exponential of the gluon field, 

(Ai, A2; x)=F exp (^-ig J ' dl (3 ■ A^^^ {lf3 + x)^ , (4.93) 

where the integration is a straight line path in the direction of the particle 4- velocity p. 
The vector potential A^^^^ is in the appropriate representation to correctly describe a 
parton of fiavour /. Such a Wilson, or eikonal, line is a special case of the more general 
Wilson line which appears many times in quantum field theory, where the path can be 
any curve C. The path ordering symbol is needed to define the Wilson line in terms of 
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a power expansion of the exponential with the SU(3) matrices ordered so that higher 
values of / stand to the left. The wide angle, soft radiation pattern which is obtained 
from replacing each of the partons in the hard event by Wilson lines exactly matches 
the wide, angle soft radiation pattern of the full event. Therefore we construct an 
operator consisting of products of Wilson lines, tied together by their colour indices. 
As a transparent example, consider two Wilson lines representing partons in a Drell- 
Yan process. 



(4.94) 



where the two Wilson lines extend from negative infinity and meet at the vertex x. 
The initial colour content of the hues is represented by the indices 6, and the colour 
exchange at the vertex (in this case singlet) is represented by the tensor 5aia2- This ver- 
tex, when dressed with higher order gluons, contains UV divergences; this is also true 
of the one-loop external leg self-interaction diagram. The sum of these is renormalised 
in the following manner. 



W(^^)W(x) = Zj^'^) (a„ /?„/?,) iy(^^)(x). 



(4.95) 



and the invariance of the renormalised operator under changes in the renormalisation 
scale implies that the renormalised quantity obeys the standard RGE, 



where we define the colour singlet anomalous dimension, 
A one-loop calculation gives 

'-2/?i-/?2 



1 d 
Ziw dji 



(4.96) 



(4.97) 



Tw = — Cp 

TT 



log 



- 1 



(4.98) 

have the interpre- 



where /3j is the 4- velocity of eikonal line i. The divergences 
tation of collinear singularities. 

We can build eikonal operators for 4-jet processes by tying together four Wilson 
hues, two for the initial state partons and two for the final state partons. The colour 
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structure is now more complicated, and the eikonal operator is constructed with a 
colour structure c connecting the colour indices of the eikonal lines. We now choose 
to span the colour space of the process with the set of colour tensors c/, and expand c 
over this basis. Therefore the set of process-dependent tensors ci link the four Wilson 
line colour indices together. Hence we define the eikonal operator, 



di 



d2d\ jdgdj^ 



(4.99) 



for the process fA{lA,rA) + fsilB^rB) + 12(^2,^2), with a corresponding 

notation, I and p, for the particle 4-velocities. The index r denotes the colour index 
of the relevant parton. The c/'s form a basis in colour space; an example would be 
singlet and octet exchange in the t-channel. The initial state Wilson lines extend from 
negative infinity to the vertex at x, and the final state Wilson lines extend from the 
vertex at x to positive infinity. This operator should reproduce the same wide angle, 
soft radiation as the full theory. In terms of the operators Wf we can define an 
eikonal soft function for fixed energy Qgap flowing into Q, 



(4.100) 



and the corresponding integrated soft function. 



(0|T [(W^f (0))y IO(e|T [wP{Q)iH}] |0), (4.101) 



which are defined as matrices in the space of the colour tensors cj. In these equa- 
tions is the energy flowing into the gap region of the intermediate state ^ and we 
sum over all intermediate states contributing to the cross section. If we demand the 
sum is over all intermediate states of maximum energy Qg into the gap region VL (i.e. 
away from the eikonal lines), the eikonal cross section is free from potential coUinear 
singularities. 
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If we define the soft function using the differential form, equation (4.100), in mo- 
ment space then, following the arguments in section 4.6, the resummation of this soft 
function is equivalent to the resummation of the integrated soft function in momentum 
space. Hence we can write the soft function in momentum space as equation (4.101), 
with the associated factorisation statement of equation (4.92) and resum the soft log- 
arithms using the resummation formalism (defined using equation (4.101)) described 
in this chapter. 

If there are no soft gluons, at tree level, then equations (4.100) and (4.101) re- 
duce to colour traces. A set of eikonal Feynman rules have been extracted from the 
definition of the eikonal line [50, 53] to allow the computation of the UV poles of the 
soft function and hence the soft anomalous dimension matrices, which were defined 
in section 4.4. The higher order diagrams which contribute to the soft anomalous 
dimension matrix are known as eikonal diagrams. The full set of eikonal Feynman 
rules are in appendix B. For a review of these ideas, see the references given above. 

4.8 Colour mixing under renormalisation 

The resummation in this chapter is driven by the eigenvalues of the soft anomalous 
dimension matrices, which are expressed as matrices in the space of the possible colour 
flows of the system. The matrices are found from the coefficients of the UV poles of 
higher order contributions to the soft function, which was defined in section 4.7; hence 
the calculation of the anomalous dimension matrices is equivalent to the calculation 
of the renormalisation constants of eikonal diagrams, discussed in section 4.7. Once 
we specify a basis of colour tensors, the higher order gluons which cause the UV diver- 
gences result in the set of colour tensors mixing into each other [55] . In other words, 
we generate a quantum mechanical mixing matrix by "dressing" a colour tensor c/ [49] 
with soft gluons and project the result along the direction of cj. This is just standard 
mixing of operators under quantum corrections. The calculation of the colour mixing 
matrices was discussed in section 2.3 and the full set of matrices appear in appendix F, 
along with the set of colour bases used in this work in appendix D. 



1 
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4.9 Conclusions 

In this chapter we have discussed the factorisation theorems of the cross section used 
in the rest of this thesis. These properties were used to write down resummed soft and 
jet functions in moment space. We then turned our attention to the apphcation of 
these results to rapidity gap processes, and we argued that, as we are only interested in 
wide angle, soft radiation into a gap region, we can neglect the colour-diagonal jets and 
focus on the soft function. The inversion of the Mellin transform is a trivial issue in 
this situation and we are able to write down a factorised cross section, equation (4.92), 
in momentum space where the soft function has the interpretation of an eikonal cross 
section. The renormalisation properties of the eikonal cross section then allow the soft 
logarithms to be resummed. 

We will apply these results to gaps-between-jets processes in chapter 6, but first 
we have to deal with a class of leading logarithms not included in the resummed soft 
function. We deal with this complication, non-global logarithms, in the next chapter. 



Chapter 5 

Non-global logarithms in interjet 
energy flow 

5.1 Introduction 

The identification of the inter- jet energy flow as an infrared safe way to study gaps- 
between-jets processes has produced extensive interest in the last few years. By con- 
sidering such obscrvables we may start to formulate a perturbative approach to the 
description of the cross section, as well as probe the interface between perturbative and 
non-perturbative physics at energy scales ~ 1 GeV. The inter-jet energy distribution 
was calculated by Sterman et al [55, 58, 59] by separating out the primary emission 
Bremsstrahlung component and calculating this quantity to all-orders for a 4 jet sys- 
tem. We discussed this resummation procedure, and its application to rapidity gap 
processes, in chapter 4 and will use the formalism in chapter 6 to calculate detailed 
predictions for energy flow processes. However Dasgupta and Salam have pointed 
out [60, 61] that this procedure does not include the effects of so-called non-global 
logarithms, a set of leading logarithms which were shown to be numerically important 
at the energy scales probed by current colliders. 

Recently the HI and ZEUS collaborations [62, 63] performed improved gaps- 
between-jets analyses in which the entire event is clustered into (possibly soft) jets 
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using the inclusive kt algorithm [64-66]. A gap event is then defined by the total 
minijct energy in the interjet region, rather than the total hadronic energy as in previ- 
ous analyses. The hope was that this would 'clean up' the edges of the gap and make 
this observable less sensitive to hadronic uncertainties. 

In this chapter we show that by demanding the gluonic final state to survive a 
clustering criterion the effect of non-global logarithms is reduced, but they are still 
numerically important, at HERA and the Tevatron. The study presented here is 
particularly interesting in the light of the recent analyses by the HI and ZEUS collab- 
orations. 

The calculation of the Bremsstrahlung component of the energy flow observable 
has resulted from extensive work in the last decade on threshold resummation, for 
examples see [48,53-55,67]. This resummation procedure was the subject of chapter 4. 
It was found that by considering the emission of soft, wide angle gluons by light- like 
quarks one may describe the soft gluon dynamics of hadronic processes by an effective 
theory known as eikonal theory. This subject is now very well developed and the 
resummation of such primary logarithms has recently made contact with high-but- 
fixed order calculations [68]. This formalism now allows detailed calculations of inter- 
jet energy flow as well as illuminating insights into the topology of colour mixing [53]. 

However, in the study of single-hemisphere observables [60] Dasgupta and Salam 
identified an important class of logarithms that were missing from the soft gluon 
calculations. These omitted contributions arise in observables that are sensitive to 
radiation in a restricted region of phase space. In [61], they showed that the omitted 
contributions are of leading logarithmic order for interjet energy flow observables, as 
confirmed in [55] , where they were neglected. 

To see the origin of these logarithms more clearly, consider soft gluonic radiation 
into a patch of phase space Q arising from a 2 jet system, where we restrict the 
total energy of radiation into Q to be less than Qq. The primary logarithms arise 
from gluons that are emitted directly into fl, as shown in figure 5.1, with energy 
vetoed down to a scale Qq, which is denoted by a crossed gluon line in the figure, 
and form the single logarithmic (SL) set (aglog-^^ , n > 1, where Q is the scale 
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Figure 5.1: A veto of primary emission going into Q. The kind of LL terms produced 
by this kind of diagram are understood from resumming independent gluon emission. 




Figure 5.2: A veto of emission into Q from a gluon radiated outside of Q. This 
configuration will produce a non-global logarithm. 



of the jet line. Now consider a gluon being emitted outside of with intermediate 
energy Qi, and then vetoing subsequent emission from this gluon into Q down to 
scale Qfi, as shown in figure 5.2. Integrating Qi up to Q then generates another SL 
set of (^as log ^ j , n > 2, which are formally the same order as the primary emission 
terms. This was studied in detail in [61,69] and in this work we consider how the 
conclusions found are modified phenomenologically when the kt clustering algorithm 
is applied to the final state. This new class of logarithm, which enter at leading 
log (LL) accuracy in energy fiow observables, are known as non-global logarithms 
(NGLs), as they arise in non-global observables, or as secondary logarithms, as they 
are produced by secondary gluonic radiation. 

The organisation of this chapter is as follows. We start in section 5.2 with a 
definition of our observable, and a discussion of how we include NGLs in our calcu- 
lations. Section 5.3 then describes the kt clustering algorithm used in the HI and 
ZEUS analyses, known as the inclusive jet clustering algorithm. We then develop the 
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algorithm, in order to derive a form suitable for use with a 2-jet system. Section 5.4.1 
then describes the order calculation of the effect of NGLs in our 2-jet system, with 
the clustering algorithm included in the calculation. We find that the asymptotic 
suppression factor with clustering is reduced with respect to the non-clustering case. 
We then proceed to calculate the non-global effect to all orders in section 5.4.2 using 
the combination of the large N^. limit and a Monte Carlo algorithm and find that the 
NGLs act as a suppression to the observable, and this suppression is reduced when we 
impose clustering on the final state. In section 5.5 we give an overview of some of the 
current research work in the area of non-global observables and finally we conclude 
with a summary and directions for further work in section 5.6. 

To summarise our conclusions, the all-orders treatment used in this chapter allows 
direct inclusion of the kt algorithm, in exactly the same way as is used experimentally, 
and we find that the clustering process reduces the magnitude of the non-global cor- 
rections to the primary suppression factor. However they are still phenomenologically 
relevant at HERA. The research work in this chapter is pubhshed in [70,71]. 

5.2 Non-global logarithms in a 2-jet system 

Following Dasgupta and Salam [61], the observable we are interested in is the total 
transverse energy Et flowing into a region of phase space Q for an event characterised 
by the hard scale Q, 



We are speciflcally interested in the cases of Q being either a slice in rapidity or a 
patch, bounded in rapidity and azimuthal angle. The quantity we shall calculate is 
called En and is defined to be the probability that E-t is less than some energy scale 



We normalise this equation to the process cross section, for example 2 jet production 
in e'^e~ annihilation, and we shall assume the strong ordering Qq -C Q- The aim of 
this work is to calculate the importance of the non-global contribution to Eq and so it 




(5.1) 



Qn, 




(5.2) 



V r\ TVT 
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is convenient to factorise this expression into a function describing primary emission 
into Q, p(t), and a function describing (secondary) emission into Q from large- angle 
soft gluons outside of S{t), 

Sn(t) =5(t)So,p(t). (5.3) 

This expression is accurate to LL level (it is helpful to think of this as a probabilistic 
expression, and the observable can be written as the probability for the region Q to 
stay empty from primary gluons times the probability for the region fl to stay empty 
from secondary gluons) and the primary emission function Ef^ p is the 2 jet analogue 
of what has been calculated for 4 jet systems. We have denoted the following integral 
of as by t, 

t{Qa,Q) = :r / i^<^sikt), (5.4) 
27r Jq^ kt 

' loJ^M). (5.5) 



where the first equality is exact, the second holds at one loop, the third assumes a 
fixed coupling and /So — (HC'a — 2n/)/(127r). The leading order contribution to S{t) 
comes in at and we shall calculate this for a 2 jet system in the next section but it 
is useful to first consider the primary emission function at first order in a^. If we do 
not restrict the phase space for gluon emission then, order by order in perturbation 
theory, we expect a complete cancellation of real and virtual soft gluon contributions 
to the primary emission form factor. However the requirement of a gap in a restricted 
region of phase space results in this cancellation being spoilt and we are left with an 
integral over the vetoed region. Hence, to order cts, 

= -iC,^ r '-^ [ (5.7) 



27r,/Q^ kt 27r 

where Aq denotes the area in {r), 0) space of the region fl. When the energy scales are 
strongly ordered, the logarithm can become large enough to overcome the smallness 



V r\ TVT 
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of as and it becomes necessary to include terms like y^s log ^ j to all orders. The 
assumption that the primary gluons are emitted independently according to a two- 
particle antenna pattern allows one to exponentiate the one loop answer and run the 
couphng to the scale kt- We can then write T,Q^p{t) to all orders, 

Sn,p(i) = exp i-iCpAnt) . (5.9) 

This equation only includes contributions from independent primary emission. We 
also need to account for secondary gluon emission, which generates NGLs, through 
the function S{t). This has the following expansion in a^, 

S{t) ^l+S2f + S3t^ + ...^l + J2Snt''. (5.10) 

n=2 

Note that S{t) has its first non-trivial term at order because we need at least two 
gluons to cause a NGL to appear. The function S{t) has been studied for a 2 jet system 
without any final state requirements [61] and the purpose of this work is to calculate 
the function, at leading and at all orders, with the requirement of kt clustering on the 
final state. Therefore, before we turn to the details of the calculation of S{t), let us 
discuss the kt clustering algorithm. 



5.3 The kt clustering algorithm 

In this section we shall review the formulation of the kt clustering algorithm which 
we will use in our calculations. We also provide a formulation that can be directly 
applied to our fixed order calculation in the next section. 

The version of the algorithm we use, which is the one used in the HERA analyses, 
is known as the inclusive kt algorithm [64-66] . The main features of importance to the 
present study are: the clustering procedure starts from the particles of lowest relative 
transverse momenta and iteratively merges them to construct pseudoparticles of higher 
transverse momentum; the decision of whether a particular pair of pseudoparticles 
are merged depends on their relative opening angle and their transverse energy (see 
below); despite this, it is possible for soft particles to be 'dragged' through relatively 
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large angles by being merged with harder particles, which are merged with even harder 

particles, and so on. Nevertheless, we do not expect our results to be qualitatively 
different from other infrared safe jet algorithms, such as the improved Legacy Cone 
algorithm of [72]. 

In the formulation of the algorithm, we represent the final state of a process 
by a set of "protojets" i with momenta pf . The algorithm works in an iterative 
way and groups pairs of protojets together to form new ones. The aim is to group 
almost-parallel protojets together so that they are part of the same protojet. Once 
certain criteria are met, a protojet is considered a jet and is not considered further. 
The algorithm therefore needs to specify what it means for two objects to be close 
together and also how to merge two protojets together. Here we follow the HI and 
ZEUS analyses closely and set the radius parameter, R, to unity. The procedure is 

1. For each protojet, define 

di = El, (5.11) 

and for each protojet pair define 

dij = mm{El„ El^)[{rji - rjj)' + (0, - <Pjf]/R'. (5.12) 

2. Find the smallest of all the di and dij, and label it dmin- 

3. If dmin is a dij, merge protojets i and j into the new protojet k with 

ET,k — ET,i + Etj, (5.13) 

and 

Vk = [ET,iVi + ET,jVj]/ET,k: 

<Pk = [ET,i<Pi + ETM/ET,k. (5.14) 

4. If (imin is a di, the corresponding protojet i is not "mergable". Remove it from 
the list of protojets and add it to the list of jets. This procedure continues 
until there are no more protojets. As it proceeds, it produces a list of jets with 
successively larger values of di — E^-. 
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This algorithm is implemented in the package KTCLUS [73], which is used for the 
all-orders calculation later in this work. 

It is necessary to use the full iterative algorithm for experimental analyses and 
for Monte Carlo applications, but when we consider a 2 gluon final state in this work 
we can reduce the algorithm to a convenient analytic form. We start by considering 
a hard jet line at scale Q which radiates a gluon with some transverse energy Et^i in 
some direction This gluon then proceeds to radiate a secondary soft gluon 

with transverse energy Et,2 in some direction (772,02)- We assume that the energies 
of the gluons are strongly ordered. 

El > E2, (5.15) 

which implies the same for the transverse energies, 

Et,i > Et,2. (5.16) 

Applying the clustering algorithm to the two gluons gives, 

di — E^i, 
6,2 — E'^2-1 

dl2 = £;|,2[(?7l-772)' + (01-02)']/i?'. (5.17) 

By considering the strong ordering of the transverse momenta, the two gluons will be 
clustered if dij < d2, so we require 

im - mf + (01 - <t>2f > (5.18) 

for the two gluons to constitute separate jets and not be merged by the algorithm. 
Therefore for a 2-gluon final state to pollute the gap and generate secondary logarithms 
we require that gluon 1 is outside the gap, gluon 2 is inside the gap and that they 
be sufficiently separated in the {r], (p) plane to avoid being merged. The result of 
this section is that the clustering condition manifests itself as a 0-function in our 
calculation, 

Q{{m-mf + 4>l-R'), (5.19) 
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where we have used our freedom to set (f)i — 0. We will use this result in the next 
section. Therefore the clustering procedure, when applied to a 2 gluon system, will 
pull the softer gluon out of the gap and merge it into the harder gluon, provided they 
are sufficiently close in the [r], 0) plane; hence we expect that clustering will reduce 
the impact of NGLs in such systems. 

5.4 Controlling the non-global logarithms with kt 
clustering 

In this section we describe our calculation of the NGLs with clustering for a 2-jet 
system. We contrast our work to the unclustered case and find that the clustering 
procedure reduces the numerical impact of the NGLs. We then use a Monte Carlo 

program in the large Nc limit to find that this conclusion persists to all orders and 
we make numerical estimates of the impact of NGL in energy flow measurements at 
HERA. 

5.4.1 The LO result 

In this section we will calculate the leading order contribution to the non-global loga- 
rithm function S{t). As has already been shown, this has its first non-trivial term at 
order and has the perturbative expansion, 

S{t) ^ 1 + S2t^ + Sst^ + . . . ^ 1 + J2Snt''. (5.20) 

n=2 

Our goal is to calculate the 5*2 piece for Q defined as a slice in rapidity of width Ar] with 
the condition that the topology of the gluon tree satisfies the kt clustering algorithm. 
We shall initially perform the calculation in terms of the polar angles of the emitted 
gluons, and then again in terms of the rapidities of the emitted gluons because the 
latter form of the calculation will turn out to have more convenient properties. A 
precise definition of 5*2 can be given by following [60, 61] and using the 2 body phase 
space of the emitted gluons. If we only keep the leading logarithmic piece of the 
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integrals over the momentum fractions, defined by 



= (5.21) 



we obtain the following expression for S2, 



/ cicos^2^^ [\2dx2 f\idx,Q (x2 - ^) Ws. (5.22) 

Therefore S2 is the coefficient of the logarithm-squared term in the expansion of S{t). 
The limits on the angular integration ensure that gluon 2 enters region fl while gluon 1 
does not. This expression for 5*2 contains the secondary part, 1^5, of the well-known 
squared-matrix element for the energy ordered emission of two gluons, which can be 
derived from standard techniques [74], 

(ab) (Ca (al) , Ca {hi) 



W = AC 



Ca\ {ah) \ 
\^ 2 J {a2){2h))^ 



(al)(16) V 2 (a2)(21) 2 (62)(21) 
= CIWp + CfCaWs, (5.23) 

where the notation {ij) denotes the dot product of the appropriate 4-momenta. This 
expression contains the primary emission piece Wp, proportional to Cf, and the piece 
that interests us, which is the part proportional to CfCa and denoted 1^5. Note 
that the last term is the dipole interference term and is absent in the large Nc limit. 
We shall begin by writing the 4-momenta in terms of the polar angles of the emitted 
gluons, 

ka = |(1, 0,0,-1), (5.24) 
k, = |(1, 0,0,1), (5.25) 
ki — xi ^(1) 0, sin6'i, cos6'i), (5.26) 
k2 — ^2^(1, sin ^2 sin 02, sin ^2 cos 02, cos ^2), (5-27) 
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where the indices a and h refer to quarks, the indices 1 and 2 refer to radiated gluons 

and we have set 0i = 0. If we evaluate the 4-momenta products we get 

{ah) = ^, (5.28) 

(la) = ^!p(i_cosei), (5.29) 

(16) = ^(i + cos^i), (5.30) 

(2a) = ^(l-cose2), (5.31) 

(26) = ^(l + cos^2), (5.32) 

(12) = — ^ — (1 — sin ^1 sin 6*2 cos 02 — cos 6*1 cos 6*2). (5.33) 

Inserting these results into equation (5.23) we obtain 
32 1 / 1 



Ws 



(2^12^2)^ \ (1 +cos^i)(l — cos^2)(l — sin ^1 sin ^2 COS 02 — cos^iCOS^2) 

1 



+ 



(1 — cos^i)(l + cos^2)(l — sin^i sin ^2 cos 02 — cos^i cos ^2) 

2 



(5.34) 



(1 + cos6'i)(l - cos6li)(l + cos6'2)(l - cos 6*2) 

We shall write the angular part of Ws (the part without the overall constant) as Ws- 
Inserting this result into equation (5.22) and then performing the straightforward 
energy fraction integrals (the ©-function ensures that sufficient energy reaches the gap 
region Q), where we have kept the leading logarithmic piece (dominant when Qn/Q 
is small), 

we arrive at an expression for 5*2 

^2 = -CfCa I Ws. (5.36) 

J angles 

Specialising to a slice in rapidity of width Ar], delimited by cos^^ = ±c and centred 
on 77 = we can exploit the symmetry of the ci integral and perform the 0i integral 
to express the integral over angles as 



i an 



r% 

— 2 I acosi / acos2 / — — . (5.37) 



angles J —1 J —c JO 
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Figure 5.3: The integration region for the polar angle integrations of the emitted 
gluons. Gluon 1 is emitted outside of the gap, and subsequently emits gluon 2 into 
the gap. 



The polar angle integration region is shown in figure 5.3. For the case of no clustering 
of the gluons it is possible to perform the azimuthal average using the result, 



2" d(j)2 1 1 



27r 1 — sin 9i sin 6*2 cos 02 — cos 9i cos 02 \ cos 02 — cos 0i \ ' 

1 



(5.38) 



cos 02 — COS 01 ' 

where we have exploited the angular ordering of the gluons in writing the second 
equality. This results in the following, azimuthal- averaged, expression for the matrix 
element, 

with two remaining cos integrations. For this case of no clustering of the gluons, these 
remaining integrations have been performed, resulting in the following expression for 
^2 [61], 



S2 — —4:CfCa 



^2 + (Arj)' - Arjlogie'^^ - 1) - -Li2(e-^^^) - -Li2(l - e^^") 



(5.40) 

which, following the literature, we have expressed in terms of Arj, which is related to c 
by the kinematical relation, 

Ari = In (^^) ■ (5.41) 
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We have used the dilogarithm function, which is defined 

Li,(y) = fML^d,. (5.42) 

We can now use the following identities (with x — exp(— 2A7^) in the first identity 
and X — exp(2A77) in the second identity), 

U^ix) = y-log(a;)log(l-a;)-Li2(l-a;), (5.43) 
Li2(l-a;) + Li2(l-a;-i) = -^log2(a;), (5.44) 



to write this result as 

-S'2 = —ACfCa 



y - ^2(6--")' 



(5.45) 



Note that as A?7 increases, S2 rapidly saturates at its asymptotic value. 



27r2 

hm ^-CpCa—, (5.46) 

Ar?— »oo o 

since we note that Li2(0) = 0. This analytic evaluation of 5*2 is not possible for 
the case of clustered gluons because, as discussed previously, the requirement that a 
gluonic final state is in a configuration that will survive the clustering algorithm can 
be written as a ©-function of all three angular integration variables, 

e{{r^,-mf-rct>l-R'). (5.47) 

In this situation we have to resort to a numerical solution, in which we integrate 
over the three remaining variables. This method has the advantage of being easily 
extendible to any gap geometry of interest. It is possible, however, to obtain a result 
for 5*2 with clustering which is expressible as a ID integral by using the rapidity of the 
emitted gluons. This has the advantage of requiring a simpler method of numerical 
solution. In analogy with equation (5.22), we can express ^'2 in terms of the gluon 
rapidities. 



-CfCa / drji^ I drj2^%- [ X2dx2 [ xidxiQ ( X2 - ) Ws 
fei^n 27r 7fc2€n ^^r lb Jo Jx2 \ Q 

2(5q/ \ 



= -^2^ {7^)+o(log(^)). (5.48) 
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Note that in this equation, xi and X2 are defined as transverse momentum fractions, 

Pt,i^x,^, (5.49) 

and not energy fractions as they are in equation (5.22), which was written in terms 
of polar angles. However, this does not affect the result. The definition of Qq has 
also correspondingly changed. We can relate the phase space in the two coordinate 

systems by 

-rrr — kt dkt dr] d(j). (5.50) 
k^ 

We begin by defining the following 4-momenta, 

ka = |(1, 0,0,-1), (5.51) 

h = §(1,0,0,1), (5.52) 
ki — Pt,i(cosh77i,sin0i,cos0i,sinhr7i), (5.53) 
^2 = Pt,2(cosh7;2,sin02,cos02,sinh772), (5.54) 

where the indices are the same as those defined for the polar angle calculation. Eval- 
uating the 4-momenta products we get 

(5.55) 
(5.56) 
(5.57) 
(5.58) 
(5.59) 

7^2) - cos(0i - ^2)) ■ (5.60) 



{ah) = 


2 ' 


(la) = 


4 exp( 771), 


(16) = 


^ exp(+7?i). 


(2a) = 


4 exp( 772), 


(26) = 


^ exp(+r72), 


(12) = 


— ^ — (cosh(r7i 



Inserting these results into the expression for Wsi equation (5.23) and using our free- 
dom to set the azimuthal angle of the first gluon to zero, 0i = 0, we obtain 

Ws = ( .r"'"' : , - 1) . (6.61) 

Q'^xfx^ \cosh(?7i - 772) - cos 02 / 
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Figure 5.4: The integration region for the rapidity integrations of the emitted gluons. 
Gluon 1 is emitted outside of the gap, and subsequently emits gluon 2 into the gap. 



The transverse momentum fraction integrals are straightforward, as before. 

Specialising to a slice in rapidity of width At] and delimited by \ri\ < Ari/2, where 
boost invariance allows us to centre the gap on = 0, we can exploit the symmetry 
of the rji integral and perform the trivial 0i integral of the first gluon to express the 
final integral over angles as 

/ =2r-*,/' d,,/'* (5.63) 

Jangles J -oo J Jo ^^T 

The rapidity integration region is shown in figure 5.4. For the case of no clustering of 
the gluons it is possible to perform the azimuthal average using the result, 
'■''^ #2 1 1 



/ 

Jo 



(5.64) 



Iq 2tt cosh(r7i - 772) - cos 02 | sinh(r7i - 7^2) | 
This can be elegantly proven by noting that the following contour integration over the 
unit circle, 

TT^^ (^-exp(?72 -77i))(^-exp(r7i -772))' ^^'^^^ 
is equivalent to the integration over the azimuthal angle of the second gluon. To see 

this, expand the cosh(?7i — 772) and cos(02) as exponentials and write z = exp(i02)- 

We get poles at 

z = e^ e-^ (5.66) 
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where y — rji — rj^- We only have to consider the poles which satisfy \z\ < 1, which 

is 2; = exp(— 1|/|), so 

I = 2m i^] lim ^ --exp(-M) 



nj z^eM-\y\) {{z- exp{-\y\)){z - exp(|y|)) 



exp(-||/|) -exp(||/|)' 
^ (5.67) 



sinh(|y|)' 

and write sinh(||/|) = | sinh(|/)|. As for the polar angle case, it is possible to do all the 
integrations for the case of no clustering. However, evaluation of S2 is not possible 
for the case of clustered gluons because, as discussed previously, the requirement that 
a gluonic final state is in a configuration that will survive a clustering algorithm 
appears as a ©-function in our calculation, and the argument of the ©-function is 
very complicated. However we can readily reduce the three-dimensional integral to a 
onc-dimensional integral using standard techniques if we consider the region of phase 
space that is vetoed by the clustering algorithm, denoted 5*2. Doing this we obtain, 

— 32CpCa 

'0 

dri, 



SI = 32CfCa 1 A77)|"2coth(77) X (5.68 

Jo 



' tanh(r^/2) ' ^ 



rf 



where we define 77 = 771—772. This result is proved in appendix C. Therefore the solution 
for 5*2 with the clustering condition imposed is obtained by subtracting S2 from the 
analytic unclustered result. Wc resort to numerical techniques to solve this equation, 
which has the advantage of being easily extendible to any gap geometry of interest. 
Our numerics were done using Monte Carlo integration methods and the routine Vegas. 
Figure 5.5 shows our result for 5*2 as a function of the gap width At; and with the radius 
parameter, R, set to 1. The solid line is S2 without clustering and the dotted line is S2 
with clustering. In [61], where the unclustered case was first obtained, the observed 
saturation was given a simple physical explanation: the dominant contribution to S2 
comes from regions where the two gluons are close together (which means that gluon 1 
is just outside the gap and gluon 2 is just inside the gap.) This occurs when 

771 ~ 772 ~ -A77/2, (5.69) 
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without clustering 

with clustering 

\ , , , 

"01234 

Ar| 

Figure 5.5: 5*2 as a function of At]. The clustered curve is obtained using a kt algorithm 
with a radius parameter of 1.0. 

and the dominant contribution arises from the coUinear region of the matrix element. 
Therefore at large Arj, S2 received no contribution from the centre of the gap, only 
the edges, and the value of —S2/{CpCa) saturates. This is the reason that NGLs 
are often referred to as an 'edge effect', and gluons closest to the boundary of the 
restricted region of phase space make the largest contribution. This effect is illustrated 
in figure 5.6. We note that S2 itself saturates when A77 ~ R. The results we have 
obtained show that when we demand the gluonic final state to survive a clustering 
algorithm, the saturation of S2 observed in [61] is still observed but the saturation 
value is reduced by 70%. In other words the value that —S2/{CfCa) saturates to is 
reduced from 6.57 to 1.81. The reason for this reduction is that we have removed 
gluons from the region of coUinear enhancement, but gluons that are still sufficiently 
separated in the (77, (p) plane to satisfy 

{m-V2f + (l>l>R^ (5.70) 

will survive the clustering process and contribute to the NGLs. Therefore we conclude 
that the saturation of the non-global contribution at fixed order is still seen when we 




Chapter 5. Non- global logarithms in interjet energy flow 



103 



o 



-3-2-10123 

Figure 5.6: Contributions to the leading order NGLs for a patch in rapidity, as a 
function of rapidity and azimuthal angle. Darker shades correspond to a larger con- 
tribution and it can be seen that the largest contribution comes from the edge of the 
patch. This figure is taken from [71]. 



demand clustering of the final state. 

5.4.2 Non-global logarithms to all orders 

The analytic treatment of the 2 gluon system, with or without clustering, is fairly 
straightforward. The extension to the many-gluon case presents considerable prob- 
lems due to the geometry and the colour structure of the 'mexican cactus' gluon 
ensemble. Therefore we have extended our calculation to all orders by working in the 
large limit and by employing numerical methods. The complicated colour algebra 
one obtains is simplified in the large limit, because the squared-matrix element 
for the gluon radiation can be broken down into a series of independent terms, each 
associated with a different colour dipole. The complicated geometry of multi-gluon 
events is handled by use of a Monte Carlo algorithm, in which the gluon cascade is 
built up from successive dipole branching. Therefore we start with a single colour 
dipole, representing the initial quark line, which emits a gluon and splits into two 



^ 3 

2 
1 

-1 
-2 
-3 



f A /~\ i IT 



Chapter 5. Non-global logarithms in interjet energy flow 



104 




0.8 



in 0.7 



0.6 



O.5L 



Slice of width 1.0 - non-global contribution 

Slice of width 1.0 - exponentiation of fixed order result 

Patch of width 2.0 - non-global contribution 

Patch of width 2.0 - exponentiation of fixed order result 



_J 

0.05 



_L_ 
0.1 



_J 

0.15 



0.2 



0.25 



Figure 5.7: The function Sif) for a slice and a patch of phase space, with the condition 
of kt clustering imposed on the gluons. These curves were obtained using a Monte 
Carlo procedure in the large A^c limit. Also shown are the curves for obtained by 
exponentiating the fixed order result. The geometry independence of the t dependence 
indicates that the buffer zone mechanism identified to exist in previous work survives 
the clustering algorithm. Note that the slice shows greater non-global suppression 
than the patch at high t, in contrast to the unclustered case [61]. 



dipoles, which themselves may go on to emit. The cascade carries on until an emitted 
gluon enters the gap region. This method also allows the kt clustering algorithm to be 
directly implemented, in exactly the same way it is used experimentally, and before 
a given gluon configuration can be deemed to have polluted the gap it must survive 
the algorithm. If the gluon which has entered the gap is pulled out by the clustering 
procedure, then the dipole cascade is allowed to continue. 

Note that there are a lot of configurations with many gluons at high t and this, 
coupled with the fact that the speed of the kt algorithm scales like N"^ (where N is 
the number of gluons), means that the Monte Carlo algorithm can be very slow in 
this region. The result is that we have large statistical errors for t > 0.3, so we do not 
show this region. 

We have verified the results obtained by Dasgupta and Salam [61], which have 
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Figure 5.8: The phenomenological implications of clustering for the case of Q being a 
slice of rapidity of width 1.0. The reduction of the non-global suppression factor when 
clustering is included can be seen. 



been calculated with no clustering requirement imposed on the gluons, and written a 
program, based on their program, to perform the all-orders calculation with cluster- 
ing. Figure 5.7 shows the function S{t) for two different geometries for Q: a slice in 
rapidity of width Ar} = 1.0 and a square patch in rapidity and azimuthal angle of side 
length Ar] — A0 = 2.0, with the requirement of clustering on the final state gluons. 
The curves labelled 'exponentiation of the fixed order result' are those obtained if 
there were a simple exponentiation of the <S'2 term. 

Firstly the figure shows that at high t {t > 0.2) the suppression increases faster 
for the full calculation than for the exponentiation of ^'2 with clustering and the 
two curves have different shapes in this region. This implies a more complex geometry 
dependence in the full result than is seen by the simple exponentiation of 82- Secondly 
we see that the t dependence of the suppression is geometry (i.e. the definition of fl) 
independent at high t; the implication of this is that the clustering process maintains a 
so-called buffer region^ of suppressed intermediate radiation around Q. Such a buffer 
^The buffer region is defined by the absence of any reconstructed jets within it. It 
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mechanism has been postulated to exist in the unclustered case, and be responsible 
for the lack of gluons in the gap region. Therefore figure 5.7 shows that the all-orders 
result with clustering is more complicated than simple exponentiation of the fixed 
order result. Our work also indicates, due to the geometry independence at high t, 
that a buffer region exists around the gap which is responsible for the suppression 
of radiation into the gap. 

Figure 5.8 shows the reduction of the phenomenological significance of the NGLs 
when we cluster the final state. The figure shows the function E(t) with only primary 
emissions and the full all-orders treatment with and without clustering effects. This 
is done with Q defined as a slice in rapidity of size Arj — 1.0. There are several points 
to note. Firstly the effect of the NGLs is a large suppression of the cross section 
relative to the primary-only result. The value of t which we can consider the realistic 
upper limit for Tevatron Run II experiments is about 0.15 so we shall take t = 0.15 as 
our reference value. We can translate this into an energy scale by using the running 
couphng definition of t, equation (5.5), and obtain Q ~ 100 GeV for Qn — 1 GeV. 
In this region the inclusion of the non-global effects without clustering increases the 
suppression relative to the primary-only result by about a factor of 1.65. When we 
include clustering effects, this difference is reduced to about 1.2. Therefore, at all 
orders, the requirement of clustering on the final state reduces the phenomenological 
significance of the non-global effects by about 70%. This reduction in magnitude of 
the effect can be seen to persist to all orders. Hence when calculating cross sections, 
if we exclude the effect of NGLs then we will overestimate the cross section by 65% 
for a non-clustered final state and by 20% for a clustered final state. At larger t 
values, the overestimation increases. For comparison, the typical errors on the HI 
gaps-between-jets data is ~ 30%. 

Finally, figure 5.9 shows the same comparison of the full result, with and without 
clustering, for a patch of size Arj = A0 = 2.0. We can see that the conclusions we 

made for figure 5.8 apply and the effect of the NGLs is of similar magnitude. 

may therefore contain gluons, provided that they get pulled out of the buffer region by the 
clustering algorithm. 
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Figure 5.9: The effect of the clustering of the gluons on the function S{t) for a patch 
of phase space of size Ar] = A0 = 2.0. 



We have also performed calculations to see how the non-global suppression is af- 
fected by varying the radius parameter, R, of the clustering procedure. Figure 5.10 
shows the full function with NGLs included and with varying R. By decreasing R 
the impact of the clustering algorithm is reduced and the effect of the NGLs is re- 
stored to the non-clustered case. In fact, we expect the magnitude of the non- global 
suppression to tend to the non-clustered case as i? ^ 0. Similarly, increasing the 
radius causes more gluons to be included in the clustering and hence the magnitude 
of the suppression to reduce. In fact, we have found that for R close to 1.5 the full 
result is almost identical to the primary result. 

In summary, non-global logarithms are important not only from the point of view 
of correctness of the leading logarithm series but also result in significant numerical 
corrections to cross sections. These corrections are reduced by about a factor of 3 if we 
cluster the final state. It is clear, therefore, that while the use of the kt algorithm has 
aided the control of the NGLs, they still have a significant numerical effect at HERA. 
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Figure 5.10: The effect of varying radius parameter on the full function Eq, for a slice 
in rapidity of width A?7 =1.0. 



5.5 Ongoing research topics in non-global observ- 



ables 



Non-global observables are a new and exciting research area in QCD. In this review 
section, we will give an brief overview of some current research directions and the 
progress made. 



5.5.1 The 4 jet system 

In this section, we will introduce the calculation of NGLs for a 4 jet system, working 
in the large Nc limit at fixed order. This is physically motivated by the wish to 
understand the effect of NGLs for hadrons to jets in, for example, rapidity gap event 
in photoproduction. The main difference between a 2 jet system and a 4 jet system is 
the greater number of dipolc radiation sources in the latter case. We model the 4 jet 
system in the large Nc limit, using non-interacting and recoilless dipoles. The initial 
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state jets have 4-momenta, 



a'' = 0,0,1), 

= |(1, 0,0,-1), (5.71) 

and the final state jets have 4-momenta, 

= p,^(cosh(A77/2),l,0,sinh(A77/2)), 

= Pi^(cosh(A77/2),-l,0,-sinh(Ar;/2)). (5.72) 

In this work we assume that the jets are produced at fixed Ar) — \r)i—r)2\- The matrix 
element squared for emitting a gluon can be written, in the large Nc limit as 

where the first term describes the probability of the gluon being emitted from the {ad) 
dipole, and the second term describes the probability of the gluon being emitted from 
the (be) dipole. There is, of course, other terms describing s and t channel dipole 
emission but we initially restrict ourselves to equation (5.73) for simplicity. In these 
equations the notation (ij) means the scalar product of the corresponding 4-momenta 
and we write the 4-momentum of the emitted gluon by 

1^ = Pt (coshr^i, COS01, sin 01, sinhr^i). (5-74) 

The matrix element squared for the first gluon to emit a second gluon is, in the same 
approximation 

^ ... jod) (Ca (al) Ca {Id) \ 

^ \al)ild) V 2 (a2)(21) 2 id2){21) J 

+ ^^^(w)Mv^MM^P)(2T)y' ^ ^ ^ 

where 

1^ =Pt(coshr72,cos02,sin02,sinhr72)- (5.76) 

By preceding in analogy to the fixed order calculation for a 2 jet system, we may 
extract the non-global suppression factor for a hadron-hadron process which pro- 
duces 2 jets at fixed A77. The physical expectation is that the greater quantity of 
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dipole sources will give a larger non-global suppression factor for a 4 jet system, rel- 
ative to the 2 jet case. The 2-jet all orders calculation can be extended in a similar 
way, to extract the all-orders 4-jet non-global suppression factor, as a function of the 
gap width and the time-energy variable t. Such a study is being performed at present. 

5.5.2 Event shape/energy flow correlations 

Another way of reducing the phenomenological impact of NGLs, in addition to the 
use of clustering algorithms, has been proposed [75], and is the study of associated 
distributions in two variables. In this work, one combines measurement of a jet shape V 
in the whole of phase space (for example thrust, V = 1 — T) and that of the transverse 
away-from-jets energy flow £^out- The former is a global measurement and the latter 
is a non-global measurement. If the observable V selects 2-jet-hke configurations, one 
measures the associated distribution. 



where Q is the hard scale. It has been shown that this distribution factorises [76] , 



where V) is the standard global distribution of V and Sout(^Q) -^out) contains the 
logarithmic distribution in E'out- This latter distribution, containing NGLs is evaluated 
at the reduced scale VQ, and hence the logarithmic terms will be {ug log{VQ/Eout))"'- 
The work of Berger, Kiics and Sterman [75] considered the region in which VQ 
and Eout were comparable, so that the NGLs give a negligible contribution. Thus, 
for a restricted subset of appropriately selected events, it is possible, to 'tune out' 
the non-global logarithmically enhanced terms in associated distributions. Hence the 
use of event shape/energy flow correlations is another method, along with clustering 
algorithms, to control the effect of NGLs in non-global observables. 



(5.77) 



E2ng(g, V, Eout) = S(g, V) ■ E,^t{VQ, E,^,), 



(5.78) 
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5.6 Conclusions 

The observation that observables that are sensitive to radiation in only a restricted part 
of phase space, so-called non-global observables, are strongly sensitive to secondary 
radiation is a new and exciting discovery. For a long time it was widely thought, 
now it seems incorrectly, that it was sufficient to consider only primary emission 
contributions to such observables. These primary, or Bremsstrahlung, contributions 
are well understood for a variety of processes. Recent measurements of gaps-between 
jets events at HI and ZEUS are exactly the class of observable that are sensitive to 
these effects and to deal with this fact, a study of NGLs in the context of HERA 
measurements is required. In this chapter we considered the NGLs generated from 
secondary radiation into a restricted region of phase space under the condition of 
the final state surviving a clustering algorithm. Such kt clustering algorithms have 
been used for the HERA measurements. Our work verified the study of Dasgupta 
and Salam, who found that inclusion of NGLs resulted in a strong suppression of the 
theoretical prediction of the non-global observable. 

We have also studied the NGLs for a 2 jet system with clustering imposed and 
found that the clustering process reduced, but did not eliminate, the non-global loga- 
rithm effect. 

Our main conclusion is that the final state specified in the HI and ZEUS analysis 
will mean that a primary emission calculation in the manner of Sterman et al [55,58, 
59], and chapter 4, of energy flow observables will overestimate the observed gaps- 
between-jets rate by around 20%. This is to be contrasted with an overestimation 
of 65% that would be found for a non-clustered final state. This result was calculated 
to all orders in the large Nc limit for a 2 jet system. The extension of both these results 
to beyond the large Nc limit would be a major step forward in the understanding of 
these effects in non-global observables and we reserve this for future work. 



Chapter 6 

Resummation of energy flow 
observables 

6.1 Introduction 

The subject of interjet energy flow [77] has attracted considerable interest ever since 
it was proposed [78, 79] as a way to study rapidity gap processes using the tools of 
perturbation theory. Rapidity gap processes are defined as processes containing two 
high pt jets with the region of rapidity between the jets containing nothing more than 
soft radiation. This region is known generically as the gap. The presence of a range of 
scales offers a chance to study the interface between the soft, non-perturbative scales 
and the hard, perturbative scales of QCD. 

In this chapter we will calculate the perturbative contribution to gaps-between- 
jets cross sections in photoproduction at HERA, which have been measured by the 
HI [62] and ZEUS [63, 80] collaborations. A feature of the recent analyses is the 
use of a clustering algorithm to define the hadronic final state and hence the gap. 
The restriction of transverse radiation in a region of phase space, defined as Q and 
directed away from the observed jets and the beam directions, produces logarithms 
at each order of QCD perturbation theory of the interjet energy fiow, Qa, over some 
hard scale, Q. The precise definition of the restricted region, or gap, is totally free 



Chapter 6. Resummation of energy flow observables 



113 



and in this work we are interested in the gap region defined by experimental rapidity 
gap analyses. The source of the large logarithms is twofold. The so-called primary 
(or global) logarithms arise from radiation emitted directly into f2; these wide-angle 
gluons decouple from the dynamics of the colour-diagonal jets and are described by 
an effective, eikonal theory [7,53,55,58,59]. The second source of leading logarithms 
arise from gluons emitted outside of the gap region, an area of phase space generically 
denoted as Q, which subsequently radiate into Q. These terms are known as non-global 
(secondary) logarithms, or NGLs [60,61,70,71], and were introduced in chapter 5. 

The primary logarithms are resummed using the formalism of Collins, Soper and 
Sterman (CSS) [6,7,48,52], which was introduced in chapter 4. In this method the 
cross section is factorised into a soft part describing the emission of soft, wide angle 
gluons up to scale Qq and a hard part, describing harder quanta. A unique feature of 
QCD is that the soft and the hard functions are expressed as matrices in the space of 
possible colour flow of the system. The scale invariance and factorisation properties 
of the cross section are then exploited to resum primary logarithms of Qn/Q. This 
resummation is driven by the ultraviolet pole parts of eikonal Feynman graphs and 
we write the resummed cross section in terms of the eigenvalues of r2-dependent soft 
anomalous dimension matrices. These matrices are known for gap definitions based 
on the cone definition of the final state [58,59] and for a gap defined as a square patch 
in rapidity and azimuthal angle [55]; here we are interested in gaps defined in terms 
of the clustering algorithms employed in the recent analyses. Hence we are required 
to calculate the corresponding anomalous dimension matrices. 

The NGLs [60, 61] are unable to be incorporated into the resummation of the 
primary logarithms, because the gluon emission patterns that produce the NGLs are 
sensitive to underlying colour fiows not included in the formalism. The effect of 
NGLs, which is a suppressive effect, on energy flow processes has been studied using 
numerical methods in the large Nc limit and overall factors describing their effect have 
been extracted for a two jet system, both without [61] and with [70] clustering. The 
latter case was described in chapter 5. This factor is not directly apphcable to the 
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4 jet systems^ relevant in the photoproduction of jets but, in the lack of a four jet 
formalism, we nevertheless use the 2-jet factor in our predictions. 

Our aim is to derive LL resummed predictions for the gap cross section, with pri- 
mary logarithms correct to all orders and secondary logarithms correct in the large Nc 
limit. The gap cross section will follow the HERA analyses and demand two hard 
jets, defined using the kt clustering algorithm [64-66], and we will closely follow the 
HI and ZEUS gap definition. The technical aspects of soft gluon resummation give a 
strong dependence on the gluon emission phase space, and hence a considerable part 
of our work will be concerned with the calculation of soft gluon effects for the specific 
detector geometry of the HI and ZEUS experiments. 

The organisation of this chapter is as follows. Section 6.2 describes, in detail, the 
energy flow analyses of HI and ZEUS. We summarise the experimental cuts employed 
and the range of measured observables. We also discuss the theoretical implementation 
of the inclusive kt algorithm employed to define the hadronic final state and the 
impact on soft gluon resummation. Section 6.3 describes the theoretical definition of 
our cross section and we employ the standard QCD factorisation theorems to write 
it as the convolution of non-perturbative parton distributions and a short-distance 
hard scattering function. We then proceed to refactorise the hard scattering function 
and exploit this factorisation to resum the large interjet logarithms. Section 6.4 then 
derives the soft anomalous dimension matrices for the kt defined final state and in 
section 6.5 we present detailed predictions of rapidity gap processes and compare 
to the HI data. Finally we draw our conclusion in section 6.6. We find that our 
description of the data is good, although the approximate treatment of NGLs results 
in a relatively large normalisation uncertainty. The research work in this chapter is 
pubhshed in [81]. 

^Note that for a two-to-two process the incoming and outgoing partons radiate, so we 
consider the process to be of "four jet" type, although only two jets are seen in the final 
state. 



Chapter 6. Resummation of energy flow observables 



115 




7] 

Figure 6.1: The experimental signature of a "jet-gap-jet" event. The dark circles 
denote two jets, and the gap region is bounded by the dashed lines. 

6.2 The HERA energy flow analyses 

In this section we will outline the experimental analyses of the photoproduction of 
gaps-between-jets processes and summarise the experimental cuts and rapidity gap 
observables. We will also describe the clustering algorithm used to define the final 
hadronic state in the more recent HI [62] and ZEUS [63] analyses. 

The data for these events were collected when HERA collided 27.6 GeV positrons^ 
with 820 GeV protons, giving a centre of mass energy of y/s ~ 300 GeV. Following 
the jet-finding phase, which we will comment on later, the total transverse energy 
flow between the two highest Et jets, denoted E^^^ , is calculated by summing the 
transverse energy of all jets produced by the kt algorithm that are not part of the 
dijets in the pseudorapidity region between the two highest jets. An event is then 
defined as a gap event if the energy is less than some energy cut E^^"^ = Qq. A gap 
fraction is then calculated by dividing the cross section at fixed E^^"^ by the inclusive 
cross section. The ZEUS collaboration performed a rapidity gap analysis several years 
ago [80] using the cone algorithm for the jet definition and presented the gap fraction 
at Qn = 0.3 GeV. We consider this value of Qq as being too small for our perturbative 
analysis and will not make any predictions for this data set. The gaps-between-jets 
experimental signature is illustrated in figure 6.1. 

^The positron energy varied a negligible amount between the two sets of analyses. 
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HI 


ZEUS 




> 6.0 GcV 


> 6.0 GeV 




> 5.0 GeV 


> 5.0 GcV 


^jctl 


< 2.65 


< 2.4 


^jct2 


< 2.65 


< 2.4 


Ar] 


2.5 < Ar^ < 4 


2 < Ar^ < 4 




N/A 


< 0.75 


y 


0.3 <y < 0.6 


0.2 <y< 0.85 




< 0.01 GeV2 


< 1 GeV2 


jet def. 


kt 


kt 


gap dcf. 


Ay = Arj 


Ay = Arj 


/? 


1.0 


1.0 


Qn 


0.5,1.0,1.5,2.0 GeV 


0.5,1.0,1.5,2.0 GeV 



Table 6.1: The experimental cuts used for the HERA analyses. 

The more recent HI and ZEUS analyses used the kt definition of the final state 
and both collaborations presented the gap fraction at four different values of Qn, as 
shown in table 6.1. We will make predictions and compare to data for the HI data 
sets and, due to the fact that the ZEUS data is still preliminary, confine ourselves 
to making predictions for the ZEUS analysis. We have summarised the cuts used in 
table 6.1. 

6.2.1 The kt algorithm 

Of special interest to those going about soft gluon calculations is the method used to 
define the hadronic final state, the reason being that this jet- finding process determines 
the phase space for soft gluon emission; the method used in the HI and ZEUS data 
sets is the kt algorithm [64-66]. Here we follow the inclusive scheme used at HI and 
ZEUS which depends on the parameter R, normally set to unity. If we assume that 
any radiation into the gap is much softer than any parent radiation, then this radiation 
with Et < -E^* will be merged into the jet (with kinematical variables (r/jet, 0jet)) if it 
satisfies 

(?7r-?7jet)' + (0r--0jet)' <i?', (6-1) 
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where we denote the kinematical variables of the radiated gluon by {r]r,(t)r)- Once 
merged, a gluon will be pulled out of the gap and can no longer produce a primary 
or secondary logarithm. The gap is defined as the interjet region minus the region of 
clustered radiation around the jets and may contain soft protojets. The gap transverse 
energy is then defined by the (scalar) sum of the protojets within the gap region, 
rji < T] < 'r]2- We have given a detailed description of how the kt algorithm works in 
chapter 5. 

The kt gap definition can be contrasted to the older ZEUS analysis [80], which 
used the well known cone definition of the final state with R — 1.0. The gap transverse 
energy is then defined as the scalar sum of the transverse energy of the hadrons within 
it, rji + R < 7] < r]2 — R. 

6.3 Factorisation, refactorisation and resummation 
of the cross section 

In this section we will exploit the standard factorisation theorems of QCD to write 
down the dijet production cross section from the interaction of a proton and a positron. 
We will then rcfactorisc the hard scattering function into the product of two matrices 
in the space of possible hard scattering colour flows, one matrix describing soft gluons 
radiated into the gap region and the other a hard scattering matrix. The renormalisa- 
tion properties of the cross section are then used to resum primary interjet logarithms, 
and write the result in terms of the eigenvalues of the matrix of counterterms used to 
renormalisc the soft function. In the following section we will calculate these matrices 
and their eigenvalues. 

6.3.1 Photoproduction cross sections 

The scattering of positrons and protons at HERA proceeds predominantly through the 
exchange of photons with very small virtuality and produces a large subset of events 
with jets of high transverse momentum, Et- The presence of this large scale allows 
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the apphcation of the perturbative methods of QCD to predict the cross section for 
multiple jet production. This process is otherwise known as jet photoproduction. 

The leading order (LO) QCD contribution can be divided into two types [59]. The 
first is the direct process in which the photon interacts directly with a parton from the 
proton and proceeds through either the Compton process, 7g — > gq, or the photon- 
gluon fusion process, — > qq. The second contribution is the resolved contribution, 
in which the virtual photon fluctuates into a hadronic state that acts as a source 
of partons, which then scatter off the partonic content of the proton. Therefore the 
reaction proceeds through standard QCD 2 — > 2 parton scattering processes. Note that 
the precise determination of the partonic content of the photon is a very open question 
and there is a relatively large error associated with the photonic parton densities. The 
spectrum of virtual photons is approximated by the Weizacker- Williams [82] formula. 



where me is the electron mass, y is the fraction of the positron's energy that is trans- 
fered to the photon, and Q^ax maximum virtuality of the photon, which is 

determined by the experimental cuts employed in the analyses. Then, by using the 
equivalent photon approximation, the cross section for the process e+p — > e^X is given 
by the convolution 



where we write d(7(7p — > X) for the cross section of 7p — > X. The centre of mass 
energy squared for the photon- proton system is W'^ — ys, where s is the centre of 
mass energy squared for the positron-proton system. At HERA, s ~ 90, 000 GeV^ 
and the values for y^im and |/max are determined by the experimental analyses. We 
can now write down the specific expression for the production of two high Et jets 
from the photon-proton system, which is written as a sum of the direct and resolved 
contributions. 




(6.2) 




(6.3) 



dae+p{s, i, Arj, as{Hr), Qo) ^ dy F^i^{y) i da^^ (s^p, £, Ar^, as{Hr), Qn) 




(6.4) 
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where we denote the 4-momentum transfer squared in the hard scattering as t and 

have temporarily suppressed the function arguments on the right-hand side. We define 
the rapidity difference and average of the two hard jets by 

Af] = |r/i-7/2|, 

Vjj = ^{Vi + V2)- (6.5) 

At this point we can appeal to the factorisation theorems of QCD and, by working in 
the 7p frame, write down factorised forms for the direct and resolved cross sections. 
The factorised direct cross section is 

—r^{s^p,i,Arj,as{f^r),Qn)^ V / dxpCpf^/pixp, Hf) 

jp,hJ2 

X {s,t,Ari,as{f^r),Qn,f^f), (6.6) 

and the factorised resolved cross section is 



^{s^p,i,Ar],asii^r),Qa)= ^ / dx^ dxp(f)f^/^{x^, Hf)(pf^/p{xp, Hf) 

,/p,/i,/2 -^^^ 



x-^{s,t,Ar],as{i^r),Qa,l^f) (6.7) 

which are written in terms of the jet rapidity, f), in the partonic centre-of-mass frame, 
and we write the factorisation scale and the renormalisation scale as /// and /ir re- 
spectively. Note that fj — Arj/2, s — XpW"^ for the direct case and s — XjXpW'^ for 
the resolved case. In these equations we denote the integration regions of the direct 
and resolved convolutions, which are defined by the experimental cuts, by Rd and Rf. 
The parton distribution for a parton of fiavour / in the photon and the proton are de- 
noted by (j)fj^{x^,iJLf) and (l)f/p{xp,iJ,f) respectively and finally ^^^^ and are the 
hard scattering functions which, at lowest order, start from the Born cross sections. 
These are the functions that will contain the logarithmic enhancements of Qn/Q, and 
hence depend on the definition of the gap Q and the gap energy fiow Qq. We assume 
that Qq is sufficiently soft that we can ignore the effects of emission on the parent 
jet, known as recoil, but large enough so that ^ ^qcd- The index / denotes the 
process /y + /p — > /i + /2 and the index /7 denotes the process 7 + /p — > /i + /2- 
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Since tfie aim of this chapter is to calculate ratios of cross sections and compare 

with data, we will take the rcnormalisation scale to equal the factorisation scale and 
set /If = = Pt, where pt is the transverse momentum of the produced jets. 

6.3.2 Refactorisation 

Following [7,55], and the arguments in chapter 4, we now refactorise the 2^2 hard 
scattering function into a hard matrix and a soft matrix. 



We introduce a factorisation scale separate to the parton distribution factorisation 
scale fif, and all dynamics at scales less than /j, are factored into Su- Therefore Hji^ 
is Qn independent, and all the Qq dependence is included in S^j. This latter function 
describes the soft gluon dynamics. The proof of this statement follows standard fac- 
torisation arguments [7] . The indices / and L label the basis of colour tensors which 
describe the possible colour exchange in the hard scattering, over which the hard and 
soft matrices are expressed. Soft, wide angle radiation decouples from the dynamics of 
the hard scattering and can be approximated by an effective cross section and in this 
effective theory the partons are treated as recoilless sources of gluonic radiation and 
replaced by eikonal hues, or path ordered exponentials of the gluon field [53]. The soft 
radiation pattern of this effective eikonal theory then mimics the radiation pattern 
of the partons participating in the hard event, or in other words the effective eikonal 
theory will contain the same logarithms of the soft scale as the full theory. The hard 
scattering function will begin at order for the resolved process and order aas for 
the direct process, and the soft function will begin at zeroth order. The lowest order 
soft function, denoted 5'£j, reduces to a set of colour traces. Note that the definition 
of the gap, and hence the soft function, depends on the jet separation A77 but we have 
suppressed this argument of the soft function for clarity. 

The construction of the soft function, and in particular its rcnormalisation prop- 
erties, have been extensively studied elsewhere [7,53], and are reviewed in chapter 4. 



dfj 




Chapter 6. Resummation of energy flow observables 



121 



A non-local operator is constructed from a product of Wilson lines, which ties four 
lines (representing the four jet process) together with a colour tensor. This operator, 
which contains ultraviolet divergences and hence requires renormalisation, is used to 
construct a so-called eikonal cross section, which serves as an effective theory for the 
soft emission. By summing over intermediate states the eikonal cross section is free 
of potential coUinear singularities. It is the ultraviolet renormalisation of the eikonal 
operator that allows colour mixing and the resummation of soft interjet logarithms. 

6.3.3 Factorisation leads to resummation of soft logarithms 

The partonic cross section, which has been factorised into a hard and a soft function, 
should not depend on the choice of the factorisation scale /i, 



It is important to point out that we have deliberately ignored the complications of 
terms in this equation arising from radiation into Cl, as in [55] and chapter 4, and 
only include radiation emitted by the soft function directly into fl. The implication 
of ignoring these non-global terms is discussed in section 6.3.4, where we also describe 
how to include their effect in a different way. Therefore we have never included the, 
technically correct, Q argument of the soft function. The matrices Tsifj, il) are process- 
dependent soft anomalous dimension matrices which depend on the details of the 
gap definition and the hard scattering. This equation is solved by transforming to a 
basis in which these matrices are diagonal and hence we require a knowledge of the 
eigenvectors and eigenvalues of the soft anomalous dimension matrices. We obtain 
the entries for Ts{r], f2) from the coefficients of the ultraviolet poles in the matrix of 
counterterms which renormalise the soft function; we can write the renormalisation 
constant as a sum over terms from different eikonal lines, each with the form of a 




(6.9) 



Following chapter 4, this means the soft function obeys 




(6.10) 



r _ 



r 
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colour factor multiplied by a scaleless integral: 

{Zs)li = = E^J;?^^^^'^- (6.11) 

The eikonal momentum integrals are process independent, and only depend on i and j, 
the eikonal lines that are connected by the virtual gluon. The colour factor is found 
from consideration of the colour flow for a given process and the basis over which 
the colour flow is to be decomposed. The result is a basis- and process-dependent 
set of colour mixing matrices, which we have listed in appendix F, together with our 
choice of bases in appendix D. The colour mixing matrices were discussed in detail in 
chapter 2 and have been obtained in [8,9,55] for all relevant subprocesses, and involves 
using SU(3) colour identities like 

mi = l^uhj - (6-12) 

for quark processes and 

dabc = 2[Tr(m^)+Tr(rtV)], (6.13) 
fabc = -2z[Tr(rt¥)-Tr(n^t^)], (6.14) 

for gluon processes, to decompose one-loop graphs over a colour basis. A more com- 
plete set of SU(3) group theoretic identities can be found in appendix A. We use the 
fact that the colour flow for a real graph is the same as the corresponding virtual 
graph, valid for primary emission. 

Therefore we need to calculate the ultraviolet divergent contribution to the mo- 
mentum function cj*^*-'^ from all contributing eikonal graphs. Working in the Feynman 
gauge there are two possible sources of divergence. The first is one loop eikonal graphs 
with a virtual gluon connecting eikonal lines i and j. From the eikonal Feynman 
rules listed in the appendix, these graphs will give a real and an imaginary contri- 
bution to r^. Note that as we are working in the Feynman gauge the self energy 
diagrams (cu*^'"-') give no contribution. The second source of ultraviolet divergences are 
the real emission diagrams, when the emitted gluon is directed out of the gap. This 
can produce an ultraviolet divergence in the eikonal graph as we only measure energy 
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Figure 6.2: The virtual eikonal graph obtained if the virtual gluon connects eikonal 
lines a and h. In the notation in the text, this diagram is denoted Iv'^^\ 



flow into the gap and are fully inclusive out of the gap. Hence the virtual graphs will 
only depend on the relative direction of the two eikonal lines and the real graphs will 
give a gap (and hence a jet algorithm) dependence. This sum over real and virtual 
eikonal graphs ensures that the soft function remains free of collinear divergences. The 
imaginary (and geometry independent) part of all our anomalous dimension matrices 
can be extracted from [8,53,55,58], and the calculation for a cone-algorithm defined 
final state has been done in [58]. For the latter case, we have re-expressed their results 
in accordance with our notation in appendix H. 

By performing the energy integral of the virtual graphs, we can combine the 
result with the corresponding real graph at the integrand level and obtain a partial 
cancellation. Therefore we start by evaluating a generic virtual eikonal graph, with 
the virtual gluon connecting eikonal lines i and j, using the eikonal Feynman rules in 
appendix B. We obtain 

= gl^AA ■ j (2vr)'^ (A;2 + ze) [5 A ' k + le) {5,(3^ ■ k + te)' ^^'^^^ 

where 6ij = ±1 and we denote the momenta of the virtual gluon by k. In this chapter 
we use the dimensional regularisation convention d = 4 — 2e. Figure 6.2 shows the 
virtual eikonal diagram li"^''. Note that in this section we will use the symbol e to 
denote a UV pole and also as a small parameter when doing contour integration. 



r 1 1 



Chapter 6. Resummation of energy flow observables 



124 



However there should be no confusion. Writing 



d'^k = (6.16) 

k^" = (6.17) 

/3f = (1,A), (6.18) 

e + ie = -\k\+ie){k° + \k\-ie), (6.19) 



we get 



7(^^) = q'^A A -3-3- I dk^- - ^- - - 

" '^'^'J (27r)<^((A;0)2_^2 + ,gj(^.[^o_^..^]+,e)(5,[A;0-/^-^]+^e) 

(6.20) 

We need to evaluate the integral using contour integration, and there are four com- 
binations of the d[s to consider; however the integrals only depend on the product diSj 
and so it is sufficient to only consider two of the four possible combinations of 5i and 5j. 
The four poles in the complex k^ plane are 

k^ = +\k\-ie, 

A;° = {SA-k-ie)/Si, 

A;0 = (5jPj-k-ie)/5j. (6.21) 
We consider the case oi 6i — Sj — —1, for which the poles are 

k^ = -\k\+ie, 

fcO = +1^1 - ie, 

k^ = j3i-k + ie, 

k^ = p..k + ie. (6.22) 

We close the contour in the lower-half plane and only pick up the pole k^ — +\k\ — ie. 
The residue of this pole is extracted from the standard result of the Laurent expansion 

••• + ^^|.=.o + ... (6.23) 

Z — Zn 
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and we get for the residue (where we have used equation (6.19) to obtain the desired 
form) 

— 2 

a_i = =; =; =; =; z; =; =; . (6-24) 

(2|^| -ie)(-\k\ + (3i-k + ie)(-\k\ + (3j -k + ie) 
Therefore, multiplying by 27ri and by —1 (for going clockwise around the contour), we 

get 

7^-^^ , (6-25) 

{27:Y-^2\k\{-P,-k){-Prk) 

as e — > 0. We note that this expression is the same as the real diagram expression 
integrated over all of phase space, and with an overall minus sign. This real diagram 
expression is readily obtained from the evaluation of an eikonal diagram with real 
gluon emission using the eikonal Feynman rules in appendix B. We can also see that, 
if we have no restrictions on the real gluon emission energy, we expect a complete can- 
cellation of real and virtual diagrams. Now consider the case of (^j = +1 and 5j — — 1. 
This configuration corresponds to virtual diagrams where the virtual gluon connects 
the incoming eikonal lines. Now we have the following arrangement of poles, 

fcO = +\k\-ie, 

k —— f3i ' k 

= p. -k + ie. (6.26) 

with two in the upper half plane, and two in the lower-half plane. We choose to 
close the contour in the lower-half plane again and pick up two poles. The residues 
are computed in exactly the same manner as the first case and we obtain, for the 
k^ = +\k\ — ie pole 

-i 1 1 , , 

a_i = ^ ^ , 6.27 

{2\k\-ie){\k\- I3i-k){-\k\+l3j-k + ie) 

and for the k^ — (5i ■ k — ie pole 

—i 

a_i = z; =i =i =; =; =; =; . (6.28) 

{{I3i-k-ief -k'^ + ie){l3j-k- I3i-k + ie) 
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Note that to get the correct form of the first residue, we need to use the form of the 
integrand containing equation (6.19), and to get the correct form for the second residue 
we need to use the k"^ = (/c°)^ — form of the integrand. Therefore the integral for 
this arrangement of poles has two contributions - one from the k^ pole and another 
from the P ■ k pole. The first gives 

and the latter gives 

" ^ ^ ^ ^ ^ ■ (6.30) 

i27ry-'((P,.ky-k^±ie){Pj-k-Pi-k + ie) ^ ' 

We will give an explanation for the s subscript later in this section. The ± piece comes 
from the unknown sign of Pi ■ k. Note that, once again, the piece arising from the k"^ 
poles gives a contribution looking like the negative of the real contribution. At first 
sight, it looks hke the real diagrams should give a, Pi • k pole if we write the Lorentz 
invariant phase space in 4 dimensions using 

where we express the on-shell condition of real gluons as 

r(p^) = 5(p2)e(pO). (6.32) 

The Pi • k pole is not picked up by the real diagram energy integral, due to this 
condition. We have determined that if we perform the energy integral of a virtual 
eikonal graph, we obtain the same form as the corresponding (same i and j) real 
graph. We shall not consider the analysis of the remaining pole configurations, for 
which we will draw the same conclusions, and describe how everything fits together. 
We are interested in the ultraviolet pole part of the real and virtual one-loop eikonal 
diagrams, in order to calculate the soft anomalous dimension matrix. The virtual loop 
diagrams contribute through their loop integrals and the real diagrams contribute in 
all regions where we are completely inclusive of gluonic radiation (everywhere except 
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any "gap" region.) Therefore 



UJ 



(ij) 



{2TiY (A;2 + ie) {6i(3i -k + ie) {6jpj ■ k + ie) 

J (27r)'' ^ 2\k\{6if3i ■ k){6j(3j ■ k) 

In writing this hne we have used the exphcit forms of the virtual and real diagrams, 

with the virtual gluon joining eikonal lines i and j. The function G(A;) = 1 when the 

real gluon is directed into the gap, and is zero otherwise. We can now do the virtual 

integral. 



d'^-^k 1 



{27:)<^-'2\k\{-(3,-k){-(5yky 
+ glAA,P..p,j£^^il-eik))-. ' 



2\k\i5A-k)i5jPj-k)' 
d'^-^k 1 



2\k\{SA-k){Sj(3j-k) 



+ hs(l-5i5j), (6.34) 

and find that its real part can be written as — 1 x the real contribution, but with no 
phase space restrictions. This part is the bit left over from picking up the k^ pole. 
The function Ig is the result of picking up the /3j • k pole, and is only present when 
5i — —1 X 5j, or when the virtual gluon connects two initial or final state eikonal lines 
(s-channel diagrams). This is the reason for the s subscript for the function Is- We 
have also performed a real/virtual cancellation and have been left with an integral over 
the vetoed gap region, plus an additional function Is, equation (6.30). To evaluate 
this latter piece we use polar coordinates in {d — 1) dimensions, 

d'^-^k = k'^-^dk d cos e d(P d'^-^fl, (6.35) 

— * — * 

where we have noted that this integral only appears in diagrams in which Pi and Pj 
are back-to-back (and so Pi ■ Pj — 2), and write 

Pi-k = kcosO, (6.36) 
0j-k ^ -kcose. (6.37) 
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Note tfiat the a general polar d- vector for k is 

k'^ — k{l; . . . , sin 9 sin (f), sin 6 cos (f), cos 9). (6.38) 

Doing this we obtain the integrand 

1 1 



(cos2(e) - l±ie)(-2cos(^) +ie) (2cos(^) - ie)(l ^ ie + l/4e2) 

(2cos(9)+ie) 
(l-cos2(e)^ie)(4^4ie + e2)' ^ ' 

where we have used partial fractions to write the second form. We can safely set e = 
in the second term, since the only poles present are at cos^ = ±1 and these are 
protected by dimensional regularisation. Therefore the second term is an odd function 
of cos^ and integrates (over the range —1 to +1) to zero. Therefore we are left with 



/_ 



_i 2cos(^)-ie 2 
Hence we find that the imaginary piece is 



dcos(9) TT , / . X 

^ ^- = +-^-0(e), (6.40) 



TT 2e' 

where we have used the fact that 



= -^i^L2AiAi-h (6.41) 
(27r)2 2e ' ^2 ' ^ ^ 

= AAAS,^'^, (6.42) 



(6.43) 



^PP k^ 26 

for d = 4 — 2e and an integration from some fixed energy to infinity. Therefore the 
momentum integrals take the form, 



CO 



(ij) 



d'^-^k 1 



2\k\{2ny-^ ' {Sipi-k){6jpj-k)' 
+ 5,5,AA-^|(1-M,), (6.44) 

which we can write in terms of the rapidity and transverse energy of the emitted gluon, 
^ -|A<A,Mpip / ^.d;c.d,ge(g) ^^^ /^--ft ^^ +/.P.. (6.46) 
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We fiave written the geometry independent imaginary part as LP.. Once we have 
obtained the momentum integrals for the kt defined final state we can construct the 
anomalous dimension matrices using the colour mixing matrices in appendix F. Con- 
sideration of the eigenvalues and eigenvectors of these matrices, together with the 
process-dependent hard and soft matrices (the full set of hard and soft matrices is 
shown in appendix E ) allows the resummed cross section to be written down, 

hT'sT exp [UKifl, ^) + MV, ^)) r -/^o«.(/x)| , (6.46) 

which follows from the diagonalisation of the soft RGE, equation (6.10). This was 
discussed in chapter 4, where we developed this resummed expression from the fac- 
torisation properties of the cross section. We denote matrices in the diagonal basis by 
barred matrices, the eigenvalues of the anomalous dimension matrices by Aj — ctg Aj and 
we write the lowest-order piece of the QCD beta function as /3o = (H-ZVc — 2n/)/(67r). 
Note that our normahsation of the eigenvalues differs by a factor of tt from that used 
in chapter 4. We will observe that, in agreement with [58], Re(A) > for all physical 
channels and hence the resummed cross sections are suppressed relative to the fully 
inclusive cross section. 



6.3.4 Non-global effects 

As we have discussed in the last section, we have deliberately ignored terms arising 
from secondary radiation into Jl, or non-global logarithms (NGLs) [60,61,70,71]. Such 
terms arise from radiation at some intermediate scale, M, being emitted outside of Q, 
i.e. into Vl, and then subsequently radiating into VL. In energy flow observables such 
effects give rise to leading logarithms. Inclusion of NGLs in the formalism of the last 
section would result in an explicit M dependence of the soft function and a sensitivity 
to more complicated, 2 — > n, colour flows for all n > 2. For further details see [55]. 
NG effects have been studied for a two-jet system by Dasgupta and Salam [60,61], by 
Appleby and Seymour with the complication of clustering [70], and in the context of 
energy flow/event shape correlations by Dokshitzer and Marchesini [76] and Berger, 
Kucs and Sterman [75] . 
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Qu [GeV] 




0.5 


'-'•-'-u_o.io 


1.0 


0.47l^i^ 


1.5 


0.65+-!^ 


2.0 


0.741H^ 



Table 6.2: The non-global emission suppression factors for the 4-jet system, obtained 
from an all-orders calculation for Q = 6 GeV. 



The effect of NGLs for four-jet kinematics has not been explicitly calculated to 
date. The best that has been managed is a two-jet calculation in the large- A^c hmit. 
The NG contributions to the gap cross section factorize into an overall suppression 
factor S^^ , making it smaller than would be predicted by the resummation of primary 
logarithms alone. In the absence of a complete calculation, we include the NGLs 
approximately, by using our all-order results in the large- A^c limit for S^^ in a two-jet 
system [70]. Since four-jet configurations are dominated, in the large- A^c limit and for 
large A?7, by colour flows in which two colour dipolcs stretch across the gap region, 
we approximate the four-jet NG suppression factor by the square of the two-jet one. 

We have reperformed our previous calculation for the kinematic range relevant to 
HERA and find that the variation of S^^ with A77 is very weak, so we neglect it. The 
variation with Qn is very strong on the other hand. S^^ is a function of 

where (5q — (llCyi — 2n/)/(67r), and is well-approximated by a Gaussian in t. Thus 
if Qq, is too close to Aqcd, t varies rapidly with it and S^'^ varies very rapidly. 

It is impossible to quantify the uncertainties in this approximation, without a 
more detailed understanding of the underlying physics. To get an idea however, we 
estimate the possible size of higher order corrections, by varying the hard scale at 
which as is evaluated. To leading logarithmic accuracy, this is equivalent to varying 
the value of as{Q) by an amount proportional to its value. Therefore we estimate the 
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uncertainty by 



as{l + as), 



(6.48) 



down 



as{l - as). 



(6.49) 



We tfierefore evaluate t, and hence 5'^*^(t)^, using our central value of as{Mz) = 
0.116, which implies as{Q = 6 GeV) = 0.196, and with raised and lowered values 
a^P(6 GeV) = 0.234 and af™°(6 GeV) = 0.158. For Qa = 1.0 GeV, for example, 
these values result in t = 0.097l^;g| and hence S^^{t)^ = 0.471^^^. We show the 
results for all relevant values of Qq in table 6.2. Note that Qn — 0.5 GeV is so low 
that the range of uncertainty in t extends beyond Aqcd and hence the estimate of S^'^ 
extends to zero. We have not shown any results for the 1995 cone-based ZEUS energy 
flow analysis [80] because the low value of Qq — 0.3 GeV means that the central value 
of the NG suppression is already zero, indicating a breakdown of our perturbative 
approach. 

The uncertainty on the secondary emission probability estimated in this way 
should be added to that on the primary emission probability, described in section 6.5. 
However, we will see that the secondary uncertainty generally dominates the two. 
This is therefore clearly an area that needs more work if more precise quantitative 
predictions are to be made. 

6.4 Soft gluon dynamics for a kt defined final state 

We now evaluate the momentum integral, uj^'^^\ over the gap region Q. The region 
of integration is determined by the experimental geometry, in which the final state is 
defined by the kt algorithm, and we shall work with the quantity 




(6.50) 



where we define k — k/kf. Therefore 



CO 



2^ 



(6.51) 
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f 



n 

Figure 6.3: The phase space regions for a kt defined final state. The shading denotes 
the regions vetoed by the algorithm, which are subtracted from the f2/ piece. Note 
that we have dropped the {ij) superscripts in this figure. 

We denote the geometry independent imaginary part by /.P., and we define the finite 
piece r*^*-'^ by 



UJ 



-s.. 



2e 



(6.52) 



We have extracted the sign function from V^''^\ 



Sij = AiAj6i6 



(6.53) 



so that 



TT 



LP.. 



(6.54) 



In this work we denote the rapidity separation of the jets by Arj and the width of 
an azimuthally symmetric rapidity gap by Ay (< Arj). Therefore the available phase 
space for soft gluon emission for a kt defined final state is given by 



- Jim^ {nf\Ay,Ar^)-n'l'>iAy,Ar^,R)-n'^'>iAy,Ar^,R) 



(6.55) 



where the first term arises from an azimuthally symmetric gap of width Ay, and we 
subtract the region around each jet which is vetoed by the kt algorithm. The regions 
of this equation are shown in figure 6.3. In these regions any soft radiation is clustered 
into the jet, and cannot form part of Q. In the first term we take Ay approaching Arj, 
and hence it contains a collinear divergence when the emitted gluon is coUinear to one 



Chapter 6. Resummation of energy flow observables 



133 



of the jets. The two subtracted pieces then remove the regions of phase space defined 

by 

{'nk-Vif + {4>k-(t>if<R\ (6.56) 

where the index i labels final state jets and k labels the emitted gluon. The coUinear 
divergences in the subtracted pieces exactly match the divergences in the first piece and 
hence the function Vf^^^^ ( is coUinear safe. Explicit definitions of the Vt functions 
are 

^ J-Ay/2 ^Jo 2nif3,-k)iPj-ky 

^ = / d0 /^^ (6.57) 

JAr,/2-R J-4>n^ {(5i ■ k){(5j ■ k) 

where we write = ^JE? — [r] — Arj/2)'^ and obtain 0,2''^ by the symmetry Og'''* = 
where the mapping i — > z is given by {a, b, 1, 2} — > {b, a, 2, 1}. If we define the 
following combinations of momentum integrals, 

7 = 56ir(''^) + 5„2r("2), (6.58) 

where we have combined classes of diagram with the same colour structure, we obtain 
the following closed form for the positive gap contributions, in the limit Ay A77, 

>i7r), (6.59) 





as 


a — 






TT 


P = 


as 






TT 




as 


7 = 





1 


Ar] 


-Ay 




1 


Ar) 


-Ay 






Ar] 


-Ay. 



(6.61) 



The subtraction pieces are straightforward to express as power series in R and e and 
we shall illustrate the calculation of the momentum integrals with an example. 
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6.4.1 Calculation of fl^f ^ 

We can write the matrix element in terms of the rapidity of the emitted gluon and 
obtain the following matrix element 



{I3a ■ k){l3i ■ k) e-''(cosh(Ar?/2 - r]) - cos(/))' 
using the 4-momenta products, 



(6.62) 



= cosh(Ar7/2) - sinh(Ar?/2), (6.63) 

13a -k = cosh(r7/2) -sinh(77/2), (6.64) 

Pi-k — cosh(A?7/2) coshTy — sinh(A?7/2) sinh77 — COS0. (6.65) 

The integrations for the function fl^^^^ are straightforward, and we obtain 
o(«i) A ^1 ^ sinh(AV2 + Ay/2) ^ 

n\ ^ = - Ay + log (^ 3.^^^^^2 - A,/2) J • ^'-''^ 



The expression for Q^"^-* is 



nf^ = / dr) 



Ay/2 f+4>nm e"^''/^ 



Ar,/2-R i-<^nm ^TT e-''(cOSh(A77/2 - Tj) - COS 0) ' 
r+Ay/2 

= / drif{ri,Ari,R), 

JAri/2-R 

= r d777(V,A77,i?), (6.67) 

JAri/2-Ay/2 

where (pn^ is defined in the previous section, we have performed the azimuthal inte- 
gration in the second step and changed variable to rj' = Ar^/2 — 77 in the third step. 
The function / can be easily obtained, but it is rather lengthy so we do not reproduce 
it here. We now note that this expression for Q^^^^ only involves jet 1 and hence, by 
Lorentz invariance, cannot depend on the other jet and so may not be a function of 
the jet separation Arj. Therefore we write 

nf^= r dv'f{r,\R). (6.68) 

JAr]/2-Ay/2 

This function f{r}', R) has a divergence as 77' ^ 0, so we add and subtract this diver- 
gence to obtain 



= r dv f /(V, R)--\+ r % (6.69) 

JAr)/2-Ay/2 \ V / J Ar]/2-Ay/2 V 
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We can rewrite the lower limit of the first, divergence free, integral as 0, and the 
collincar divergence is now contained in the second term. Therefore we have used Ay 
as a cut-off for the divergence, and we can write 

qJ"^) = fij"^) + \og2R - log(A77 - Ay) + 0{Ari - Ay). (6.70) 

We will always denote the divergence free angular integrals, which result from such 
subtractions, as barred quantities. We can now rescale the Q^"^^ integral, using 

77 = ry'/i?, (6.71) 

to obtain 

^^i'^= fdJu-giv^R)-^). (6.72) 



We have denoted the rescaled version of / by g'. This quantity, which is only a function 
of R, can now be expressed as a power series in R and the integrals done on a term- 
by-term basis. Doing this wc obtain the rapidly converging series, 

Q = - log(2 \ + + + . . . . (6.73 

^ ^ TT 8 IStt 576 54007r 5292007r 4147200 ^ ^ 

To calculate Vi'^^^ we use the parity symmetry mapping, {a, b, 1, 2} {b, a, 2, 1}, and 
obtain the expression, 

dn / ^ . (6.74) 

7-0,^ 27re''(cosh(AV2 + r7) + cos0) 

We now perform similar manipulations to the case of Vt^^^\ However, as Og"^-* is a 
function of both final state jets, the resulting expression must be a function of A77 and 
we also note that Jlg"^-* is not divergent. We hence obtain the expression 

^t^- f dfj{R-f{fj,Arj,R)), (6.75) 

which we can expand as a power series in the variables R and z = exp(— Ar/), and 
perform the remaining integrations term-by-term. 

The pole arising in the subtraction term Q,^^^^ now cancels against an equivalent 
pole in the function Q^"^\ when we expand the latter in Ay around the point A77, 

lim nf '^^ Arj- log( At? - Ay) + log(2 sinh A77) . (6.76) 

Ay—*Ar] ■' 
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Therefore we find the final, divergence free, form of Q^" ' as 

n'j^^ = - A77 + log(2 sinh A77) - log(2i?) - fiS"^^ - Q^"^^ . (6.77) 

We have presented the full set of series expansions in appendix G and these, together 
with equations (6.59-6.61), are sufficient to compute the set of kt defined momentum 
integrals and hence the corresponding anomalous dimension matrix. It is worth noting 
that, although the off-diagonal terms for the kt anomalous dimension matrices are no 
longer pure imaginary, as in the cone case, their real parts still vanish for large A77. 
Indeed for A77 = 2, the real part is more than two orders of magnitude smaller than 
the imaginary part. We have listed the closed- form momentum integrals for the cone 
defined final state using our notation in appendix H. 

6.5 Results 

We now have the tools we need to calculate resummed cross sections at HERA, which 
correctly include primary emission to all orders and secondary emission approximately 
in the large N^. limit. The colour bases used for the contributing partonic cross sections 
are presented in the appendix, along with the decomposed hard and soft matrices. We 
also present the complete colour mixing matrices and the correct sign structure for 
the three classes of diagram. Therefore we can use the eigenvectors and eigenvalues of 
the soft anomalous dimension matrices, together with the hard and soft matrices, to 
calculate the primary resummed cross section using equation (6.46), for either a kt or 
a cone defined final state. The differential cross section, in Arj, can then be computed 
using the cuts given in section 6.2, both for the totally inclusive cross section (no 
gap) and for the gap cross section at fixed Qq. The gap fraction is then found by 
dividing the latter quantity by the former. All our results are computed using GRV-G 
(LO) photon parton densities [83] and the (HERWIG [17] default) MRST-LO proton 
parton densities [84]. We have included an estimate of the theoretical uncertainty 
in the primary resummation by varying the hard scale in the evaluation of as, while 
keeping the ratio of the hard and soft scales fixed. 
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2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 



Figure 6.4: The cross sections for the HI data (left) and the ZEUS data (right), 
which was defined using the kt algorithm with R — 1.0. On both plots the solid 
line is the total inclusive cross section, the dashed line is the gap cross section for 
= 1 GeV with only primary emission included, and the dotted lines indicate the 
range of theoretical uncertainty in the prediction. 



6.5.1 Totally inclusive ep cross section and the gap cross sec- 
tion 

The left hand side of figure 6.4 shows the totally inclusive dijet cross section for 
the HI analysis and the gap cross section for = 1-0 GeV. We have not shown 
further values of Qq as all the plots show qualitatively the same behaviour. We have 
cross-checked our total inclusive cross section against the Monte Carlo event generator 
HERWIG [16, 17] and we obtained complete agreement for the HI and both the ZEUS 
sets of cuts. In figure 6.4 the solid curve is the total inclusive cross section, the 
dashed line is the cross section with the primary inter jet logarithms resummed and the 
dotted lines show the theoretical uncertainty of the primary resummation, estimated 
by varying ag as described above. The inclusion of the primary gap logarithms gives 
a substantial suppression of the cross section; our analysis confirms simple "area of 
phase space" arguments which say that the kt defined final state will have greater 
soft gluon suppression than a cone defined final state due to the increased gap area in 
the (77, 0) plane. This plot for the ZEUS analysis is shown in the right hand side of 
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Figure 6.5: The gap fractions for the HI analysis with a kt defined final state {R = 1.0), 
at varying Q^. = 0.5, 1.0, 1.5, 2.0 GeV for plots (a), (b), (c) and (d) respectively. 
The HI data is also shown. The solid line includes the effects of primary emission 
and the secondary emission suppression factor. The overall theoretical uncertainty, 
including the primary uncertainty and the secondary uncertainty, is shown by the 
dotted lines. The dashed hne indicates the gap fraction obtained by including only 
primary emission. 
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Figure 6.6: Tfie gap fractions for the ZEUS analysis with a kt defined final state {R — 
1.0), at varying Qq. Qq = 0.5, 1.0, 1.5, 2.0 GeV for plots (a), (b), (c) and (d) re- 
spectively. The solid line includes the effects of primary emission and the secondary 
emission suppression factor. The overall theoretical uncertainty, including the primary 
uncertainty and the secondary uncertainty, is shown by the dotted lines. The dashed 
line indicates the gap fraction obtained by including only primary emission. 
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figure 6.4. 

6.5.2 Gap fractions 

The gap fraction is defined as the gap cross section, at fixed Qq, divided by the 
total inclusive cross section. Figure 6.5 shows the gap fraction for the HI cuts at the 
four experimentally measured values of Qq and figure 6.6 shows the gap fractions for 
the ZEUS analysis. The sohd hue is the gap fraction curve obtained by including 
the primary emission and the NG suppression factors of table 6.2 in the prediction 
for the gap cross section. The dotted lines show the theoretical uncertainty of both 
the primary and secondary emission probabilities, and the dashed line shows the gap 
fraction obtained by including only the primary emission contribution. We find that 
our gap fraction is consistent with the HI values for the measured Qq,. The large 
uncertainty in the gap fraction predictions comes from an approximate treatment of 
the NG suppression and from using perturbation theory at ■--'1 GeV. Nonetheless, 
this uncertainty is principally in the normalisation of the curves and we expect our 
resummation to describe accurately the shape of the gap fraction curves. 

6.6 Conclusions 

In this chapter we have computed resummed predictions for rapidity gap processes 
at HERA. We include primary logarithms using the soft gluon techniques of GSS, as 
surveyed in chapter 4, and include the effects of NGLs using an overall suppression 
factor computed from an extension of our earlier work in chapter 5. The kt definition 
of a hadronic final state determines the phase space available for soft primary emission 
and we have computed a set of anomalous dimension matrices specific to the geometry 
of the HI and ZEUS analyses. Of course this method can be used for any definition 
of the gap, provided is directed away from all hard jets. We then compared our 
predictions with gaps-between-jets data from the HI collaboration and found a consis- 
tent agreement. The theoretical uncertainty of our predictions is relatively large, and 
generally dominated by the secondary emission uncertainty. However our resummed 
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predictions correctly predict the shape of the HI data, and the normalisation agrees 
within errors. There is a suggestion that the dependence is not quite right, with 
our central Qq = 0.5 GeV prediction below data and our central = 2.0 GeV 
prediction above data, although all are within our uncertainty. It is possible that 
a more complete treatment of the perturbative/non-perturbative interface would im- 
prove this. We expect that calculation of primary emission will be correct if Qq is not 
too large, so that we can neglect jet recoil. However our calculation is of sufficient 
accuracy in the region of phase space probed at HERA. 

Our treatment of the NGLs is very approximate. For a fuller treatment, it is 
necessary to extend the extraction of the suppression factor to beyond the large Nc 
limit and overcome the inherent disadvantages of the numerical methods used. For 
the current application, consideration of the four jet system described in chapter 5 is 
also necessary. We reserve the latter extension, in the large Nc limit, for future work. 
Our calculation has not included power corrections [85]. The inclusion of such non- 
perturbative effects is required for a full and correct comparison to the experimental 
data. Again, we reserve this for future work. 

Our calculation involves a numerical integration over all kinematic variables, so 
it would be straightforward to calculate the dependence of the gap fraction on, for 
example, the fraction of the photon's momentum participating in the hard process, x^. 

In conclusion, in this chapter we have shown that the calculation of primary and 
secondary gluon emission patterns, using the tools of chapters 4 and 5, can give a 
good description of rapidity gap data at HERA. A fuller treatment would refine our 
approximation of NGLs and include power corrections. 



Chapter 7 



Conclusions 



In this thesis we have studied rapidity gap processes at modern coUiders. Using Quan- 
tum Chromodynamics, we have made detailed calculations and made comparisons to 
the newest experimental analyses; we have found that the predictions of the theory 
agree well with the available data. 

In the study of gap-jet-gap processes, we were able to show that the factorised 
model of Ingelman and Schlein is able to describe the diffractive observations at the 
Tevatron. To achieve this we used pomeron parton densities obtained from HERA 
experiments and a gap survival factor which is consistent with theoretical estimations. 
This issue of gap survival is not very well understood at the present time but it is 
encouraging that we are able to use the HERA parton densities at the Tevatron and 
have a framework in which to understand hadronic diffractive processes. 

We then turned our attention to the total gap energy flow in rapidity gap pro- 
cesses, which is an exciting way study these processes in a perturbative way. The 
contribution to these cross sections from primary emission, meaning gluons that arc 
emitted directly into the gap, can be resummed if we replace the partons in the hard 
scattering by ordered exponentials, or eikonal lines. This exponentiation is driven by 
a set of anomalous dimension matrices, which are sensitive to the geometry of the soft 
gluon emission phase space. We were led to compute a set of these matrices for the 
geometric definition of the gap cross section used in recent HI and ZEUS analyses, 
which we then used to calculate resummed primary emission gap cross sections. 
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However, the leading logarithmic behaviour in such events also includes emission 
into the gap by harder ghions, which arc emitted out of the gap. Such secondary, or 
non-global, logarithms have a significant numerical importance in rapidity gap cross 
sections. We were unable to include the effects of secondary emission in our primary 
emission calculations, and we performed numerical calculations in the large Nc limit 
to extract an overall scale and gap dependent non-global suppressive factor. 

The cross sections we computed, which included both primary and secondary 
emission, were then compared to HI gaps-between-jets data. The predictions were 
consistent with the data within errors, although the theoretical uncertainty was rel- 
atively large. The description of the shape of the data was good. The inclusion of 
power corrections, and a refined treatment of non-global logarithms, would improve 
the fit to data. This is reserved for future work. 

In summary, the research work in this thesis is based on the calculation of rapidity 
gap processes using QCD and the successful comparison to experimental data. Further 
study of these topics will shed light on the mechanisms of QCD and the interplay 
between short and long distance physics. 



Appendix A 



Useful identities from group theory 



There are many SU(3) group theoretic identities used in this thesis. In this appendix, 
based on [86], we will summarise the most frequently occuring identities. The funda- 
mental quark indentity is 

which allows the expansion of quark-gluon vertices. For gluons we have 



dabc = 2[Tr(m'')+Tr(nV)], (A.2) 
Uc = -2z[Tr(rt'=)-Tr(W)], (A.3) 

which allow the expansion of triple gluon vertices into their colour content. The first 
one is the symmetric form and the second is the antisymmetric form. Only the latter 
appears in the QCD Feynman rules. We can also expand a closed quark-gluon loop 
using 

itH% = (A.4) 
where Cp = 4/3, and multiply two SU(3) matrices using 



.a.b _ }_ 



J^SabSik + {dabc + ifabc) ^Ik 



(A.5) 



The antisymmetric structure constants are antisymmetric under the interchange of 
two indices and obey the Jacobi relation, 

fabefecd ~l~ fcbefaed ~l~ fdbeface — 0. (A. 6) 
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The summetric structure constants are symmetric under the exchange of any two 
indices, obey dabb = and satisfy 

fabedecd + fcbedaed + fdbedace — 0. (A. 7) 

These constants also satisfy 

- 4 

dacddbcd — ^ab, (^-8) 

and the antisymmetric and the symmetric structure constants are related 

i*} = + (A.9) 

The following identies are also useful, 

dijkfijk = 0, (A. 10) 

fijkfijk = NJii, (A.ll) 
iVc — 4 

dijkdijk = — — — Sii, (A- 12) 

fmigfgjnfnkm — Q^fijki (A. 13) 

dmigf gjnfnkm Q^dijkj (-'^•l^) 

- 4 

dmigdgjnfnkm — 2iV •^''•i^'' (A. 15) 

N"^ — 12 

dmigdgjndnkm — 2J^ dijk- (A. 16) 



Appendix B 

Eikonal Feynman rules 



The Feynman rules for the eikonal graphs can be extracted from the definition of the 
eikonal cross section. These rules have appeared several times in the literature and we 
list them here for completeness. A propagator for a quark, antiquark or gluon eikonal 
line is, 

ov ■ q + ie 

where 5 — +1(— 1) for the loop momentum flowing in the same (opposite) direction 
as the vector v'^. A propagator for a gluon is given by 

The interaction vertex for a virtual gluon (with colour index c) with a quark or an 
antiquark eikonal line is 

-igstty^: (B.3) 
where is an SU(3) matrix in the fundamental representation and A is +1(-1) for a 
quark (antiquark) eikonal line. The interaction vertex for a virtual gluon with a gluon 
eikonal line is 

-gsF'^'v^A, (B.4) 

where f"^'^ is a SU(3) structure constant and A is +1(-1) for a gluon located above 
(below) the eikonal line. We agree to read the gluon colour indices in an anticlockwise 
direction. This last convention is necessary because of a minus-sign ambiguity in the 
gluon/gluon vertices. 



Appendix C 

The derivation of the subtraction 
formula (5.69) 

The goal of this appendix is to derive the subtraction formula (5.69), for 5*2 with 
clustering. Therefore we wish to integrate over the phase space which is vetoed by 
the clustering procedure, where the two gluons are merged by the algorithm. The 
full result for S2 with clustering is then found by subtracting this result from the 
unclustered result. We begin by writing down the general integral, 

r = r dr^, r dr^2 f ^fiVi - V2, </>), (C.l) 

where we integrate over all possible phase space, and the function /(r^i — 1]2,4>) is 
obtained from the squared-matrix element, equation (5.23). We can use the fact that 
the integrand is only a function of rji — 772 (boost invariant) to centre the edge of the 
gap on ?7 = 0. Therefore, ensuring that gluon one is outside the gap and gluon two is 
inside the gap, we obtain the following integration region, expressed as 9-functions, 

e(77i)e(-772)e(772 + Ar))e{R'' - {rji + 772)' - 0'). (C.2) 

The final 6-function ensures we only integrate over a region where the two gluons 
would be clustered i.e. within in the (77, 0) plane. Now we make a variable change 

V = V1 + V2, (C.3) 
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which has the Jacobian 



drji dri2 = -drj drj. (C.4) 



The ©-functions become 

e»(^(^ + ^))6>(-^(^ - - r?) + Ar^)e(i?2 - r^^ - 0^), (C.5) 

which can be trivially rewritten 

6(77 > -^7)6(7; > 77)6(77 > 77 - 2Ari)e{R^ > rj^ - 0^), (C.6) 

where we have adopted a non-standard notation which is better to use in this appendix. 
We now expand the first and the third 6-functions into a sum of two terms with definite 
orderings: the first with the ordering 77 > —77 > 77 — 2A77 and the second with the 
ordering 77 > 77 — 2A77 > —77, 

6(77 > -77)6(77 > 77 - 2A77) = 6(77 > -77)6(-77 > 77 - 2A77) 

+ 6(77 > 77 - 2A77)6(r7 - 2A77 > -77). (C.7) 

We can now deduce that 

6(77 > -77)6(77 < 77) =^ 6(77 > 0), 
6(-77 > 77 - 2A77) ^ Q{7] < At]), 
6(77 > 77 - 2A77)6(77 < 77) =^ 6(A77 > 0), 

6(77 - 2A77 > -77) 6(77 > A77). (C.8) 

Therefore we find that 

1 r°° r d6 r°° 
= 2 y "^"^J £j ^^/(^»[e>(A77- 77)6(77)6(77 -77)6(77 + 77) 

+6(77 - A77)e(77 - 77)6(77 - 77 + 2A77)]6(i?2 - 77^ - (f)'^). (C.9) 

The integrand is not a function of 77 and so we can do that integral, with different 
limits for the two terms. The 6-functions for the first term give 

rdfj^2rj, (C.IO) 

J-r] 
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and the second term gives 

/ dr] = 2Ar]. (C.ll) 

Jr)-2Ar] 

The remaining ©-functions now allow us to cover the whole of the integration region 
(the first term goes from to Arj and the second term goes from Ar) to oo) and we 
get 

I" = dv f ^fiv, <P)min{rj, Arj)e{R' - rj' - 0^), (C.12) 
The expression for /(jy, 0) is derived in section 5.4.1, 

fM) = s( J°f'''' (C.13) 

\cosh(?7) — cos(0) / 

and we note it is even in 0. We arrive at 



Jo Jo 27r Vcosh(r/) - cos(0) / 

where we have changed the upper limit of the r} integration because of the 6-function 
and the requirement that be a real number. We can do the final integration to 
obtain 

o fR ( Aanf ^^ )\ \ 

7" = - y d77min(77, Ar]) I 2coth77 arctan I ^^^^D I " - I > (C.15) 

and we obtain our final result from 



SI = -ACfCaP. 



(C.16) 



Appendix D 



Colour bases 



In this section we present the colour bases used in chapter 6 of this thesis. All the 
bases in this section have appeared in [8,9,55,58]. 

The process qq qq 



C2 = --^SalSb2 + ^SabSi2- (D.l) 



The process qq qq 



Cl = SalSb2, 

C-2 = --;^5al8b2+ -j5a25bl- (D-2) 



The process qg — > qg 



Cl = 5al5b2, 

C2 = db2ciTp)ia, 

C3 = ifmdmia- (D.3) 
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The processes qq gg and gg qq 

The process gg — > qq has the basis, 

Cl = SabSl2, 

C2 — dabc(Tp)l2, 

C3 = tfabcinU- (D.4) 

To find the basis for qq ^ gg, we interchange a 2 and 6 1. 

The process gg — > gg 

The complete basis is 

{ci, C2, C3, Pi, Pss, -Psa) -Pio©To' -P27} , (D-5) 



where 



2 

Cl = - [fabcdl2c — dabcfl2c] , 

C2 = ^ [/afec'^12c + '^a6c/l2c] ) 

C3 = ^ [/alcC^62c + C^alc/62c] , 

Pi = o^al562, 
o 

3 

Pss = —dalcdb2c, 
Psa — -^falcfb2c, 

PlO0TO = - ((506512 — (^a2(^6l) — -/alc/62c) 

1 13 
-P27 = 2 ('^»&'^12 + (5a2(56l) - g(5al(562 " -^dalcdb2c- (D-6) 



The direct processes 

Since there is only one colour structure, these are basis independent. 



Appendix E 



The hard and soft matrices 



We now show the complete set of hard and soft matrices used in 6. These matrices 
have appeared in a variety of forms in [8,9,55,58]. In all these equations we have set 
Nc — S and have written the coupling scale as /j,. Note that all our hard matrices 
differ from the normalisation used in [8,58] by a factor of 2s/7r and from that used 
in [55] by a factor of Atu/s^^ while they agree with that used in [9]. 



The process qq — > qq 

The hard matrix has, in the basis D.l, the form 

16 




^(1) _ 1 «^M7r / siXi 27X2 j 
9 5 



where we define 



Xi 



X2 = 3- 



^9 ) 



St ' 

= ^ + 9^-3^ ^^-'^ 
The unequal flavour process qq' qq' is found by dropping the s-channel terms from 
these equations, and the unequal flavour process qq q'q' is found by dropping the t- 
channel terms. The hard matrix for qq — > qq is found using the transformation t u. 
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The corresponding soft matrix for all these processes is 



m-^) 



= I ° I . (E.3) 



The process qq qq 

The hard matrix has, in the basis D.2, the form 

^(1) ^ 1 [ §Xi ^X2 

9 s 4 



where we define 



Xi 



X2 = 3^ - 



27 



X2 X3 



(E.4) 



- ^ + 9^-3^^- ^^-^^ 

For the process qq' — > ^g' only keep the t-channel terms. The corresponding soft 
matrix is 

= I ° I . (E.6) 





The process qg — > qg 



The hard matrix has, in the basis D.3, the form 



where we define 



^(1) ^ 1 ^sMtT 

24 2s 



^ IsXl 6-^1 



Xi 

X2 
X3 



2 
1 

3 



1X2^ 



ixi |xi X2 

\ |x2 X2 X3 y 



St/, ' 

2 SM st 
^su 1 P 
~ 2M' 



(E.7) 



(E.8) 
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The hard matrix for the process qg — > gq is found by the transformation t u. The 
corresponding soft matrix is 



/ N^(N^-l) ^ 

2fe(^c-4)(iVe'-l) 

V 1N,{N^-1) J 



The processes qq gg and gg qq 



(E.9) 



In the basis D.4 the hard matrix for these processes has the form 



^a|(/i)7r 
A 2s 



qXi 2X1 

V ^2 1x2 ks J 



where we define 



Xi 

X2 
X3 



iu 

1 + yjx. 

4tM\ 

i--^jxi. 



(E.IO) 



(E.ll) 



The constant A = 9 for the process qq gg and A = 64 for the process gg qq. 
The matrix for the process gg — > qq is found from the transformation t u. The soft 
matrix is 



C{0) ^ 

2Nr 



2m 











\ 



N^-A 











(E.12) 



The process gg gg 



The hard matrix, in the basis D.5 has the block-diagonal form 



-'3x3 '-'3x5 
-'5x3 -"5x5 



(E.13) 
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where the matrix i^sxs has the form 



H, 



(1) _ 1 "s(^) 



TT 



5x5 



16 2s 



/ 


9X1 


ixi 


1X2 





-3xi 


\ 




fxi 


fxi 


1X2 





-fxi 








9 

1X2 


X3 





-|X2 























V 


-3xi 


-fxi 


-1X2 





Xi 


/ 



and we write 



Xi 

X2 

Xa 



= 1- 


iu 

¥~ 


si 
- — + 


su ' 


si 

V? 


tit 


u? 

+ - 

St 


52 

tu 


27 




' su 


liu 



- 1 St 



+ 



9 



u 



32 



+ 



2 V st tu 



1^ 

2 srt 



For this process the soft matrix is 



v 















000000 20 

000000 27 



5 

5 

5 

1 

8 

8 



The direct processes 



(E.14) 



(E.15) 



(E.16) 



For both these processes the zeroth order soft factor is unity and the hard functions 
are 



u-(i)/ -X fsr^ 2\«sQ;em7r UU U t\ 



2 GL qGLqt^^TX ^ify Hi 



3 ^ 2s 



-U S 

— H 

s —u 



(E.17) 
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where is the electric charge of quark flavour q, in units of the electron charge. Note 

that if the sum for 'yg — qq is taken to be over four flavours, then this gives a factor 
of 10/9. 



Appendix F 

Colour decomposition matrices 



We now give the full set of colour decomposition matrices used in chapter 6 of this 
thesis, and also the sign function S, defined by equation (6.53), for a, (5 and 7, defined 

by 

7 = 56ir(")+5„2r('^2)_ (p_;l) 



The process qq — > qq 



The signs are 



- - I CfP ^(a + 7) , ^ ^ 

a + j Cpa- 2^{a + (3 + 2j) 



Sa = +1, (F.3) 
Sp = +1, (F.4) 
= -1. (F.5) 
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The process qq qq 



The signs are 



a + 7 Cf7- 2^(2q; + /3 + 7) 



Sa - -1, (F.7) 
.5^ = +1, (F.8) 
= +1. (F.9) 



The process qg qg 



-"ig^qg 



\ 





-(a + 7) 



jVg-4 



-kia + j) \ 

(a + 7) X J 



The signs are 



(F.IO) 



and we define 



Sa = +1, 

Sp = +1, 
— — 1, 



X 



(F.ll) 
(F.12) 
(F.13) 

(F.14) 



The processes qq gg and gg qq 



For qq — > gg we have 





(/3 + 7) 







f (/9 + 7) 



iV|-4 
4iVc 



(/? + 7) xl ) 



:f.i5) 
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The signs are 



Sa = +1, 
Sp = +1, 
— — 1, 



and we define 



X 



For gg — > qq we have 

^ Cppi^ + CAr(«^) i(/3 + 7) ^ 

V (/5 + 7) 



f(/? + 7) 



The signs are 



and we define 



AT, 



5a = +1, 

Sp = +1, 
— —1, 



(F.16) 

(F.17) 
(F.18) 

(F.19) 



(F.20) 



(F.21) 

(F.22) 
(F.23) 

(F.24) 



The process gg gg 



Q9g-*99 — 



M3X3 0,3x5 



05x3 A^5x5 



where the matrix Alsxs is 



M3X3 = 



/ f{a + P) 







^(Qi-7) 

f(/3-7), y 



(F.25) 



(F.26) 



Appendix F. Colour decomposition matrices 



160 



and the matrix Alsxs is 
/ 



M5X5 = 



3/3 




(a 



7) 








K« + 2/3-7) 



for Nc — 3. The signs are 



(a 



7) 



a + 7) 




S{a + 7) 

i(« + 7) 

(a + 2/3 - 7) 



|(a + 7) 



i(« + 7) 





>-7) 






i(« + 7) 

nT(" + 7) 



|(a + 7) 2a-/3-27 y 

(F.27) 



»5a 



+ 1, 
+ 1, 
-1. 



(F.28) 
(F.29) 
(F.30) 



The direct processes 



This processes has no matrix structure, and so we present colour decomposition func- 
tions. 

(^jq-^gq ^ 3_r(''2) + ^ (r^^^) + r^^^)) . (F-3i) 



Appendix G 

The r^^-^) series expansions 



We have not found a closed form for these integrals, but they are straightforward to 
express as power series in R and e~^^ (by Lorentz invariance, only the contributions 
from dipoles containing jet 2 are Ajy-dependent) , 



(al) ^ D , 1 1 1 



At] — Ay 

0.31831 R + 0.06250 i?^ - 0.00884 R^ + 0.00087 i?^ 



-0.00003 (G.2) 
^(61) ^ ^(".liog^. liog^ — 0.31831 i?- 0.06250 i?^ (G.3) 

- 0.00884 R^ - 0.00087 R^ - 0.00003 i?^) , (G.4) 
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vr V2 




Arj — Ay 



1 



+ 



(G.5) 



(+0.31831 R + 0.06250 + 0.00884 + 0.00087 i?^ + 0.00003 R^) 
+(-0.08616 - 0.03383 - 0.01197 i?^ - 0.00282 R^ 
-0.00039 i?^ + 0.00001 i?^)e-(^^-^) + 
(+0.01166 R + 0.00916 R^ + 0.00551 R^ + 0.00305 i?^ 
+0.00122 R^ + 0.00038 R^ + 0.00011 R'^ + 0.00003 R^)e-^^^'^-^^ 
+(-0.00158 R - 0.00186 R^ - 0.00162 R^ - 0.00139 i?^ 
-0.00088 i?^ - 0.00041 i?^ - 0.00017 i?^ - 0.00007 i?^ 
-0.00002 i?9)e-^(^^-^) + 

(+0.00021 R + 0.00034 R^ + 0.00039 R^ + 0.00045 i?^ 

+0.00039 R^ + 0.00024 R^ + 0.00013 R'^ + 0.00007 i?^ 

+0.00003 R^ + 0.00001 i?^°)e-^(^^-^) + 

(-0.00003 R - 0.00006 i?^ - 0.00008 i?^ - 0.00012 i?^ 

-0.00013 i?^ - 0.00010 i?^ - 0.00007 i?^ - 0.00004 i?^ 

-0.00003 i?^ - 0.00001 i?^°)e-^(^^-^) + 

(+0.00001 R"^ + 0.00002 R^ + 0.00003 R'^ + 0.00004 i?^ 

+0.00003 R^ + 0.00003 i?^ + 0.00002 R^ + 0.00001 R^)e-^^^''-^A ,(G.6) 



Appendix G. The r(^^) series expansions 



163 



Q{b2) ^ (+0.00458 + 0.00389 + 0.00229 + 0.00104 



+0.00038 + 0.00012 + 0.00003 R^)e-^^^'T^^ + 
(-0.00124 i?^ - 0.00158 i?^ - 0.00124 i?^ - 0.00079 i?^ 
-0.00041 i?^ - 0.00018 i?^ - 0.00007 i?^ - 0.00002 R^)e-^^^''~^^ 
+ (+0.00025 R^ + 0.00043 R^ + 0.00042 i?^ + 0.00034 R^ 
+0.00024 R^ + 0.00014 R^ + 0.00007 R^ + 0.00003 R^ 
+0.00001 i?^°)e-^(^''-') + 

(-0.00005 R^ - 0.00010 R^ - 0.00011 i?^ - 0.00011 R^ 
-0.00010 i?^ - 0.00007 i?^ - 0.00004 i?^ - 0.00002 R^ 
-0.00001 i?^°)e-^(^''-2) + 

(+0.00002 R^ + 0.00003 i?^ + 0.00003 R^ + 0.00003 R^ 



71 




(G.7) 
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nif) ^ ^[-\r'+ (g.8) 

(+0.06767 + 0.02872 i?''^ - 0.00096 i?^ + 0.00002 i?^)e-(^''-2) 
+(-0.01374 - 0.01166 - 0.00229 i?^ - 0.00038 i?^ 
-0.00021 i?^ - 0.00003 i?^)e-2(^^-^) + 
(+0.00248 + 0.00316 + 0.00124 i?^ + 0.00032 R^ 
+0.00041 R^ + 0.00027 R^ + 0.00007 R^ + 0.00001 R^)e-^^^''-'^^ + 
(-0.00042 R^ - 0.00071 R^ - 0.00042 i?^ - 0.00019 R^ 
-0.00024 i?^ - 0.00019 i?^ - 0.00007 i?^ - 0.00002 R^ 
-0.00001 i?^'^)e-^(^^-') + 

(+0.00007 + 0.00014 i?''^ + 0.00011 i?^ + 0.00007 i?^ 
+0.00010 R^ + 0.00009 R'^ + 0.00004 R^ + 0.00002 R^ 
+0.00001 i?^°)e-^(^''-2) + 

(-0.00001 R^ - 0.00003 R^ - 0.00003 i?^ - 0.00002 R^ - 0.00003 i?^ 
-0.00004 i?^ - 0.00002 R^ - 0.00001 i^^je-^^^""^)) , (G.9) 

where all coefficients larger than 10~^ are shown (recall that we are mainly interested 
in the case R — 1, Arj > 2). By symmetry, we have Qg''^ = where the mapping 
i — > z is given by {a, b, 1, 2} — > {b, a, 2, 1}. 



Appendix H 

The ^^j*^^ angular integrals for a 
cone geometry 



We present these results as they have not appeared previously in this form. They have 
the expression, 

' ~ J -Ay/2 Vo 2^{P.-kKPrky ^""-'^ 

where the integrand is found from the appropriate 4-momenta, and the phase space 
is taken to be of width Ay and azimuthally symmetric. Note that these expressions 
do not include the sign factors. We obtain 

of) = 2Ay, 



f 

(12) _ / sinh(AV2 + A|//2) \ 

^ Vsmh(Ar//2 - Ay/2) J ' 
o(ai) . , ^_^ sinh(AV2 + Ay/2) >^ 

= -^^ + ^"g^ sinh(AV2-Ay/2) j' 
sinh(A77/2 + A?//2)\ 
sinh(Ar//2 - Ay/2) ) ' 
n^^'^ - Ay I i,,/ «i"MA77/2 + Ay/2) \ 



of) = -Ay + log(^ 



- A. + logf -^t^^f + y,ol V (H-2) 
^ ^ Vsinh(A77/2 - Ay/2)y ^ ^ 



-I /-> f 



Appendix I 

The ZEUS energy flow 
measurements 

In this appendix, wc present our theoretical gap fraction predictions for the ZEUS 
analysis, figure 6.6, together with a comparison to the (brand new) preliminary ZEUS 
data points. This comparison was not possible in chapter 6 due to the unavailability 
of data when this thesis was written, and we perform the comparison in this extra 
appendix. Figure I.l shows the gap fraction curves and the ZEUS data, together with 
the theoretical uncertainty. This uncertainty, in agreement with figure 6.6, includes 
both the primary and secondary emission uncertainty. 

We find that our gap fraction predictions are consistent with the ZEUS values for 
the measured E'ff'^ although, as for the HI data, there is a large theoretical uncertainty. 
However, this uncertainty is principally in the normalisation of the curves and we can 
see that our resummation describes accurately the shape of the gap fraction curves. 



inn 
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0.25 
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Figure I.l: The gap fractions for the ZEUS analysis with a kt defined final state, 
at varying E^'^{= Qn). The (preliminary) ZEUS data is also shown. The solid 
line includes the effects of primary emission and the secondary emission suppression 
factor. The overall theoretical uncertainty, including the primary uncertainty and the 
secondary uncertainty, is shown by the shaded region. 
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